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A  Lagrangian  analysis  of  the  turbulent  dispersion 
process  was  examined  by  observing  the  motion  of  suspended 
particles  in  a  turbulent  flowing  liquid.   The  liquid,  tri- 
chl oroethy 1 ene  ,  was  pressure-driven  through  a  0.0254  meter 
(1  inch)  I.D.  glass  pipe.   The  average  velocity  was  0.969 
meter/sec  (3.18  feet/sec)  which  corresponded  to  a  Reynolds' 
number  of  110,000.   A  photographing  technique  was  developed 
which  enabled  the  photographing  of  particle  motion  over 
large  distances.   Particle  positions  in  time  were  computed 
from  the  photographic  images, 

Lagrangian  time  correlations  were  computed  from  the 
displacement  and  velocity  deviations  obtained  from  the 
particle  trajectories.   The  time  dependence  of  the  various 
correlations  was  obvious  from  an  examination  of  the  plotted 
results.   The  effect  of  a  confining  wall  on  the  radial  spread 


of  the  particles  was  apparent  for  large  times.   It  appeared 
that  the  assumption  of  a  constant  radial  eddy  diffusivity 
coefficient  could  be  used  as  a  fair  approximation  of  the 
time-dependent  coefficient.   Within  the  experimental  accur- 
acy of  the  data,  the  resultant  correlations  agreed  well  with 
the  theory  and  published  data. 
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CHAPTER  1 
INTRODUCTIOr 


Turbulence  is  a  common  occurrence  in  both  natu- 
rally and  artificially  generated  flow  systems.   The  rapid 
spread  of  smoke  from  chimneys  and  the  vapor  trails  of 
high-altitude  jet  aircraft  are  examples  of  turbulent 
effects  in  the  atmosphere.   The  diffusion  of  mass  and 
heat  in  turbulent  fluids  in  numerous  chemical  engineering 
processes  is  another  example.   Few  physical  phenomena 
have  attracted  more  interest  in  the  scientific  and  en- 
gineering fields  than  turbulence,  which  attests  to  its 
general  occurrence  and  importance.   A  great  deal  of  work 
has  yielded  only  a  little  understanding  of  turbulent  fluid 
motion,  indicating  its  difficult  and  complex  nature. 

The  Lagrangian  characteristic  of  turbulent  liquid 
flow  is  the  object  of  study  in  the  present  work.   To  in- 
vestigate the  turbulence,  paths  of  small  particles  were 
observed  as  approximations  of  the  paths  of  fluid  points. 
A  distinction  is  made  here  between  a  "particle,"  a  solid, 
rigid  piece  of  matter  immersed  in  the  continuum,  and  a 
"fluid  point,"  a  mathematical  point  moving  with  the  con- 
tinuum.  A  fluid  point  is  a  volume  of  fluid  so  small  that 
in  the  context  of  the  continuum  it  may  be  considered  as 
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an  individual  point  moving  along  with  the  fluid.   The 
system  selected  for  study  is  fully  developed  pressure- 
driven  flow  in  a  pipe.   Under  these  conditions  the  fluid 
is  in  shear  flow,  and  the  turbulence  is  nonhomogeneous 
and  anisotropic. 


CHAPTER  2 
LITERATURE  REVIEW 


The  concept  of  turbulence  is  widely  discussed, 
and  in  a  general  sense,  its  effects  and  meaning  are 
understood.   Turbulence  was  discussed  as  early  as  the 
fifteenth  century  in  the  writings  of  Leonardo  da  Vinci; 
however,  it  wasn't  until  the  nineteenth  century  when 
quantitative  data  were  reported.   The  first  reportings 
on  turbulence  were  attempts  to  present  the  characteris- 
tics.  Hagen  and  Poiseuille  independently  and  qualita- 
tively observed  two  modes  of  fluid  flow  in  pipes,  non- 
turbulent  and  turbulent  (Rouse  and  Ince,  1963).   The 
concept  of  an  eddy  diffusivity  was  introduced  by  Bous- 
sinesq  (1877)  to  allow  for  turbulent  effects  analogous 
to  Newton's  law  of  fluid  viscosity.   With  the  aid  of  the 
Reynolds  number  advanced  by  Reynolds  (1883)  through  an 
experimentally  determined  d imensi onl ess  ratio  of  fluid 
and  system  parameters,  it  was  possible  to  tell  when  to 
expect  turbulence  in  a  flow  system.   Turbulence  is  said 
to  exist  in  flow  systems  for  which  the  Reynolds  number 
has  exceeded  a  critical  value.   Reynolds  (1895)  also 
recognized  the  statistical  nature  of  turbulence.   By 


expressing  the  dynamic  variables  of  the  equation  of  motion 
as  a  mean  value  plus  a  deviation,  and  by  averaging  the 
terms  formed  in  the  equation  of  motion,  the  turbulence 
term 


J\y    =  -p  u!  u'- 


(2.1) 


was  generated.   This  term  can  be  interpreted  as  the  turbu- 
lent stresses  on  an  element  of  the  fluid  in  addition  to 
the  stresses  of  pressure  and  viscous  stress.   Because 
Reynolds  was  the  first  to  give  the  equation  for  turbulent 
flow,  the  turbulence  stresses,  T--  ,  are  often  called 
Reynolds '  stresses . 

It  was  in  the  early  twentieth  century  that  useful 
and  workable  theories  and  models  of  turbulence  appeared. 
It  is  obvious  that,  in  the  early  days  of  turbulence  in- 
vestigations, when  knowledge  about  and  insight  into  the 
mechanism  of  turbulent  fluid  flow  were  rather  poor,  the 
turbulent  processes  could  be  studied  only  in  a  rather 
rough  way.   To  describe  the  behavior  of  turbulent  motion 
of  a  fluid  near  a  wall,  Prandtl  (1904)  proposed  the 
boundary  layer  theory.   The  flow  far  from  the  wall  is 
assumed  uniform.   Also  among  the  theories  arising  from 
the  early  studies  of  turbulent  flow  are  those  based  upon 
the  concept  of  a  "mixing  length."   These  theories  have 


been  the  most  fruitful,  not  so  much  because  they  describe 
the  mechanism  correctly,  but  because  they  have  resulted 
in  useful  and  practical  semi empi r ical  relationships.   Taylor 
(1915),  Schmidt  (1917),  Prandtl  (1925),  and  Von  Karman  (1934 
proposed  mixing  length  theories  for  turbulence.   The  dif- 
fusive action  of  turbulence  is  considered  to  result  in  an 
eddy  viscosity  or  eddy  heat  conductivity  from  which  the 
distribution  of  mean  values  can  be  calculated,  just  as,  in 
the  kinetic  theory  of  gases,  mol ecu  1 ar- transport  processes 
result  in  a  viscosity  and  heat  conductivity.   Basically, 
it  is  assumed  in  the  mixing  length  theories,  that  each  lump 
of  fluid  that  is  subjected  to  the  turbulent  motion  may  be 
considered  as  an  individual  entity.   The  properties  of  the 
lump  are  conserved  during  a  certain  time,  i.e.  over  a  cer- 
tain distance.   At  the  end  of  this  time,  the  lump  is  as- 
sumed to  mix  with  the  surrounding  fluid.   These  theories 
are  purely  phenomenol ogi cal  .   Later  studies  have  shown  that 
the  physical  picture  based  upon  the  concept  of  a  mixing 
length  cannot  be  correct  in  all  details.   The  mixing  length 
theories  still  prove  to  be  very  useful  to  engineers. 

The  statistical  theory  of  turbulence  also  began 
its  existence  in  the  early  twentieth  century.   The  concepts 
of  the  correlation  coefficient,  spectrum  function,  and 
local  similarity  which  have  proved  so  helpful  in  revealing 
the  structure  of  turbulence  are  all  a  part  of  this  theory. 


Taylor  (1921)  advanced  the  concept  of  the  Lagrangian 
correlation  coefficient  which  provided  a  theoretical 
basis  for  turbulent  diffusion.   Richardson  (1920,  1926) 
postulated  the  turbulence  consists  of  a  hierarchy  of 
eddies.   The  energy  is  cascaded  from  the  largest  eddies 
to  smaller  and  smaller  eddies  until  the  molecular  level 
is  reached  at  which  point  the  energy  is  dissipated  by 
viscous  effects.   Taylor  (1 935a  ,b  ,c  ,d  )  again  considered 
the  statistical  nature  of  turbulence  in  a  revised  form 
and  defined  many  of  the  properties  of  isotropic  turbu- 
lence.  This  was  followed  by  the  Von  Karman  and  Howarth 
(1938)  demonstration  of  the  tensor  properties  of  the 
Eulerian  correlation  coefficient  and  their  formulation 
of  the  relationship  between  the  double  and  triple  coef- 
ficients through  the  Navi er-Stokes  equations.   At  about 
the  same  time,  Taylor  (1938)  showed  the  transform  relation- 
ship between  the  correlation  coefficient  and  the  energy 
spectrum  function.   The  Taylor  concept  of  the  energy  spec- 
trum has  dominated  turbulence  research  to  the  present.   The 
first  theory  for  the  small  scale  structure  of  turbulence, 
presented  by  Kolmogorov  (1941)  and  Obukhov  (1941),  served 
as  the  starting  point  of  many  investigations  and  led  to 
general  theories  of  spectral  similarity  such  as  those  by 
Von  Karman  and  Lin  (1949)  and  Lin  (1948).   This  theory 
predicts  for  any  highly  turbulent  motion,  the  small  scale 


structure  is  isotropic  which  was  empirically  shown 
earlier  by  Richardson  (1926). 

The  following  publications,  among  others,  present 
turbulence  in  general  or  some  of  its  particular  aspects: 
Prandtl  and  Tietjens  (1934a, b),  Goldstein  (1938),  Agostini 
and  Bass  (1950),  Hopf  (1952,  1957),  Batchelor  (1953), 
Frenkiel  (1953),  Sutton  (1953),  Batchelor  and  Townsend 
(1956),  Corcoran  et  al .  (1956),  Opfell  and  Sage  (1956), 
Townsend  (1956),  Hinze  (1959),  Lin  (1959),  Priestly  (1959), 
Schlichting  (1960),  Lin  (1961),  Schubauer  and  Tchen  (1961), 
Tatarski  (1961),  Favre  (1962),  Pasquill  (1962),  Lumley  and 
Panofski  (1964),  Brodkey  (1967),  and  Monin  and  Yaglom  (1971 
Papers  considered  classics  on  turbulence  are  collected  by 
Goering  (1958)  and  Friedlander  and  Topper   (1961).   The 
textbooks  by  Bird  et  al .  (1960)  and  Brodkey  (1967)  also 
present  sections  on  turbulence. 


2 . 1   Dispersion  of  a  Scalar  Property  in  Turbulent  Flow 

Since  turbulent  fluid  motion  has  an  intrinsic  ran- 
dom fluid  motion,  an  inherent  characteristic  of  turbulent 
flow  is  its  dispersive  action;  accordingly,  the  dispersion 
process  is  more  active  in  turbulent  than  in  nonturbulent 
flow.   It  is  not  unusual  then  that  most  of  the  work  on 
turbulence  centers  around  its  transport  or  dispersive 
characteristic.   One  of  the  most  important  analytical 


inroads  into  the  problems  of  turbulent  dispersion  was 
made  by  Taylor  (1921).   Taylor's  work  not  only  laid  the 
groundwork  for  the  Lagrangian  study  of  turbulent  disper- 
sion but  it  also  represented  an  initial  step  in  the  ideas 
essential  to  the  development  of  a  general  statistical 
theory  of  turbulence.   From  this  base,  Richardson  (1926), 
Kampd  de  Fe"riet  (1939),  Batchelor  (1949,  1952,  1957),  Tay- 
lor (1954),  and  Monin  and  Yaglom  (1971)  have  successively 
expanded  and  deepened  the  theoretical  foundation  for  under- 
standing the  statistical  nature  of  turbulence  as  well  as 
its  transport  processes. 

2.1.1   Isotropic  homogeneous  flow 

Consider  a  condition  in  which  the  turbulence  in  a 
fluid  with  no  mean  flow  is  uniformly  distributed  so  that 
the  average  conditions  of  every    point  in  the  fluid  are  the 
same.   The  above  fluid  flow  condition  is  that  investigated 
by  Taylor  (1921).   Within  this  flow  field,  Taylor  analyzed 
the  dispersion  of  a  scalar  property  which  may  be  associated 
with  all  or  a  fraction  of  the  fluid  particles.   Because  of 
the  isotropy  of  the  problem,  only  one  direction  of  disper- 
sion need  be  considered,  say  the  x  direction.   By  defini- 
tion, the  velocity  of  a  fluid  particle  at  time  t,  V(t),  is 


V(t)  =  — x(t 
dt 


(2.2) 


and  its  position  at  time  t,  x(t 


t)  =  a  +  ^tV(t1)dt1 


(2.3) 


there  a  =  position  of  dispersing  particle  at  t  -  0 


Let  X'(t)  be  the  displacement  of  the  particle  from  a,  then 


X'  (t)  =  x(t)  -  a. 


(2.4) 


Taylor  further  supposes  that  the  statistical  properties 
of  V(t)  are  known,  in  particular  the  standard  deviation 

o    i 

of  V(t),  <V>,   and  the  Langrangian  correlation  coeffi- 
cient, R^(t).   R^L^(tn)  is  the  correlation  coefficient 
between  the  value  of  V  for  a  particle  at  any  instant,  and 
the  value  of  V  for  the  same  particle  after  an  interval  of 
time  t-,  ,  i.e. 


<V(t)V(t  +  t] )> 

RV"Ul)  =  ^i^72<v2(t  +  t1)>171 


(L) 


(2.5 


It  is  also  assumed  the  standard  deviation  is  constant  with 
time,  so  that 


:V2(t)>   =  <VZ(t  +  t1)>   =  <V  > 


(2.6 


]The  <•>  indicates  a  statistical  ensemble  average 


10 


then 


<L>u 


<V(t)V(t  +  t,)> 


2.7) 


One  measure  of  the  degree  of  dispersion  is  given 

2 
by  the  variance  of  the  displacement,  <X '  (t)>.   The  rate 

of  change  in  the  variance  can  be  expressed  in  terms  of 

<V2>  and  R^1-)  (t]  )  as  follows: 


First  ~   <X|2(t)>  -  2<X'(t)~X'(t)> 


2<X  '(t)V(t)>  (2.8) 


and  si  nee   X ' ( t ) 


=  /  V(t  +  t]  )dt] 


(2.9) 


then  dT<x'2(t)>  =  2</  V(t  +  t1)v(t)dt1> 


2  /<V(t  +  t, )V(t)>dt, 
-t        '         ' 


(2.10) 
(2.11) 

=  2<V2>/  R^UJdt,  .  (2.12) 

-t       '    ' 

2     ( L ) 
For  constant  <V  >,  Rv   (t-,)  is  symmetric  about  zero,  and 

therefore 


^<X'2(t)>  =  2<V2>  /  R(L)(t1)dt].  (2.13) 
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Integrating  Equation  (2.13)  yields 


<X,2(t)> 


9      t    t2    ( \  \ 
2<MC>  I   I     Rv  ;(t] )dt,dt2 
oo 


2.14) 


Equation  (2.14)  reduces  the  problem  of  dispersion,  in  a 
simplified  type  of  turbulent  motion,  to  the  consideration 
of  a  single  quantity,  the  Langrangian  correlation  coeffi- 
cient, R(L)(t). 

As  mentioned  above,  an  indication  of  the  rate  of 
dispersion  is  the  rate  of  change  in  the  variance.   Taylor 
(1935a)  explicitly  recognized  <V(t)X'(t)>  as  the  eddy  dif- 
fusivity,  E(t),  a  turbulent  transport  property  similar  to 
the  molecular  diffusivity.   The  eddy  dif fusivity  is  also 
related  to  the  Lagrangian  correlation  coefficient.   From 
Equations  (2,8)  and  (2.13) 


E(t)  =  2<V2>  /  R(L)(t,)dt,. 

n         I     I 


(2.15) 


Now  consider  the  physical  implications  of  (2.14) 
and  (2.15).   When  t  is  so  small  that  R^L'(t)  does  not 
differ  appreciably  from  1,  then 


<X,2(t)>  =  <V2>t2 


(2.16) 


and 


E(t)  -  2<VSt 


(2.17) 


i;:; 


Looking  further,  as  t  becomes  large 


R(L)(t)  ->  0  , 


<X,£(t)>  a  t 


(2.18) 
(2.19) 


and 


E(  t)  ->  constant . 


2.20) 


Defining  a  Lagrangian  integral  time  scale  as 


T,  =  f]   R(L)(t,)dt1 
L     o  I' 


(2.21) 


where  R^(t)  =  0  for  t  >  T-j  ,  the  eddy  diffusivity  for  large 
time  can  be  expressed  as 


E   =  <VC>T, 


(2.22) 


T.  is  considered  a  measure  of  the  average  longest  time  dur- 
ing which  a  particle  persists  in  a  motion  in  a  given  direc- 
tion.  Also  a  Lagrangian  integral  length  scale  may  be 
defined  as 


AL  -  <V2>1'2TL 


which  leads  to 

E 


<V2>1/2A 


L- 


(2.23) 


(2.24) 
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This  result  is  analogous  to  that  obtained  from  the  kinetic 
theory  of  gases  for  the  molecular  diffusivity.   It  is  also 
consistent  with  Prandtl's  mixing-length  theory  (Prandtl, 
1925);  however,  it  is  obtained  without  having  to  depend  on 
the  idea  of  mixture  by  subdivision  and  ultimate  molecular 
diffusion. 

A  Lagrangian  micro  time  scale  may  also  be  defined 
(Taylor,  1935a)  as 


x2 


1  (TD(L) 


2  dt2 


P/U'(t 


2.25) 


0 


Physically,  this  time  scale  is  taken  as  a  measure  of  the 
time  duration  of  the  smallest  eddies  which  are  responsible 
for  energy  dissipation. 


2.1.2   Anisotropic  homogeneous  flow 

2.1.2,1   Single- particle  dispersion 

Batchelor  (1949)  investigated  the  three-dimensional 
homogeneous  turbulent  flow  field.   Even  though  Batchelor 
started  with  a  generalization  of  Taylor's  (1921)  single- 
particle  dispersion  which  was  a  Lagrangian  analysis, 
Batchelor's  final  results  were  Eulerian.   The  Lagrangian 
to  Eulerian  transformation  was  through  probability  argu- 
ments and  some  simplifying  assumptions  about  the  statistical 
nature  of  homogeneous  turbulence. 


1  n 


To  generalize  Taylor's  analysis  to  three  dimensions, 
let  X!(t)  be  the  i-th  component  of  the  vector  displacement 
of  a  fluid  particle.   The  time  rate  of  change  of  the  covari- 
ance  of  two  displacement  components  X!(t)  and  X'.(t)  is 


written  as 


^<Xj(t)Xj(t)>  =  <X!(t)Vj(t)>  +  <V!(t)Xj(t)> 


(2.26) 


or 


^<X!(t)Xj(t)>  =  /<V!(t1)Vj(t)>dt1+  /<V!(t)Vj(t1)>dt1      (2.27) 

Using  tensor  notation,  define  the  second  order  tensor  S^  as 
f ol 1 ows  : 


s.^t-t^  =  <v:(t1)vj(t)> 


2.28) 


If  the  process  is  assumed  to  be  a  time-stationary  process 
then  the  elements  of  S  depend  only  on  the  elapsed  time 
£  =  (t  -  t,).   Now  Equation  (2.27)  can  be  written 


^<X!(t)Xj(t)>  -     I   (S1d(0  +  Sj.(c))d5 


2.29) 


and  when  integrated  yields 


t  t. 


<X!(t)X;(t)>  =  /  /  '(S  .At,)    +  S   (C))dCdt1   (2.30) 
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As  shown  by  Kampe"  de  Fe>iet  (1939)  in  a  study  of  stationary 
random  functions,  a  change  of  variable  in  Equation  (2.30) 
and  a  partial  integration  can  be  made  to  give 


<X!(t)X'.(t)>  =  /  (t  -  ?  )(S.  .(O  +  S.,  (?))(£.(  2. 31) 


i     J 


lj 


Jl 


The  above  results  can  be  expressed  in  terms  of  La- 
grangian  velocity  correlation  coefficients.   Defining  a 
generalized  coefficient  as 


R<L>(t  -  t, 


<v;(t1)v'.(t)> 


<V!2(t1)>1/2<V'.2(t)>1/2 


S.  .(t  -  t,) 

_LL      ' 


<V!2(t1)>1/2<Vj2(t)>1/2 


(2.32) 


(2.33) 


and  since  the  turbulence  is  homogeneous  and  stationary 


<V.2(t)>  =  <V  !2(t,  )>  -  <V!2> 


(2.34) 


then  Equations  (2.29)  and  (2.31)  become 


-d_<x:(t)x'.(t)>  =  <v:2>1/2<v'.2>1/2  w\\]U)  +  R(,-)(d)dc        (2.35) 

dt       i         j  i  J  „    u  J' 


and 


cx:(t)X'.(t)>  =  <v.2>1/2<v.2>1/2   fit  -  0(k\\](0  +  R^tO)^.     (2.36) 
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The  probability  distribution  of  X..(t)  is  intimately 
connected  with  the  spatial  di  stri  buti  on  of  themean  concentra- 
tion of  the  property  transported  by  the  fluid  particle. 
From  this  viewpoint  Batchelor  considers  a  volume  of  marked 
fluid  particles  with  a  probability  that  a  fluid  point  de- 
fined by  the  position  vector  X_  lies  within  the  marked  fluid 
at  time  t.   Note  that  P(X_,0)  is  unity  or  zero  depending  on 
whether  the  point  X  lies  within  or  outside  the  initial  volume 
of  marked  fluid. 

Now  by  using  the  probability  associated  with  the 
dispersion  of  a  fluid  particle,  the  probability  density 
function,  P(X,t),  can  be  related  to  its  initial  p.d.f., 
P(X,0).   Each  particle  diffuses  independently  of  its  neigh- 
bors.  Representing  by  Q(X_,t)  the  p.d.f.  of  the  vector  dis- 
placement X  of  any  fluid  particle  during  a  time  t,  then 
Q(X,t)dv(X)  is  the  probability  that  a  fluid  particle  lies 
within  a  volume  element  dv(X)  at  time  t.   Furthermore,  the 
p.d.f.,  P(X,t),  for  a  finite  volume  of  marked  fluid  is 


P(X,t)  =  /P(X',0)Q(X-X',t)dv(X'). 


(2.37) 


If  the  form  of  Q(X,t)  is  known,  the  solution  for 
P(X,t)  can  be  obtained  from  (2.37).   The  experimental  re- 
sults of  Schubauer  (1935),  Col  1 i  s  (1948),  and  Kalinske 
and  Pien  (1944)  indicate  that  the  separate  p.d.f. 's  of  the 
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X '.  s  are  normal  for  all  values  of  t.   The  above  data  are 
i 

not  sufficient  to  define  the  p.d.f.  Q  ( X_ ,  t )  in  all  its 
generality.   However,  in  order  to  make  the  analysis  as 
general  as  possible,  Batchelor  assumed  the  p.d.f. 's  of 
the  X '.  s  to  be  jointly  as  well  as  separately  normal.   In 
this  case  it  is  known  that  the  appropriate  form  of  Q(X_,t) 
i  s 


Q(X,t) 


1/3 


co   X  X 
rs  r  s 


(8t^)17 


jjj     e  x  p ( -   2  ft 


(2.38 


where 


<x^>   <x1 x2><x1 x3> 


<X?X-,>  <x2>  <x2Xo> 


<x3x]>  <x3x2>  <x^> 


(2.39 


oo  is  the  co-factor  of  the  typical  element  Q  of  the 
determinant  Q,  and  repeated  suffixes  implies  summation 
over  values  1,  2,  and  3.   Thus  Equation  (2.37)  becomes 


"&t)  -  ■£$%  f  exP(-  ^rr(v;n^;))Mr)  ■ 

[on  U) 


(2.40 


Furthermore,  Equation  (2.40)  is  a  solution  to  the  diffusion 
equati  on  , 

,2, 


P(X,t) 


P(X,t) 
ij    8X.9X.   ' 


(2.41 


provided  the  diffusion  coefficient  tensor  E_  is  such  as  to 
satisfy  the  equations 


ij  U    2  dt 


E  .  .  w  .  ui . 
ij  ir  js 


!  o  d(w  /n) 

—  0  — 

2  "     dt 


(2.42 
(2.43 


Remembering  that  by  definition 


uir  jr   "°ij 


where  6 •  ■  is  the  delta  function,  then  Equation  (2.43)  can 
be  reduced  to 


t  dn. .   1   , 

E.  ■  =  I    -rl1    =    Jp  ~  <X.(t)X.(t)> 

ij    2   dt    2  dt    i     J 


(2.44 


which  is  also  consistent  with  (2.42). 

Thus  the  diffusion  Equation  (2.41)  provides  a 
description  of  the  turbulent  dispersion  process  when 
Eulerian  analysis  is  used,  provided  the  p.d.f.'s  of  X,, 
Xp,  and  X~  are  jointly  as  well  as  separately  normal.   At 
present,  sufficient  information  is  not  available  to  sup- 
port the  above  normality  hypothesis  rigorously;  however, 
it  is  a  fortunate  discovery  that  a  single  differential 
equation  describes  the  homogeneous  turbulent  dispersion 
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process  with  high  accuracy.   Even  though  this  hypothesis 
is  presently  (Monin  and  Yaglom,  1971)  thought  to  be  exact 
or  nearly  exact,  Equation  (2.44)  should  be  treated  as  a 
phenomenol ogi cal  and  not  a  fundamental  law  of  turbulent 
dispersion. 

It  is  important  to  notice  that  the  tensor  coeffi- 
cient E..  is  a  function  of  time,  which  is  in  contrast  to 
the  molecular  diffusion  coefficient,  and  that  Equation 
(2.44)  is  valid  at  all  times.   The  tensor  E  is  a  generali- 
zation of  the  eddy  diffusivity  coefficient.   From  Equations 
(2.29)  and  (2.35)  the  generalized  eddy  diffusivity  coeffi- 
cient can  be  expressed  as 


E..(t)  ■  ^(S1  (O  +  Sj1(5))d5 


(2.45) 


or  in  terms  of  the  Lagrangian  velocity  correlation  coeffi 
c  i  e  n  t 


E.  .(t 


i<v:2>1/2<v.2>1/2 

2    l       3 


^(R^U)  +  R^U))^ 


(2.46) 


In  addition,  a  generalized  Lagrangian  integral  time  scale 
may  be  defined  as 


tu  -rO +  "$y<«»d«- 


(2.47) 


2  0 


In  view  of  (2.46),  for  large  times  the  generalized  eddy 
diffusivity  coefficient  can  be  expressed 


(E  )..  ■  '<v:2>'/2<v:2>1/2  t.. 

V  »'lj     2    1  J  1  J 


(2.48; 


and    one   may  easily  generalize  Boussinesq's    phenomenol  ogi 

cal    hypothesis    as 


p(t)    =    _e    • vr 


(2.49) 


where  r  is  any  scalar  and  P  '  is  the  turbulent  flux. 

The  limiting  cases  for  the  fluid  element  displace- 
ment covariance  and  eddy  diffusivity  are  given  below.  For 
smal 1  t 


<x'.(t)x'.(t)>  -  <v:2>1/2  <v;2>1/2t2 


2.50) 


and    E..    (t)    is    proportional    to    t.       For    large    t 


<X  !  (t)X'.(t)>    =    <V!2  >1/2  <V'2  >1/2  T.  .t    +    constant      (2.51) 
l  j  i  J  "i  J 


and    E .  .  ( t )    is    constant . 
ij 

2.1.2.2      Two-particle    dispersion 

To    thoroughly    describe    turbulent    dispersion    in    a 
homogeneous    flow    field,    all    n-particle    interactions    should 
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be  included  in  the  statistical  analysis.   The  single- 
particle  dispersion  problem  has  been  discussed  in  the 
previous  section,  and  the  two-particle  dispersion  is  dis- 
cussed in  this  section.   Due  to  increasing  mathematical 
difficulties,  the  analysis  of  the  other  n-particle  inter- 
actions have  not  been  made.   Even  though  a  complete  sta- 
tistical determination  of  turbulent  dispersion  is  desirable, 
significant  and  workable  models  based  on  the  presently 
available  works  have  been  obtained. 

As  in  the  single -particle  studies,  the  two-particle 
analysis  was  first  made  for  one-dimensional  isotropic  case. 
This  work  was  done  by  Richardson  (1926).   Batchelor  (1952) 
extended  the  two-particle  analysis  to  the  general  homo- 
geneous situation  by  investigating  the  relative  displace- 
ment and  relative  velocity  of  two  fluid  particles.   Batch- 
elor shows  that  even  though  the  flow  field  may  be  assumed 
stationary,  the  relative  velocity  covariance,  <v.(t)v.(t)>, 
is  not  stationary,  it  depends  on  the  actual  elapsed  time 
t-t0.   The  asymptotic  form  of  the  covariance  is  related  to 
Eulerian  correlations  as  follows 


<vi(t)vj(t)>  +  2<u!(t)uj(t)>, 


(2.52) 


and  the  relative  displacement  asymptotic  form  is  related 
to  the  single- particle  displacement  as 


<Yi(t)Yj(t)>  -  2<Xi(t)Xj(t)>. 


(2.53) 
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Furthermore,  under  the  normality  assumption  for  the  joint 
probability  distribution  of  the  components  of  Y(t),  the 
conditional  p.d.f.  of  the  separation  of  two  fluid  particles 
satisfies  the  diffusion  equation, 


frd-tiVV  =  ^ij^^r/^iVt, 


(2.54) 


fhere    the    relative    dispersion    tensor    is    determined    by 


V  ■  ■    t 


^(<Yi(t)Yj(t)>    -    <Y.(t0)>    <Yj(tQ)>)  (2.55) 


or 


■V-  ■  (t 
i  J 


}{    (<v1(t)vj(t1)>    +    <v.(t1)vj(t)>)dt1      .  (2.56) 


Unfortunately,  the  diffusion  equation  does  not  apply  to  the 
p.d.f.  of  the  relative  concentration. 


2.1.3   Shear  flow 

The  number  of  contributions  in  which  there  is  a 
theory  developed  on  shear  turbulence  is  very  small  when 
compared  to  those  contributions  to  isotropic  turbulence. 
Obviously  this  is  because  of  the  extreme  complexity  of 
the  problems  encountered.   Though  the  results  are  still 
very  meager,  a  few  features  of  actual  shear  flow  have  been 
studied  and  give  some  understanding  of  the  nature  of  this 
type  turbulence. 
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Batchelor  (1957)  did  some  work  in  the  area  to  free 
turbulent  shear  flows.   He  extended  the  single-particle 
analysis  to  the  dispersion  of  a  scalar  to  nonhomogeneous , 
self-preserving  turbulent  flows  such  as  exist  in  steady 
turbulent  jets,  wakes,  and  mixing  layers.   Flow  fields, 
which  are  self-preserving,  retain  their  structure  with 
elapsed  time  since  only  changes  in  the  length,  time,  and 
velocity  scales  occur.   Using  this  property  of  the  flow 
field,  Batchelor  transformed  this  nonstati onary  dispersion 
problem  into  a  stationary  one  by  an  appropriate  transfor- 
mation of  the  Lagrangian  velocity  of  the  fluid  particle. 
An  important  result  is  that  the  lateral  and  longitudinal 
dispersions,  as  functions  of  time,  are  proportional  to  the 
average  width  of  the  shear  layers  at  identical  times.   Also, 
the  time  rate  of  change  of  the  dispersion  becomes  constant 
for  large  times.   No  convenient  diffusion  type  equation  was 
found  applicable  as  for  the  homogeneous  turbulence  case. 

The  study  of  turbulent  shear  flow  has  been  mainly 
directed  at  correlations  between  velocities  and  their 
derivatives.   One  of  the  most  reported  or  discussed  studies 
of  this  type  is  that  of  Laufer  (1954).   Laufer  studied  the 
fully  developed  turbulent  shear  flow  in  a  large  pipe.   He 
not  only  investigated  the  second  and  third  order  correla- 
tions but  some  fourth  order  correlations  as  well.   The 
importance  of  a  detailed  knowledge  of  flow  conditions  near 
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the  wall  is  clearly  suggested.   The  various  turbulent 
energy  rates,  such  as  production,  diffusion,  and  dissipa- 
tion, were  shown  to  reach  a  maximum  at  the  edge  of  the 
viscous  sublayer.   It  was  also  found  that  kinetic  energy 
is  transferred  away  from  the  edge  of  the  viscous  sublayer 
while  an  equally  strong  movement  of  pressure  energy  is 
transferred  toward  the  sublayer. 


2.1.4   Interaction  of  molecular  and  turbulent  dispersion 

The  effects  of  turbulent  dispersion  and  molecular 
diffusion  are  not  easily  separated.   Taylor  (1935)  suggested 
that  the  two  are  statistically  independent,  so  that  the  vari- 
ances due  to  these  two  effects  are  additive.   Townsend  (1954 
and  Batchelor  and  Townsend  (1956)  showed  that  the  two  effects 
are  synergistic  in  a  homogeneous  turbulent  field.   They 
concluded  that  for  short  times,  the  dispersion  of  a  scalar 
is  accelerated  since  the  stretching  and  rotational  effects 
of  the  turbulence  would  enhance  the  molecular  diffusion. 
It  was  decided  that  the  total  transport  of  a  scalar  depends 
on  the  interaction  of  the  two  modes  of  transport  in  a  de- 
pendent fashion.   Townsend  confirmed  his  predictions  experi- 
mentally.  He  studied  the  heat  transport  behind  a  source  in 
homogeneous  turbulence,  and  in  addition,  indicated  adjust- 
ments were  needed  on  Uberoi  and  Corrsin's  (1953)  studies 
of  a  similar  transport  problem.   The  adjustments  had  no 
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noticeable  effect  on  the  asymptotic  eddy  diffusivity  ob- 
tained by  Uberoi  and  Corrsin.   Mickelsen  (1960)  also  found 
the  accelerated  molecular  diffusion  is  negligible  for  large 
times  in  his  experimental  measurements  of  mass  transport  in 
the  f 1 ow  field. 

Later,  Saffman  (1960,  1962)  added  to  the  analysis 
of  Townsend  and  Batchelor.   He  pointed  out  that  the  molecu- 
lar diffusion  reduced  the  turbulent  dispersion  since  the 
scalar  within  the  original  fluid  particle  is  diluted  by 
the  molecular  diffusion.   However  the  two  effects,  when 
not  negligibly  small,  tend  to  balance  themselves.   Saffman 
(1960)  and  Okubo  (1967)  have  shown  that  the  interaction  re- 
duces the  dispersion  relative  to  the  origin  from  the  value 
it  would  have  had  if  the  two  processes  had  been  independent 
and  additive.   Saffman  claimed  that  his  predictions  were  in 
qualitative  agreement  with  the  measurements  of  Mickelsen 
(1960).   Both  Saffman  and  Mickelsen  suggested  that  further 
experimentation  was  needed  to  resolve  this  problem  for 
large  times.   In  general  the  small  absolute  value  of  the 
molecular  diffusion  relative  to  the  turbulent  dispersion 
allows  one  to  neglect  the  purely  molecular  diffusion. 
Monin  and  Yaglom  (1971)  present  an  excellent  survey  plus 
supplementary  results  covering  this  area  of  study. 
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2 . 2   Lagrangian  Characteristics  of  Turbulent  Flow 

The  fundamental  analytical  approaches  of  Taylor 
(1921)  and  Batchelor  (1949)  are  basically  derived  from 
the  Lagrangian  statistics  of  a  turbulent  homogeneous  flow 
field.   Both  experimental  and  numerical  efforts  have  been 
made  to  obtain  the  Lagrangian  characteristics  pertinent  to 
these  studies.   Ideally,  the  statistics  should  be  generated 
from  numerous  individual  fluid  point  paths.   Numerically, 
the  generation  of  fluid  point  paths  is  straightforward  once 
the  statistical  parameters  and  model  have  been  decided  upon. 
However,  to  experimentally  track  or  follow  fluid  particles 
would  be  a  very  difficult  if  not  impossible  task.   Approxi- 
mate Lagrangian  statistics  have  been  obtained  from  experi- 
ments which  make  possible  the  visualization  of  the  paths  of 
small  particles  suspended  in  a  turbulent  flow  field.   These 
small  particles  consist  of  a  separate  phase  and  retain  their 
identity,  and  their  paths  are  thought  to  closely  approximate 
the  paths  of  fluid  points.   In  addition,  attempts  have  been 
made  to  extract  Lagrangian  statistics  from  the  concentration 
and  temperature  profiles  generated  by  the  turbulent  disper- 
sion of  mass  and  heat,  respectively. 


2.2.1   Lagrangian  numerical  simulations 

Several  different  stochastic  models  of  turbulent 
dispersion  have  been  used  to  obtain  the  different  Lagrangian 
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statistics  which  are  of  interest.   Because  of  the  mathe- 
matical difficulties  and  the  number  of  repetitive  calcula- 
tions involved  in  a  numerical  simulation,  the  digital 
computer  is  required  to  expedite  the  solution.   The  capacity 
and  speed  of  the  computer  used  often  limits  the  complexity 
and  size  of  the  simulation  to  be  made. 

Using  a  Monte  Carlo  procedure,  Kirmse  (1964)  simu- 
lated turbulent  dispersion  in  pressure-driven  nonhomogeneou s 
pipe  flow.   A  proportionality  parameter  was  used  to  relate 
the  variance  of  the  Lagrangian  acceleration,  a.,  to  measur- 
able Eulerian  variances  in  the  expression 


a .    =  do  o  . 
v     u  u 


(2.57) 


where  o   =  variance  of  Eulerian  velocity 

u 

o-  =  variance  of  Eulerian  acceleration 
u 


The  parameter  "B"  is  similar  to  Mickelsen's  (1955)  Eulerian 
to  Lagrangian  transformation  parameter.  By  trial  and  error 
the  parameter  "B"  was  adjusted  until  concentration  profiles 
obtained  in  the  simulation  coincided  with  available  experi- 
mental data.  From  a  simple  random  walk  model,  Patterson 
and  Corrsin  (1965)  simulated  a  one-dimensional  dispersion 
process.  A  homogeneous  velocity  field  was  assumed  and  both 
Lagrangian  and  Eulerian  temporal  velocity  correlations  were 
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obtained.   The  model  generated  results  which  appear  to 
represent  the  expected  turbulent  behavior. 

Single-particle  paths  and  their  Lagrangian  sta- 
tistics were  simulated  by  Kraichnan  (1970)  using  a  normal 
multivariate  velocity  distribution.   An  incompressible, 
stationary,  isotropic  flow  field  was  assumed.   The  velocity 
field  was  generated  from  a  set  of  Fourier  components  ob- 
tained from  Eulerian  correlations  and  spectra.   Dispersions, 
velocity  correlations,  and  eddy  d i f f usi vi ti es  were  computed 
using  the  simulation.   The  results  were  consistent  with 
Taylor's  earlier  work. 

Deardorff  and  Peskin  (1970)  generated  single-  and 
double-particle  trajectories  using  a  three-dimensional 
numerical  model  of  shear  flow  within  a  channel.   The  study 
was  limited  to  regions  outside  of  the  viscous  sublayer  for 
a  large  Reynolds'  number.   The  model  used  was  developed  by 
Deardorff  (1970)  and  was  based  on  averaging  the  Navier- 
Stokes  equations  over  grid  volumes  by  the  method  of  Reynolds 
(Reynolds,  1895;  Lamb,  1932).   After  generating  a  fully  de- 
veloped velocity  field,  forty-eight  fluid  particles  at  dif- 
ferent initial  points  within  the  turbulent  core  were  selected 
and  followed  for  two  hundred  time  steps.   A  new  set  of  fluid 
particles  was  then  selected  with  the  same  initial  starting 
points  and  the  procedure  repeated.   An  ensemble  of  ten  such 
sets  of  forty-eight  particles  was  used  to  compute  the 
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Lagrangian  temporal  mean  square  displacements  and  velocity 
autocorrelation  coefficients  for  single-particle  tracks. 
The  decay  of  the  two-particle  velocity  correlation   coef- 
ficients toward  zero  was  very   slow,  and  consequently,  the 
two-particle  mean  square  separation  did  not  reach  the 
asymptotic  value  of  twice  the  single-particle  mean  square 
displacement  predicted  by  Batchelor  (1952).   The  Lagrangian 
time  macroscales  of  the  two-particle  correlation  were  five 
to  twenty-five  times  greater  than  the  corresponding  single- 
particle  macroscales. 

A  numerical  stochastic  model  depending  on  only 
three  parameters  was  proposed  by  Sullivan  (1971)  for  fully 
developed  turbulent  flow  in  an  open  channel.   The  two-dimen- 
sional model  uses  only  the  bulk  velocity,  ub>  the  local 
Eulerian  space  macroscale,  Ar,  and  the  standard  deviation 
of  the  local  lateral  velocity  fluctuations,  a      .   With  an 

uy 

initial  lateral  velocity  equal  to  the  local  standard  devi- 
ation, a  fluid  particle  is  allowed  to  retain  its  velocity 
for  a  time  period  of 


au  (y  ) 
y  o' 


(2.58) 


The  longitudinal  and  lateral  displacements  at  the  end  of 
the  time  period  t  are 
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Ar     Y(t)  +  t) 

Ax  =  IT-  f(f\    %dy 

vy   y(V 


(2.59) 


and 


Ay  =  +— — x  =  +  Ac 

-   \      -  E 


(2.60) 


respectively,  where 


tk  =  (k  -  l)x 


2.61) 


From  the  new  position,  the  values  of  uh,  Ar,  and  a   are 

d    t         Uy 

corrected,  and  the  sign  of  the  next  lateral  displacement 
is  randomly  assigned.   By  repeating  this  procedure,  particle 
trajectories  of  any  time  length  could  be  generated.   Sulli- 
van generated  one  thousand  particle  tracks  for  ten  initial 
lateral  points  distributed  evenly  on  the  height  of  the  channel. 
Ensemble  averages  of  the  longitudinal  concentration  profiles 
for  given  elapsed  times  and  the  temporal  Lagrangian  velocity 
autocorrelation  coefficients  were  obtained.   In  view  of  the 
simplicity  of  the  simulation,  the  agreement  with  experimental 
data  is  surprisingly  good, 

Gielow  (1972)  numer i cally  s tudi ed  shear  flow  between 
parallel  plates  using  a  stochastic  model  based  on  Kirmse's 
(1964)  model.   Also  by  using  the  mean  Lagrangian  accelera- 
tion obtained  from  the  Reynolds'  stresses,  deterministic 
velocity  and  position  loci  were  found.   By  imposing  the 
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condition  that  these  deterministic  loci  are  closed  curves 
in  a  coordinate  system  connected  by  the  local  mean  velocity, 
relationships  were  obtained  to  supplement  and  to  check  the 
consistency  of  available  experimental  data  used  in  the 
stochastic  model.   In  the  model,  the  time  derivative  of 
the  velocity  signal  was  chosen  from  a  probability  distri- 
bution and  integrated  to  give  the  velocity.   Statistical 
and  spectral  analyses  of  the  simulated  velocity  records 
were  made  to  show  the  suitability  of  the  model.   Lagrangian 
time  correlations  of  the  nonhomogeneous  pressure-drive  flow 
were  computed  and  compared  to  the  same  correlations  for 
homogeneous  flow  computed  from  the  stochastic  model.   Be- 
cause of  the  confining  walls  the  eddy  diffusivity  obtained 
from  the  Lagrangian  correlations  was  clearly  time-dependent. 
The  model  was  also  used  to  simulate  the  turbulent  transport 
of  a  scalar  such  as  heat  or  mass.   The  results  for  the 
steady  state  simulations  were  consistent  with  available 
experimental  data , 


2.2.2   Lagrangian  experimental  studies 

The  Lagrangian  characteristics  of  fluid  flow  have 
been  investigated  experimentally  by  many  researchers.   The 
experiments  performed  in  this  area  used  "tagged"  fluid 
particles  to  observe  both  microscopic  and  macroscopic 
turbulent  activity.   "Tagging"  is  accomplished  by  the 
release  of  a  visible  foreign  material  into  the  flowing 
fluid.   The  bulk  of  the  experimental  investigations  have 
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been  macroscopic  in  nature,  i.e.  due  to  the  type  experi- 
ment or  objectives,  the  attention  was  focused  on  the  over- 
all effect  of  the  large  scale  nature  of  the  turbulent 
activity.   The  microscopic  type  of  experiments  deal  mainly 
with  the  individual  fluid  particles  and  their  paths  through 
the  f 1 ow  field, 

Hagen,  in  1854,  appears  to  be  the  first  to  report 
the  use  of  this  "tagging"  technique  to  study  fluid  particle 
motion  (Rouse  and  Ince,  1963).   By  observing  the  behavior 
of  suspended  particles,  Hagen  discovered  the  existence  of 
two  different  modes  of  flow.   Reynolds  (1883)  also,  in  his 
well-known  experiments  using  dye  injection,  demonstrated 
vividly  the  gross  Lagrangian  characteristic  of  a  flowing 
fluid.   Prandtl  and  Tietjens  (1934b)  discussed  and  pre- 
sented results  of  their  investigations  using  visual  tech- 
niques for  liquid  and  gas  flows. 

In  confined  flows,  the  turbulent  activity  has  been 
visually  observed  in  all  regions  of  the  flow  field,  the 
core  region  as  well  as  the  wall  region.   Turbulence  has 
even  been  observed  inside  the  viscous  sublayer  using  the 
particle  "tagging"  technique.   Fage  and  Townend  (1932) 
and  Fage  (1936),  using  a  microscope  to  view  suspended 
particles  in  a  confined  turbulent  fluid,  observed  signifi- 
cant turbulence  in  the  viscous  sublayer.   In  addition  to 
observing  the  qualitative  nature  of  the  turbulence,  Fage 
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and  Townend  made  quantitative  computations  of  local  turbu- 
lent intensities.   More  recently,  Corino  and  Brodkey  (1969) 
visually  observed  the  turbulent  fluid  motions  very  near  a 
pipe  wall,  including  the  viscous  sublayer,  using  a  high- 
speed movie  camera  moving  with  the  flow.   This  study  was 
made  in  an  effort  to  establish  a  physical  picture  which 
related  the  turbulent  activities  near  the  wall  to  the 
generation  of  turbulence  and  the  turbulent  transport  pro- 
cesses.  With  an  apparent  consistent  average  period,  fluid 
elements  were  observed  to  be  ejected  outward  into  the 
turbulent  core  from  a  thin  region  adjacent  to  the  viscous 
sublayer.   These  ejections  and  resulting  fluctuations  were 
the  most  important  feature  of  the  wall  region  and  were 
believed  to  be  a  factor  in  the  generation  and  maintenance 
of  turbulence.   Other  studies  of  this  type  have  been  done 
by  Kline  and  Runstadler  (1959),  Nedderman  (1961),  Kim  et. 
a]_.  (1971  )  ,  and  Grass  (1971  )  . 

Very  few  Lagrangian  experiments  of  the  microscopic 
nature  have  been  performed,   The  reason  for  the  lack  of 
experimental  work  in  this  area  is  most  likely  because  of 
the  laborious  and  tedious  chore  of  following  the  particles 
Even  of  the  experiments  done,  the  main  interest  was  in  the 
Lagrangian  correlation  of  the  velocity  of  a  solid  particle 
The  limiting  case  of  the  particle  velocity  correlation  is 
the  Lagrangian  correlation  function  discussed  by  Taylor 
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(1921).   Studies  of  the  movements  of  neutrally  buoyant 
particles  were  made  by  Vanoni  and  Brooks  (1955)  and  later 
extended  and  improved  by  Frenzen  (1953).   Droplets  of 
nitrobenzene  and  olive  oil  were  injected  into  the  grid- 
generated  turbulence  of  a  water  tunnel.   They  made  mea- 
surements of  successive  displacements  of  the  particles 
and  adjusted  the  velocity  data  to  account  for  decay  of 
the  turbulent  intensity.   The  Lagrangian  velocity  auto- 
correlation function  was  computed  from  the  corrected  data. 
Kennedy  (1965)  measured  the  dispersion  of  1250-micron  soap 
bubbles  and  700-micron  and  900-micron  polystyrene  beads  in 
the  grid-generated  turbulence  of  a  vertical  wind  tunnel. 
Kennedy  did  not  measure  particle  velocities  directly,  but 
inferred  them  from  the  initial  slope  of  the  dispersion 
curve.   Snyder  and  Lumley  (1971)  made  some  experimental 
measurements  of  the  particle  velocity  autocorrelation 
function.   Several  types  of  particles  were  used  in  dis- 
persion experiments  within  a  wind  tunnel  designed  after 
Kennedy's.   Within  the  experimental  error,  the  particle 
velocity  correlations  coincided  with  the  Eulerian  spatial 
correlations,  when  the  separation  was  made  d imensi onl ess 
by  dividing  by  the  integral  scales.   Data  generated  using 
hollow  glass  beads,  were  thought  to  be  a  good  estimate  of 
the  Lagrangian  fluid  properties.   The  Lagrangian  time 
integral  scale  was  found  to  be  approximated  byA[:/u',  where 
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A   is  the  Eulerian  integral  scale  and  u'  is  the  turbulent 
i  ntens  i  ty . 

2.2.3   Eulerian-Lagrangian  transformation 

A  fundamentally  correct  theory  of  turbulent  disper- 
sion of  fluid  points  is  that  of  Taylor  (1921).   The  major 
reason  this  theory  has  not  been  put  to  common  use  is  that 
it  is  formulated  in  Lagrangian  terms  rather  than  the  much 
more  easily  measured  Eulerian  variables.   As  was  indicated 
in  the  previous  section,  quantitative  measurements  of  the 
Lagrangian  characteristics  are  laborious  and  tedious,  par- 
ticularly since  reliable  ensemble  averages  required  that 
numerous  individual  particle  paths  be  investigated.   Under 
these  circumstances  it  is  advantageous  to  develop  a  the- 
oretical, or  empirical  if  need  be,  relationship  between 
the  Lagrangian  variables  and  the  more  easily  measured 
Eulerian  variables.   As  a  start  in  this  direction,  Taylor 
(1938)  proposed  that  the  changes  at  a  fixed  Eulerian  point 
are  simply  due  to  the  passage  of  an  unchanging  pattern  of 
turbulent  motion  through  that  point.   This  concept  has  been 
designated  as  frozen  turbulence  or  Taylor's  hypothesis.   An 
immediate  result  of  this  hypothesis  is  that  time  is  given 
by  the  displacement  divided  by  the  mean  velocity,  -x/  U  , 
and  in  derivative  form 


J 3_ 

<u>  at  ' 


(2.62) 


36 


Taylor  considered  the  hypothesis  applicable  only  to  uniform 
low  intensity  turbulence.   Lin  (1953)  and  Uberoi  and  Corrsin 
(1953)  tested  more  fully  the  conditions  under  which  Taylor's 
hypothesis  is  useful.   For  isotropic  turbulence  and  large 
Reynolds'  numbers,  Lin  found  the  accuracy  of  Taylor's  hypo- 
thesis can  be  estimated  by 


<(d£)  > 
Mdt; 
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For  shear  flow,  the  validity  of  the  concept  is  less  clear. 
Laufer  (1954)  and  others  indicated  experimentally  the  hy- 
pothesis can  be  extended  to  fully  developed  shear  flow 
outside  the  viscous  sublayer.   An  important  application  of 
Taylor's  hypothesis  permits  the  analysis  of  Eulerian  data 
obtained  from  a  truly  Lagrangian  type  dispersion  experiment. 
Numerous  Lagrangian  studies  have  been  made  of  turbulent 
dispersion  from  continuous  sources  in  uniform  turbulence. 
By  applying  Taylor's  hypothesis  to  the  above  type  experi- 
ments in  which  Eulerian  measurements  are  obtained  at  fixed 
points  downstream  of  the  source,  temporal  mean  functions 
can  be  approximated.   Each  fixed  position  is  assumed  to 
correspond  to  a  certain  elapsed  time  of  a  fluid  particle 
flight  from  the  originating  source. 

Taylor's  (1921)  theory  reduces  the  dispersion  pro- 
blem to  that  of  determining  the  Lagrangian  autocorrelation 
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function.   The  Lagrangian  velocity  correlation  may  be 
extracted  from  the  mean  square  displacement  function  by  a 
double  differentiation  as  suggested  by  Equation  (2.14). 
In  experiments  in  which  the  fluid  particles  were  not  followed 
individually,  the  autocorrelation  function,  if  obtained, 
was  found  through  the  questionable  step  of  doubly  differen- 
tiating the  second  moment  of  the  distribution  of  particles. 
Graphical  differentiation  is  always  difficult,  but  the 
double-differentiation  of  an  empirical  curve,  even  one  in 
which  the  scatter  of  the  data  is  seemingly  small,  can  give 
results  which  vary  widely  from  one  analyst  to  another. 
Burgers  (1951)  proposed  that  the  Lagrangian  velocity  corre- 
lation may  be  approximated  from  an  easily  measured  Eulerian 
correlation  as  follows 


(^)2>  =  <u>2  «§>  > 
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It  has  been  suggested  by  several  authors  (Mickelson, 
1955;  Baldwin  and  Walsh,  1961;  Hay  and  Pasquill,  1959; 
Kirmse,  1964)  that  Lagrangian  and  Eulerian  velocity  corre- 
lations are  similar  in  shape  and  one  can  be  used  to  generate 
the  other  by  scaling  their  arguments.   For  homogeneous  tur- 
bulence, Mickelsen  (1955)  investigated  the  similarity  between 
the  Lagrangian  velocity  autocorrelation  function  and  the 
Eulerian  longitudinal  velocity  correlation.   Assuming  a  pro- 
portionality between  the  Lagrangian  and  the  Eulerian  integral 
length  scales;  i.e. 
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A, 
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and  through  the  definitions  of  the  two  length  scales,  it  was 
derived  that 


C 


,<U'2>1/2T 


(2.66) 


where 


z,    is  the  separation  distance  between  velocity  vectors 

in  the  Eulerian  correlation,  and 
x  is  the  elapsed  time  for  the  Lagrangian  correlation. 
Kirmse  (1964)  further  deduced  that  the  above  implies 


a .  =  Boo, 
v     u  0 


(2.67) 


The  value  of  "B"  was  determined  by  Mickelsen  to  be  approxi- 
mately 0.6  for  relative  turbulent  intensities  of  about  0.03. 
This  value  of  "B"  was  determined  by  comparing  the  doubly 
integrated  Eulerian  longitudinal  correlation  evaluated  from 
hot-wire  anemometry  data  with  the  corresponding  value  of  the 
Lagrangian  turbulent  spreading  coefficient  evaluated  from 
data  collected  on  the  dispersion  of  helium  into  air.   Baldwin 
and  Mickelsen  (1963)  used  the  same  value  of  "B"  to  fit  heat 
dispersion  data.   Hay  and  Pasquill  (1959)  related  the  La- 
grangian velocity  correlation  to  the  Eulerian  time  correlation 
by  assuming  a  proportionality  of  the  time  arguments,  i.e. 

t  -  Bt  (2-68) 

where  t  is  the  Lagrangian  time  and 

t  is   the  Eulerian  time. 
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The  constant  g  was  evaluated  from  experimental  data  of  the 
integral  of  the  Eulerian  time  correlation  and  the  Lagrangian 
dispersion  in  medium  scale  atmospheric  turbulence.   At  a 
relative  turbulent  intensity  of  about  0.14,  6  varied  from 
1.1  to  8.5  in  an  uncorrelated  fashion.   Baldwin  and  Walsh 
(1961)  also  investigated  the  similarity  in  the  Eulerian  and 
Lagrangian  correlation  coefficients.   They  assumed  the  re- 
lationship 

2  1/2 

-I  (2.69) 


<u'^> 


BW 


By  an  incremental  integration  technique  applied  to  both  heat 
and  mass  dispersion  data,  a  twofold  variation  of  the  empiri- 
cal factor  Km,  with  t  and  turbulent  intensity  was  observed. 
BW 

Applying  Taylor's  hypothesis,  3  can  be  related  to  Kgw  as 
fol 1 ows 


<U> 
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1/2  INBW- 


(2.70) 


Philip  (1967)  theoretically  developed  a  relationship  for 
the  ratio  of  the  Eulerian  and  Lagrangian  integral  time 
scales  for  homogeneous,  isotropic,  stationary  turbulence. 
The  experimental  results  already  mentioned  as  well  as 
that  of  Angell  (1964a,  b)  agree  with   Philip's  results. 

A  general  treatment  and  review  of  the  theoretical 
considerations  of  the  relationship  between  the  Eulerian 
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and  the  Lagrangian  velocity  correlation  functionscan  be 
found  in  Hinze  (1959),  Brodkey  (1967),  and  Monin  and 
Yaglom  (1971  )  . 


2  .  3   Experimental  Determination  of  Eddy  Pi f f usi v i ti es 

From  the  work  of  Taylor  (1921)  and  Batchelor  (1949), 
it  was  shown  that  turbulent  dispersion  has  its  grass  roots 
implanted  in  the  Lagrangian  aspects  of  the  flow  field. 
Unfortunately,  direct  measurement  of  the  physically  appeal- 
ing Lagrangian  characteristics  is  a  laborious  and  awkward 
task.   The  overwhelming  majority  of  the  experimental  data 
available  on  turbulent  dispersion  is  of  the  more  easily 
measured  Eulerian  type.   As  was  seen  in  the  previous  section 
there  is  no  completely  satisfactory  Eulerian  to  Lagrangian 
transformation  scheme.   At  this  juncture  there  are  two  princi- 
pal methods  that  may  be  used  to  analyze  the  Eulerian  measure- 
ments of  turbulent  dispersion.   One  method  is  to  assume  the 
frozen  turbulence  hypothesis  and  use  Taylor's  analysis.   The 
second  method  is  to  retain  the  Eulerian  description  of  the 
phenomenon  with  the  aid  of  Boussinesq's  hypothesis.   In 
either  case  the  eddy  diffusivity  seems  to  be  the  common 
denominator  of  the  two  approaches.   The  concept  of  the  eddy 
diffusivity  is  a  natural  extension  of  the  more  fundamental 
concept  of  the  molecular  diffusivity  and  was  introduced  as 
a  consequence  of  the  phenomenol ogi cal  approach  to  the  dispersion 
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of  a  scalar  in  turbulent  fluids.   Boussinesq  (1877)  pre- 
sented an  essay  on  the  theory  of  flowing  water,  followed 
by  a  book  (Boussinesq,  1897),  in  which  he  discusses  the 
flux  of  momentum  in  open  channels.   Generalizing  Newton's 
law  of  viscosity,  he  developed  the  concept  of  an  effective 
momentum  diffusivity,  later  to  be  referred  to  as  the  eddy 
diffusivity.   Numerous  measurements  of  the  eddy  diffusivity 
for  heat  and  mass  transfer  have  been  made  for  systems  hav- 
ing various  geometries  and  a  wide  range  of  Reynolds  numbers 
Subbotin  et  al .  (1966)  correlated  numerous  data  on  the 
radial  component  of  the  thermal  eddy  diffusivity  and  found 
an  empirical  relationship  that  fits  very  well  the  data  for 
Prandtl  numbers  from  zero  to  10,  and  for  Reynolds  numbers 
ranging  from  10,000  to  500,000.   The  relationship  depends 
only  on  the  Reynolds  number  and  the  dimensionless  radial 
position.   Simplifications  can  be  imposed  for  the  core 
region  of  the  flow.   General  discussions  of  the  eddy  dif- 
fusivity can  be  found  in  Hinze  (1959),  Brodkey  (1967),  and 
Monin  and  Yaglom  (1971). 


2.3.1   Determination  by  Lagrangian  analysis 

Many  Lagrangian  type  experiments  have  been  performed 
in  which  only  Eulerian  type  measurements  have  been  made. 
In  these  experiments  heat  or  mass  was  dispersed  from  dif- 
ferent source  configurations  by  the  turbulent  fluid  action. 
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Subsequent   temperature  and  concentration  profiles,  re- 
spectively, were  measured  at  several  stations  located 
downstream  from  the  source.   Then  by  assuming  the  frozen 
turbulence  hypothesis  and  applying  Taylor's  (1921)  La- 
grangian  analysis  to  the  mean  square  displacements  obtained 
from  the  profiles,  certain  Lagrangian  statistics  were  com- 
puted . 

The  simplest  source  configuration  is  of  course  the 
point  source.   This  configuration  is  most  adaptable  to 
mass  transfer  experiments.   Towle  and  Sherwood  (1939) 
studied  the  turbulent  dispersion  of  carbon  dioxide  and 
hydrogen  from  an  injection  point  within  the  turbulent 
core  of  pipe  flow  using  air  as  the  carrier  gas.   By  mea- 
suring concentrations  at  several  points  downstream  from  the 
injection  point,  the  eddy  di f f usi vi ti es  were  computed  for  a 
Reynolds'  number  range  from  12,000  to  180,000.   The  eddy  di f 
fusivity  was  found  to  increase  with  Reynolds'  number  and  to 
asymptotically  approach  a  constant  value.   Kalinske  and 
Pien  (1944)  obtained  eddy  di f f us i vi ti es  using  Taylor's 
analysis.   They  investigated  the  dispersion  of  liquids 
within  water  flowing  in  an  open  channel.   It  was  revealed 
in  this  study  that  the  scale  of  the  turbulence  enters  di- 
rectly into  the  eddy  diffusion  relationship  and  it  must  be 
measured  or  estimated  if  dispersion  in  turbulent  fluids  is 
to  be  predicted  accurately.   The  dispersion  from  a  point 


43 


source  within  packed  and  fluidized  beds  has  been  investi- 
gated by  Hanratty  et  al .  (1956).   These  experimenters  in- 
jected dye  into  particle  beds  having  various  bulk  densities 
ranging  from  the  pack  bed  through  the  fluidized  bed  condi- 
tion to  the  case  where  no  particles  were  present.   Eddy 
di ff usi vi ti es  were  computed  as  functions  of  elapsed  time 
using  Taylor's  analysis.   A  nonlinear  variation  of  the  eddy 
diffusivity  with  particle  Reynolds   number  was  observed. 
With  a  constant  Reynolds   number,  the  eddy  diffusivity  was 
found  to  vary  also  with  the  particle  size.   Finally,  for  a 
certain  particle  size  and  a  constant  Reynolds   number,  the 
eddy  diffusivity  increased  with  particle  density  to  a  maxi- 
mum at  a  0.35  solid  fraction  then  decreased  with  particle 
density  until  a  fixed  bed  condition  was  reached. 

Some  heat  transfer  experiments  were  made  using  a 
line  source  generated  by  stretching  a  wire  across  the  flow 
field  and  heating  it  by  an  electrical  current.   Schubauer 
(1935)  appears  to  be  first  in  the  use  of  this  technique.   He 
measured  the  lateral  spread  of  heat  from  a  line  source  in 
small  scale  isotropic  air  flow.   Thermal  eddy  di f f usi vi ti es 
were  computed  from  measured  temperature  profiles.   A  thor- 
ough investigation  of  the  dispersion  of  heat  from  a  line 
source  was  made  by  Uberoi  and  Corrsin  (1953)  using  grid- 
produced  isotropic  turbulence  within  a  wind  tunnel.   The 
mean  temperature  distributions  were  measured  for  systematic 
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variations  in  wind  speed,  size  of  turbulence-producing  grid, 
and  downstream  location  from  heat  source.   Temperatures  and 
velocity  fluctuations  were  also  measured.   A  comparison  of 
Lagrangian  and  Eulerian  analysis  for  dispersion  in  the 
non-decaying  turbulence  was  performed.   The  ratio  of  Eulerian 
to  Lagrangian  microscales  was  determined  theoretically  and 
was  shown  to  roughly  agree  with  the  experimental  results. 
Townsend  (1954)  carried  out  a  similar  experimental  study. 
Experimental  work  by  Brodkey  (1967)  supports  the  Eulerian- 
Lagrangian  relationships  developed  by  Uberoi  and  Corrsin. 
Crum  and  Hanratty  (1965)  measured  mean  temperature  profiles, 
temperature  fluctuations,  and  spatial  correlations  of  the 
temperature  fluctuations  in  the  wake  of  a  heated  wire.   The 
flow  system  was  fully-developed  turbulence  inside  of  a  pipe. 
Traverses  of  the  pipe  were  made  at  fixed  stations  down- 
stream from  the  source  to  obtain  the  measurements.   The 
mean  square  of  the  temperature  fluctuations  was  found  to 
have  a  relative  minimum  at  the  centerline  and  to  have  a 
maximum  at  about  two  thirds'  the  radius.   As  the  wall  is 
approached,  the  mean  square  fluctuation  tends  to  zero.   For 
increasing  distances  downstream  of  the  source,  the  maximum 
of  the  mean  square  temperature  fluctuation  increased  to  an 
apparent  asymptotic  value  at  large  distances. 

Other  source  configurations  have  been  studied  using 
Lagrangian  analysis.   Taylor  ( 1  954)  i  nvesti  gated  the  dispersion 
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of  a  "plug"  of  foreign  material  instantaneously  inserted 
into  a  fluid  flowing  in  a  pipe.   The  dispersion  of  the 
"plug"  for  long  elapsed  times  was  the  principal  point  of 
interest  in  his  study.   Taylor  adapted  his  one-dimensional 
single-particle  analysis  to  pipe  flow.   In  spite  of  the 
simplicity  of  the  resulting  model,  the  predictions  for 
long  elapsed  times  of  the  longitudinal  concentrations  were 
in  excellent  agreement  with  experimental  results  (Taylor, 
1954;  Batchelor  et  al .  1955).   Hanratty  (1956)  theoretically 
studied  the  transfer  of  heat  in  isotropic  turbulent  flow 
between  a  line  source-sink  pair  by  superpositioning  of 
solutions  for  point  source-sink  pairs.   Using  Taylor's 
exponential  form  of  the  Lagrangian  velocity  autocorrelation 
function  and  an  experimental  thermal  eddy  diffusivity,  Han- 
ratty calculated  temperature  profiles  between  the  source 
and  sink.   The  results  agree  well  with  experimental  data 
for  heat  transfer  between  parallel  plates. 


2.3.2   Determination  by  Eulerian  analysis 

The  transformation  of  Eulerian  data  to  Lagrangian 
data  is  not  necessary  if  a  complete  Eulerian  analysis  is 
used.   In  the  Eulerian  approach  it  is  assumed  that  the 
mean-values  of  velocity,  temperature,  and  concentration 
can  be  calculated  with  the  aid  of  a  suitably  chosen  eddy 
diffusion  coefficient  using  Boussineq's  hypothesis  or  the 
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mixing-length  theories.   In  these  theories  the  rate  of 
transfer  of  a  property  is  taken  to  be  proportional  to  the 
spatial  rate  of  change  of  the  mean  value  of  the  property. 
By  application  of  the  mass,  energy,  or  Navier-Stokes  equa- 
tions to  experimental  data,  the  eddy  diffusivity  parameter 
can  be  evaluated  and  used  in  numerical  solutions  to  these 
equations . 

As  in  Lagrangian  studies,  the  point  source  has  also 
been  used  in  Eulerian  analysis  of  turbulent  dispersion. 
Various  atmospheric  scales  of  turbulence  have  been  measured 
and  analyzed  using  one-dimensional  Fickian  type  models 
(Akerblom,  1908;  Taylor,  1915;  Schmidt,  1917;  and  Defant, 
1921).   Richardson  (1926)  summarized  the  results  of  these 
and  additional  workers.   Roberts  (1923)  broadened  the  analy- 
sis to  multidimensional  problems  by  using  a  constant  value 
of  the  mass  diffusivity  for  each  principal  direction.   Then 
Sutton  (1932)  allowed  for  the  components  to  be  positionally 
dependent  as  well.   This  model  incorporated  the  experimental 
facts  then  known.   Sutton  also  proposed  functional  forms 
of  the  Lagrangian  velocity  correlation  and  obtained  dif- 
fusivities  for  a  nonhomogenous  atmosphere  which  depended 
only  on  the  dispersion  time.   Sutton's  model  despite  its 
limited  accuracy,  is  still  one  of  the  few  theoretical 
analyses  that  accounts  for  many  atmospheric  observations 
(Monin  and  Yaglom,  1971).   Recently,  measurements  of  medium 
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and  large  scale  atmospheric  turbulence  have  been  made  by 
Hay  and  Pasquill  (1957,  1959),  Smith  (1961),  and  Angell 
(1964a, b).   Pasquill  (1962)  discusses  the  experimental 
resul ts  up  to  1 962 . 

Based  on  experimental  data,  Pasquill  (1961)  and 
Gifford  (1961)  presented  semi -empi r i cal  nomograms  which 
could  be  used  to  predict  scalar  dispersion  in  the  atmos- 
phere.  Meteorological  services  have  made  extensive  prac- 
tical use  of  the  Pasqui  1 1 -Gi f f ord  nomograms  (Pasquill, 
1971  ). 

Further  atmospheric  studies  in  the  early  fifties 
by  Lettau  (1951)  and  Davies  (1954)  have  measured  the  com- 
plete eddy  diffusivity  tensor  including  the  off-diagonal 
components  of  the  tensor. 

Experiments  using  point  sources  have  also  been 
performed  to  investigate  turbulent  dispersion  in  nonhomo- 
geneous  confined  flows.   These  studies  have  been  carried 
out  in  pipe  flow  as  well  as  in  packed  bed  flow.   Bernard 
and  Wilhelm  (1950)  studied  the  dispersion  of  a  tracer  in- 
jected in  a  gas  flowing  through  a  cylindrical  packed  bed. 
Eddy  d if fusi vi ti es  were  found  by  fitting  the  experimental 
data  to  analytical  solutions  of  the  mass  equation.   The 
analytical  solution  was  obtained  by  assuming  the  average 
axial  and  radial  diffusivity  coefficients  to  be  equal  and 
then  solving  the  system  of  equations  by  separation  of 
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variables.   The  di ff usi vi ti es  were  found  to  increase  with 

Reynolds'  number  for  the  particle  Reynolds'  number  of  5 

to  2,400  studied.   Considerable  work  has  been  done  in  this 

area  by  Fahien  and  co-workers  using  a  procedure  developed 

by  Fahien  (1954)  for  a  point  source.   The  axial  component 

of  the  eddy  diffusivity  tensor  was  neglected  on  grounds  of 

small  gradients  in  that  direction.   This  assumption  enabled 

an  analytical  solution  of  the  mass  equation  for  the  axial 

direction  to  be  obtained  using  the  separation  of  variables 

technique.   The  radial  component  of  the  eddy  diffusivity 

tensor  is  allowed  to  vary,  and  the  mass  equation  for  the 

radial  direction  is  solved  numerically. 

This  solution  technique  has  been  applied  to  both 

fully-developed  pipe  flow  and  packed  bed  flow.   Fahien 

and  Smith  (1955)  studied  dispersion  in  a  packed  bed  for 

gas  flows  by  injecting  a  tracer.   The  radial  dependence 

of  the  radial  component  of  the  eddy  diffusivity  tensor; 

E   ,  was  determined  for  flow  conditions  having  a  particle 

Reynolds'  number  ranging  from  12  to  15,000.   This  component 

was  found  to  have  a  relative  minimum  at  the  centerline,  a 

maximum  near  the  wall,  and  tended  toward  zero  at  the  wall. 

The  value  of  E   observed  at  the  centerline  was  in  good 
rr 

agreement  with  the  asymptotic  values  obtained  by  Hanratty 
et  al  ,  (1956)  for  packed  beds.   Dorweiler  and  Fahien  (1959) 
confirmed  these  results  and  extended  the  interpretation  of 
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the  data.   Using  a  modification  of  Fahien's  procedure, 
Seagrave  and  Fahien  (1961)  investigated  dye  injection 
into  water  flowing  in  a  pipe.   They  found  that  the  turbu- 
lent eddy  diffusivity  was  greater  than  the  molecular 

5 
diffusivity  by  a  factor  of  10   for  Reynolds'  numbers  of 

about  7,500.   For  the  lower  Reynolds'  numbers,  axial  dis- 
persion appears  to  increase  in  its  importance.   Similar 
behavior  was  reported  by  Roley  and  Fahien  (1960)  for  gaseous 
flows.   Additional  data  and  eddy  diffusivity  analysis  have 
been  obtained  by  Frandolig  and  Fahien  (1964)  and  Konopik 
and  Fahien  (1964)  . 

Eulerian  studies  of  heat  transfer  between  parallel 
plates  have  been  conducted  by  a  group  of  workers  between 
the  years  1947  and  1970  under  the  1  eadershi  p  of  Sage  and 
Corcoran  (Corcoran,  1948;  Corcoran  et  a! .  ,  1947,  1952, 
1956;  Cavers  et  al . ,  1953;  Schlinger  et  al . ,  1953a, b; 
Schlinger  and  Sage,  1953,  Hsu  et  al . ,  1956;  Sage,  1959; 
Venezian  and  Sage,  1961;  Chia  and  Sage,  1970).   For  Rey- 
nolds' numbers  ranging  from  8,900  to  100,000,  the  mean 
temperature  profile  of  air  was  measured  between  a  heated 
source  wall  and  the  parallel  sink  wall.   Only  the  lateral 
heat  flux  was  considered,  thereby   allowing  only  the 
lateral  component  of  the  thermal  eddy  diffusivity  tensor 
to  be  calculated.   The  dependence  of  this  component  on 
lateral  position  was  found  to  be  similar  to   dependence 
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found  for  E    in  pipe  flow  by  Fahien  and  co-workers.   The 
value  of  the  thermal  eddy  diffusivity  increases  with  in- 
creasing Reynolds   number  for  a  given  lateral  position. 
It  was  generally  greater  than  the  eddy  viscosity.   It  was 
found  in  the  latter  phase  of  the  work  that  the  lateral 
heat  transfer  was  influenced  significantly  by  changes  in 
viscous  dissipation.   Since  the  work  of  Venezian  and  Sage 
(1961),  the  eddy  diffusivity  results  have  been  corrected. 
The  viscous  dissipation  effect  (Chia  and  Sage,  1970)  mani- 
fests itself  as  a  sigmoid  increase  of  the  lateral  heat 
flux  away  from  the  source  wall  and  maximizes  at  the  sink 
wall.   The  correction  required  increases  with  Reynolds 
number,  being  about  four  per  cent  at  a  Reynolds   number 
of  40,000  and  thirteen  per  cent  at  a  Reynolds   number  of 
100,000.   When  the  viscous  dissipation  effect  is  accounted 
for,  the  thermal  eddy  diffusivity  is  found  to  be  symmetric 
about  the  centerline.   Chia  and  Sage  (1970)  summarize  the 
results  of  this  excellent  series  of  works. 

Heat  transfer  within  a  uniformly  heated  pipe  has 
been  studied  in  the  fully  developed  flow  of  air  and  water. 
The  principal  investigations  of  air  flow  systems  were  done 
by  Deissler  and  Eian  (1952),  Nunner  (1956),  Schleicher 
(1958),  Abbrecht  and  Churchill  (1960),  Tanimoto  and  Han- 
ratty  (1963),  Ibragimov  et  al  .  (1969,  1971),  and  Bourke  and 
Pulling  (1970).   Both  the  thermal  entrance  region  and  the 


isothermal  region  have  been  investigated,   Mean  tempera- 
ture profiles,  mean  square  temperature  fluctuations,  and 
the  radial  thermal  eddy  diffusivity  have  been  measured. 
Axial  heat  transfer  was  neglected,  an  assumption  acceptable 
for  a  Peclet  number  greater  than  100  (Schneider,  1957). 
The  behavior  of  the  eddy  diffusivity  profile  for  the  axial 
and  radial  direction  was  similar  to  that  observed  in  above 
mentioned  studies  of  the  diffusivity.   In  these  studies 
Reynolds   numbers  ranged  from  7,000  to  71,000,  except  for 
the  work  of  Ibragimov  et  al . ;  a  Reynolds   number  of  32,500 
and  260,000  was  used.   The  correlation  between  velocity  and 
temperature  fluctuations  was   measured  by  Ibragimov  et  al .  , 
and  they  also  determined  the  axial  heat  flux.   The  ratio  be- 
tween the  local  radial  and  axial  turbulent  flux  components 
is  always  smaller  than  one,  with  a  maximum  of  about  0.70 
midway  between  the  centerline  and  the  wall.   Beckwith  and 
Fahien  (1963),  Truchasson  (1964),  Rust  and  Sesonske  (1966), 
and  Smith  et  al  .  (1966)  have  investigated  liquid  flow  systems. 
Reynolds   numbers  varied  from  5,000  to  243,000.   All  used 
water  as  the  liquid  except  Rust  and  Sesonske;  they  used 
mercury  and  ethylene  glycol.   The  corresponding  Prandtl 
numbers  are  0.0018,  10,  and  44  for  mercury,  water,  and 
ethylene  glycol,  respectively. 

Analytical  and  semi numer i cal  Eulerian  investigations 
of  the  mass  and  energy  equations  have  been  made  for  fully 
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developed  turbulent  flow  within  a  pipe  and  between  parallel 
plates.   Beckers  (1956),  Schleicher  and  Tribus  (1957),  Spar- 
row et  al  .  (1957)  studied  heat  transfer  for  pipe  flow  using 
empirical  correlations  for  the  thermal  eddy  diffusivity 
and  considered  only  the  radial  dependence.   Gielow  (1965) 
considered  fully  developed  flow  in  a  pipe  having  a  constant 
temperature  at  the  wall.   He  performed  a  numerical  analysis 
of  the  energy  equation  to  obtain  local  mean  temperature  pro- 
files.  Kakac  and  Paykoc  (1968)  did  a  similar  study  for  flow 
between  parallel  plates  with  constant  heat  flux  and  constant 
temperature  at  the  walls.   Again,  since  experimental  data 
are  available  for  only  the  radial  direction,  only  the  radial 
dependence  was  considered.   Semi  numerical  investigations 
of  mass  transfer  was  made  by  Russo  (1965).   Using  experi- 
mentally available  data  for  various  components  of  the  eddy 
diffusivity  tensor,  Russo  obtained  solutions  to  the  mass 
equation . 


CHAPTER  3 
EXPERIMENT  AND  DATA  ANALYSIS 

3  . 1   Experiment 

The  experiment's  objective  required  the  generation 
of  fully  developed  turbulent  flow  with  subsequent  observa- 
tion of  individual  fluid  point  paths.   To  accomplish  this 
objective,  liquid  flowing  in  a  glass  pipe  was  chosen  as  the 
flow  system.   Small,  spherical  particles  were  suspended  in 
the  liquid  and  illuminated  so  their  paths  were  made  visible. 
Ideally,  the  particles  should  be  small  in  size  and  of  the 
same  density  as  the  fluid.   However,  the  requirement  of 
visibility  and  availability  necessitated  a  compromise.   The 
particles  selected  were  hollow  glass  micro  spheres  with 
diameters  in  the  range  of  10  to  100  microns.   Particle 
density  was  about  0.34  grams  per  cubic  centimeter.   These 
solid  particles  were  assumed  to  approximate  fluid  points. 

The  apparatus  used  to  photograph  the  particle  motion 
consisted  of  three  basic  systems;  (1)  the  flow  loop  in  which 
the  turbulent  motion  to  be  observed  was  generated,  (2)  a 
lighting  system  for  illuminating  small  suspended  particles 
depicting  the  fluid  motion,  and  (3)  a  unique  camera  for 
photographing  the  particle  motions  over  extended  distances. 
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Figure  3.1  shows  the  integrated  arrangement  of  these  three 
systems . 


3.1.1   Flow  loop 

To  insure  fully  developed  turbulent  flow,  several 
precautions  were  taken  in  the  design  of  the  flow  loop. 
These  included  (1)  a  constant  head  fluid  source,  (2)  a 
baffled  entrance  section  to  reduce  disturbances  caused  by 
piping,  and  (3)  a  sufficiently  long  entrance  region  to  the 
test  section  to  allow  the  flow  to  develop  fully. 

A  constant  head  source  was  necessary  to  eliminate 
the  pressure  fluctuations  created  by  the  pump.   Process 
fluid,  tri chl oroethyl ene  (TCE),  was  pumped  from  a  receiv- 
ing and  storage  drum,  through  a  filter,  and  to  an  elevated 
drum  for  a  constant  head  supply  of  fluid.   The  pump  was  a 
centrifugal  type  equipped  with  a  mechanical  seal  which 
provided  an  excellent  seal  for  the  TCE  and  reduced  the 
potential  of  contaminating  the  process  fluid  with  foreign 
oils  and  greases  present  in  other  types  of  sealing  arrange- 
ments.  A  pot  type  filter  assembly  was  used  to  remove 
solids  and  contained  five  0.304  meter  (1  foot)  long  Fram 
Corporation,  CF10EIH  filter  elements.    These  filter  ele- 
ments removed  all  particles  above  5  microns  in  diameter. 
Both  the  receiving  and  elevated  drums  were  constructed 
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from  standard  55-gallon  drums  with  appropriate  pipe  coup- 
lings welded  into  their  walls  to  accommodate  piping.   An 
internal  standpipe  was  installed  inside  the  elevated  drum. 
By  pumping  the  TCE  to  this  drum  at  a  faster  rate  than  was 
required  for  the  test  section,  the  excess  fluid  flowed 
down  the  standpipe  maintaining  a  constant  level  in  the 
drum.   The  receiving  drum  was  vented  to  the  elevated  drum 
so  that  vapors  entrained  in  the  downcoming  fluid  would 
seek  their  level  in  the  upper  drum  without  hindering  liquid 
flow  down  the  standpipe.   Sight  glasses  in  both  drums  pro- 
vided a  convenient  method  for  examining  the  quantity  of 
fluid  in  each  drum. 

From  the  bottom  of  the  elevated  drum,  the  TCE  flowed 
down  a  0,0508  meter  (2  inch)  diameter  aluminum  pipe  to  a 
reducing  section  before  flowing  into  the  horizontal  test 
section,  a  glass  pipe  4.56  meters  (15  feet)  long  with  an 
inside  diameter  of  0.0254  meter  (1  inch).   The  reducing 
section  provided  a  reduction  of  pipe  size  as  well  as  damp- 
ing of  abnormal  fluid  motion  created  by  the  piping  configur- 
ation.  Internal  straightening  vanes  were  installed  just 
before  the  reducer.   The  flow  cross  section  was  divided  into 
axial  quadrants  by  the  vanes.   This  was  thought  adequate  to 
dissipate  any  large  scale  disturbances  generated  by  the 
piping  configuration.   Two  fused  glass  joints  were  present 
in  the  test  section,  one  at  the  0.304  meter  (1  foot)  position 
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and  the  other  at  the  1.52  meter  (5  foot)  position.   Obser- 
vations were  made  over  1.52  meters  of  pipe  starting  1.22 
meters  (4  feet)  from  the  last  fused  glass  joint.   This 
arrangement  provided  48  diameters  of  straight  unobstructed 
pipe  for  flow  development.   The  first  1.52  meters  of  pipe 
also  provided  a  flow  developing  region  even  though  two 
fused  glass  joints  were  present.   Following  the  observa- 
tion section,  the  remaining  length  of  pipe  reduced  the 
influence  of  the  exit  region.   The  full  length  of  the  glass 
pipe  was  supported  at  several  points  to  facilitate  leveling 
and  eliminate  sagging  of  the  pipe. 

To  minimize  distortion  of  the  particle  visualiza- 
tion, an   observation  section  was  constructed  so  that  the 
effects  of  the  pipe's  curvature  were  eliminated.   Approxi- 
mately 1.83  meters  (6  feet)  of  the  glass  pipe  were  enclosed 
within  two  glass  plates  and  a  supporting  frame.   Care  was 
taken  in  the  construction  of  the  enclosure  to  make  sure  the 
angle  between  the  glass  plates  was  90  degrees.   The  space 
within  the  enclosure  surrounding  the  flow  tube  was  filled 
with  glycerin  to  yield  two  undistorted  perpendicular  views. 
Glycerin  has  approximately  the  same  refractive  index  as 
glass  and  TCE,  and  thereby  circumvents  the  distortion  which 
is  present  when  viewing  a  bare  circular  glass  pipe.   It  is 
through  the  two  perpendicular  glass  plates  that  the  particles 
are  photographed,   See  Figure  3.2. 
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The  process  fluid  was  recycled  to  minimize  the 
amount  of  TCE  needed  for  proper  circulation  in  the  loop. 
The  TCE  exits  the  glass  section  through  the  side  arm  of 
a  tee  into  an  expansion  section.   The  tee  allows  passage 
of  fluid  through  the  flow  loop  while  enabling  the  flow 
tube  to  be  illuminated  axially  through  a  glass  port 
mounted  to  a  second  arm  of  the  tee.   A  0. 050 8_ meter  (2- 
inch)  aluminum  pipe  was  used  to  return  the  TCE  through  a 
rotameter  to  the  receiving  and  storage  drum.   The  rota- 
meter was  calibrated  with  TCE  and  had  a  capacity  of  1.26 
liters/sec  (20  G.p.m.),  yielding  a  maximum  Reynolds   number 
for  the  glass  pipe  of  200,000.   The  temperature  of  the  pro- 
cess fluid  was  measured  just  before  returning  the  fluid  to 
the  receiving  drum. 

A  small  stream  of  the  process  fluid  was  passed 
through  a  silica  gel  trap  to  remove  water  from  the  TCE. 
This  minimized  corrosion  of  iron  materials  used  in  the 
construction  of  the  flow  loop. 


3.1.2.   Lighting  system 

The  suspended  particles  were  illuminated  using  a 
high  intensity  carbon  arc  lamp.   Two  0.0634-meter  (2.5-inch) 
diameter  condensing  lenses  were  used  to  colli  mate    the 
light  beam  and  enhance  the  dark  field  illumination.   By 
rotating  a  four  bladed,  0.304-meter  (12-inch)  diameter  disk 
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perpendicular  to  the  axis  of  the  light  path  and  at  the 
focal  point  between  the  two  condensing  lenses,  the  illumi- 
nation of  the  suspended  particles  was  interrupted  at  a 
constant  frequency.   By  interrupting  the  light,  a  time 
scale  was  imposed  on  the  photographic  images  obtained. 

The  chopping  frequency  was  sensed  with  the  photo- 
cell circuit  shown  in  Figure  3.3.   When  one  of  the  blades 
of  the  disk  came  between  the  arc  and  the  photo  cell,  the 
illumination  of  the  cell  was  reduced  and  its  resistance 
to  electrical  current  flow  increased.   This  increase  in 
resistance  caused  a  corresponding  increase  in  the  voltage 
drop  across  the  photo  cell.   The  reverse  occu red  when  the 
cell  was  illuminated.   The  voltage  pulses  generated  in 
this  manner  were  counted  by  a  General  Radio  model  1191 
counter  and  the  corresponding  period  displayed. 


3.1.3   Camera 

The  camera  was  unique  in  that  it  took  80  pictures 
simul taneously  --  forty  pictures  for  each  view,  along  a 
1.52-meter  (5-foot)  length  of  the  flow  tube.   This  length 
corresponded  to  60  pipe  diameters.   The  camera  construc- 
tion is  shown  in  Figure  3.4  and  is  described  below.   Forty 
0.0159-meter  (5/8-inch)  holes  were  drilled  and  tapped  in 
two  0.0761-meter  (3-inch)  wide  aluminum  channels.   Lenses, 
0.012  meter  in  diameter  and  having  a  0.138-meter  focal 
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length,  were  mounted  in  copper  bushings  and  inserted  in 
the  drilled  and  tapped  holes.   In  each  lense  mount,  a 
disk  with  a  central  0.00159-meter  (1/16-inch)  diameter 
hole  was  used  as  an  aperture.   This  size  aperture  pro- 
vided sufficient  light  to  expose  the  film  and  resulted 
in  good  image  definition  for  a  field  depth  of  about  0.0254 
meter  (1  inch).   The  two  aluminum  channels  were  attached 
to  a  0.127-meter  (5-inch)  angle  iron  to  establish  a  right 
angle  between  the  two  views  of  the  flow  tube  and  give  a 
rigid  camera  construction.   Each  lense  and  aperture  assembly 
focused  a  0.038-meter  (1  1/2-inch)  section  of  the  flow  tube 
onto  a  corresponding  frame  of  3 5- millimeter  photographic 
film.   Each  sequential  section  of  the  flow  tube  was  photo- 
graphed on  a  corresponding  sequential  frame  of  the  film. 
Partitioned  compartments  for  each  lense  prevented  double 
exposure  on  adjacent  frames.   The  film  was  mounted  on  two 
0.0761-meter  (3-inch)  wide  aluminum  channels,  one  for  each 
view  and  positioned  0.70  meter  from  the  lenses,  see  Figure 
3.5.   The  distance  from  the  glass-air  interface  to  the 
lenses  was  approximately  0.7 2- meter.   The  camera  assembly 
was  positioned  over  the  flow  tube  and  aligned  such  that 
the  projected  image  of  the  tube  was  centrally  located  on 
the  film.   Each  viewed  direction  was  at  approximately  45 
degrees  with  the  vertical,  see  Figure  3.2. 
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A  dark  room,  2,43  meters  (8  feet)  square  by  2.43 
meters  (8  feet)  high  was  constructed  around  the  camera 
area  of  the  flow  loop  so  as  to  yield  a  darkened  volume 
suitable  for  exposing  the  film.   One  shutter  was  installed 
to  admit  a  light  beam  along  the  axis  of  the  flow  tube  into 
the  darkened  room  enclosing  the  camera. 


3.1.4   Procedure 

The  procedure  used  to  photograph  the  suspended 
particles  is  given  below: 

(1)  The  system  is  charged  with  sufficient  TCE  to 
avoid  cavitation  in  the  pump,  approximately 
303  liters  (80  gal  Ions)  . 

(2)  The  pump  is  turned  on,  and  the  system  is  al- 
lowed to  reach  thermal  equilibrium  at  the 
desired  flow  rate. 

(3)  A  small  10-micron  filter  assembly  is  charged 
with  particles  by  allowing  TCE  to  flow  from 
the  pot  filter  through  the  small  filter 
assembly.   In  this  step,  particles  collected 
in  the  large  pot  filter  are  flushed  to  the 
smaller  filter  and  reused. 

(4)  By  reversing  the  flow  through  the  10-micron 
filter  assembly,  particles  are  entrained  and 
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directed  to  the  glass  tube  for  observation. 
The  flow  rate  through  the  small  filter 
assembly  is  adjusted  so  that  approximately 
10  to  20  particles  are  entrained  per  10- 
second  interval, 

(5)  Photographic  film  is  then  clamped  into  place 
on  the  aluminum  channel  pieces  and  inserted 
into  the  camera. 

(6)  With  the  shutter  closed  the  carbon  arc  lamp 
and  light  chopper  are  turned  on.   When  the 
period  at  which  the  light  beam  is  interrupted 
becomes  constant,  the  shutter  is  opened  for 
approximately  7  seconds  and  closed  again. 

(7)  The  flow  rate,  temperature,  and  chopping  period 
are  noted  and  recorded. 

(8)  The  light  source  is  turned  off,  and  light  chopper 
i  s  turned  off . 

(9)  The  film  holders  are  removed  from  the  camera, 
and  the  film  removed  for  processing  as  described 
in  the  next  section. 


3.1.5   Processing  of  photographic  film 

High  speed  photographic  film  was  used  to  photograph 
the  particles.   Kodak,  35-millimeter,  2475  recording  film 
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was  selected  because  of  its  availability  and  inexpensive 
cost.   The  ASA  number  rating  for  this  film  is  1600.   After 
exposure  of  the  film  as  described  earlier,  the  film  was 
removed  from  the  film  holder  and  processed  in  a  dark  room. 
The  film  was  processed  for  30  minutes  at  65  degrees  Fahren- 
heit in  Kodak  DK-50  developer.   The  developer  solution  was 
prepared  in  the  manner  and  in  the  proportions  given  in  the 
directions  which  came  with  the  chemical.   A  30-second  stop 
bath  was  used  to  neutralize  the  developing  chemicals.   Then 
a  10-minute  fixer  step  was  used,  followed  by  a  30-minute 
water  wash  step.   The  fixer  solution  was  prepared  as  pre- 
scribed by  Kodak  directions.   To  reduce  the  possibility 
of  drying  stains,  the  film  was  flushed  with  a  solution  of 
Kodak  Photo  Flo  200.   The  film  was  then  allowed  to  dry  for 
several  hours . 

3  .  2   Data  Analysis 

Before  the  Lagrangian  correlations  could  be  computed, 
two  transformations  of  the  data  were  necessary.   First,  the 
two-dimensional   coordinates  of  two  perpendicular  views  of 
the  same  point  in  the  flow  tube  were  transformed  into  a 
three-dimensional  Cartesian  coordinate  point.   Secondly,  the 
Cartesian  coordinates  were  transformed  into  cylindrical 
coordinates.   The  resulting  particle  points  were  then  smoothed 
and  the  Lagrangian  correlations  computed. 
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The  photographed  particle  tracks  were  scanned 
using  a  PDP-11  minicomputer  in  conjunction  with  a  X-Y 
recorder.   The  track  images  were  projected  into  the  plane 
of  the  X-Y  recorder  by  using  a  photo  enlarger.   To  obtain 
the  coordinates  of  a  particle's  position,  both  parallelax 
and  the  refraction  at  the  glass-air  interface  were  consid- 
ered.  The  algorithm  used  to  collect  the  two-dimensional 
coordinates  of  the  perpendicular  views  is  given  in  Appendix 
A. 


3.2.1   Transformation  of  photographic  data 

Using  the  laws  of  geometry  and  optics  as  applied  to 
the  photographing  system,  transformations  of  the  image  can 
be  made  to  account  for  parallelax  and  refraction.   In  the 
analysis  that  follows,  the  difference  in  refractive  index 
of  the  glass,  glycerin,  and  the  TCE  is  neglected.   Figure 
3.6  is  used  as  a  guide  to  the  following  discussion  and  a 
means  of  defining  the  nomenclature.   Note  that  the  coordi- 
nate perpendicular  to  D,  and  D~  is  measured  from  the  center- 
line  of  the  pipe  in  the  direction  of  "d",  i.e.  perpendicular 
to  the  glass-air  interface.   For  the  present  time  consider 
"d"  to  be  the  distance  from  the  glass-air  interface  to  the 
centerline  of  the  pipe.   The  analysis  is  started  by  first 
relating  the  image  in  the  photographing  plane  to  the  virtual 
image  in  the  glass-air  interface.   From  similar  triangles 
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G  =  H 

v    u 


H  =  G^  . 
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(3.1) 


Also  using  the  law  of  refraction, 


s  i  n 


(H2  ±  u2)"2 


st  n 


T^hT 


((D-H)2  +  d2)1/2 


(3.2) 


where  "n"  is  the  refractive  index  of  the  glass.   From  Equa 
tion  (3.2),  the  ratio  between  D  and  H  is  found  as 


a  =  ( i  + i 

H        ((n2  -  1)H2  *    n2.2)1'2 


From  Figure  3.6,  it  is  obvious  that 


H  i       H  n  H 
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Using  Equation  (3.1),  Equation  (3.3)  can  be  written  for 
coordinate  direction  m 


M   =  Hi 


1  + 


^((n2  -  l)Gfc  +  n 


2^2  2,1/2 


or 
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JH  =  1L 


n   -  1 


r2  .   2  2x1/2 

G   +  n  v  ) 


(3.4) 


The  length  of  "G  ",  i.e.  the  image  height  produced  on  the 
film,  is  projected  by  the  photo  enlarger  with  a  magnifica- 
tion factor  of  "M";  therefore,  the  distance  "Gm"  on  the 
film  is  magnified  to  a  length  "C^"  according  to 


C   =  MG  . 
m     m 


(3.5) 


Using  Equation  (3.4)  and  (3.5),  now  the  coordinate  posi- 
tions in  the  flow  tube,  (  D  -,  ,  D  ?  )  ,  can  be  related  to  the 
coordinates  measured  by  the  computer  in  the  plane  of  the 
X-Y  recorder,  (C,,C?),  as 


_LL  + 
Mv 


ii    2    n,2  J   2M2  2J/2 
(  (  n   -  1  )  C   +  n  M  v  ) 


(3.6) 


The  analysis  given  above  was  concerned  with  only 
one  view  of  the  particle  position;  however,  to  fix  the 
three  Cartesian  coordinates  of  the  particle  position,  two 
such  perpendicular  views  were  necessary.   Let  D-,  be  the 
coordinate  direction  corresponding  the  flow  direction,  i.e. 
the  common  coordinate  of  two  perpendicular  views.   Introducing 
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an  index  "j",  j  =  1  or  2,  to  denote  the  two  views,  Equa 
tion  (3.6)  can  be  written  for  view  "j"  as 


D  .      y  . 

C  .     M  •  v  . 
jm    j  j 


d .  -  D?0 
J     J2 


2    ,  ,rZ    J       2M2  2JJ2 
n   -l)C.+nM.v.) 


3.7 


where  the  "j"  index  denotes  the  complement  of  view  "j". 
Note 


dji 


ji 


or 
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In  order  to  align  the  images  projected  and  to  evaluate  the 

yi 
terms  ..  J   and  d.  of  Equation  (3.7),  several  datum  points 

Vj         J 

were  recorded  for  each  view  of  the  flow  tube.   An  alignment 
and  scaling  tape  with  inch  markers  on  it  was  also  placed 
within  the  enclosure.   This  tape  could  be  photographed  from 
both  directions  and  thereby  yield    an  axial  scale  as  well 
as  known  positions  within  two  corresponding  views  of  a  tube 
section.   The  two  corresponding  views  of  a  flow  tube  sec- 
tion were  projected  simultaneously  into  the  plane  of  the 
X-Y  recorder.   Figure  3.7  shows  the  arrangement  of  the 
datum  points  taken.   Adding  another  index,  "k",  to  Equation 
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Figure  3.7    Arrangement  of  Datum  Points 


(3.7)  to  denote  the  datum  point  within  a  view,  Equation 
(3.7)  becomes 


D  ■    -  u, 
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yi 
The  value  of  "n"  is  known;  however,  the  values  of  ..     , 

MJVJ 
d.,  and  D o „ .  were  evaluated  from  the  datum  points  recorded 

J         J  <-  K 

from  each  view.   The  method  used  to  compute  these  constants 
is  discussed  in  Appendix  B.   Using  Equation  (3.8),  all 
points  on  a  track  image  projected  in  the  plane  of  the  X-Y 
recorder  can  be  transformed  to  their  corresponding  Cartesian 
coordinates  within  the  flow  tube.   The  second  transformation 
from  Cartesian  to  cylindrical  coordinates  is  straightforward 
and  can  be  obtained  from  any  elementary  algebra  text. 


3.2.2   Computation  of  particle  velocity 

The  velocity  of  the  suspended  particle  was  obtained 
from  smoothed  particle  position  data.  A  fourth-order  poly- 
nomial was  used  to  represent  the  particle's  trajectory  over 
a  short  distance.  This  order  polynomial  allows  for  second 
order  variations  in  the  particle's  acceleration.  The  poly- 
nomial form  is 


*i  =  Cl  +  c2*1  +  =  3*1  +  c4l1  +  Vi 


(3.9) 
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where  i  designates  the  time  coordinate  number.   Different!' 
ating  Equation  (3.9)  gives  the  velocity  as 


V.  -  c0  +  2c, t.  +  3c„t2  +  4c, t3  . 
l     2      3  i      4  l      5  l 


(3.10 


The  coefficients  of  Equation  (3.9)  can  be  evaluated  using 
a  least-squares  fit  over  a  short  length  of  the  particle 
track  which  spans  the  point  being  smoothed.   The  coeffi- 
cients found  in  this  manner  can  now  be  used  to  compute 
the  velocity  via  Equation  (3.10). 


3.2.3   Computation  of  Lagrangian  correlations 

The  statistical  properties  of  the  random  motion 
of  a  single  fluid  particle  can  be  used  to  effectively 
evaluate  the  dispersion  of  a  scalar  property  within  a 
turbulent  flow  field.   Particle  position  and  velocity 
with  respect  to  time  is  a  starting  point  from  which  cer- 
tain statistical  correlations  can  be  generated  that  help 
describe  turbulent  dispersion.   The  flow  field  of  par- 
ticular interest  in  this  study  is  that  of  turbulent  shear 
pipe  flow.   Consider  a  fluid  particle  starting  at  time  t 

from  a  point  x   within  the  flow  field.   At  any  time  during 
r     —o 

the  particle's  flight,  its  velocity,  V_(  1 1  t  ,_x  )  ,  is  equal 
to  the  Eulerian  velocity,  U(t,x),  at  the  point  x_  which  is 
the  point  at  which  the  fluid  particle  is  located  at  time  t 
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Therefore 


v(t|to,x.0)  =  u(t,*J 


(3.11) 


or,  expressing  the  velocity  as  a  mean  plus  a  deviation 


V(t|  t  ,x.0)  =  <U(t,x:)>  +  U_'  (t,x_), 


(3.12) 


Note  that  the  deviation  in  the  fluid  particle's  velocity 
from  the  local  mean  velocity  of  the  field  is 


V'(t|to,xo)  =  V(t|tQ,x0)  -  <U_(t,xJ>=  u'(t,xj.  (3.13) 


Let 


X'(t|t  ,x  J  =  /  V'(tJ  t  ,X  )dt 


o  -o 


'o 


11  "o'^o^0! 


(3.14) 


or 


X'(t|to,i0)  =  /  V(t1|to.x_o)dt1 


/  <U(t1  ,x_)>dt] 


(3.15) 


be  the  deviation  of  the  fluid  particle's  position  after  an 
elapsed  time  t-t   from  the  position  it  would  have  had  had 
it  been  moved  only  by  the  local  mean  velocity. 

The  transport  of  a  scalar  quantity  within  a  steady 
nonhomogeneous  flow  field  has  not  been  theoretically  derived 
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at  the  present  time.   However,  from  the  studies  of  iso- 
tropic homogeneous  turbulent  flow  by  Taylor  (1921)  and 
Batchelor  (1949)  respectively,  a  functional  relationship 
between  the  eddy  diffusivity  tensor  and  the  Lagrangian 
correlations  apparently  exists.  From  the  results  of  these 
studies,  it  was  shown  that  the  eddy  diffusivity  tensor  is 


E..(t)  =  !^<x:(t)x'.(t)> 


(3.16 


E.  .(t 


£(<x:(t)vj(t)> 


:V!(t)Xj(t): 


3.17 


From  the  statistical  analysis  of  single-particle  trajectories, 
the  above  Lagrangian  correlations  can  be  obtained,  i.e. 
<X  !  (t)X  .'(t)>  and  <X  !  ( t)V  '•  ( t)>  .   Other  correlations  such  as 
<V^.  (tQ)V  '•  (t)>  and  <V  j  ( t)  V  1  ( t )  >  can  be  computed  easily  from 
the  data  on  the  particle's  position  and  velocity  at  discrete 
time  intervals. 

In  the  case  of  cylindrical  coordinates,  the  computa- 
tion of  displacements  and  velocities  needs  further  explanation 
The  application  of  Equations  (3.13)  and  (3.15)  is  straight- 
forward for  the  axial  and  radial  directions,  i.e. 


Vz(t|W  =  Uz(t|to'^o)  "  <u2(t>^)>        (3-18) 


V  '  ( 1 1 t  ,x  )  =  v  (t|t  ,x  ) 
r      '  o  — o     r v  '  o  —o 


(3.19) 
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z  '  ( t  1 1  ,x 

1  o  — 0 


=  /  U _(t,  |t  ,x  )dt,  -  /  <U  (tn,x;)>dt 


'  1  • '  o  '  — 0 


to 


'1 "-' 


"1 


3.20 


and 


r '  (t  It  x  )  =  /  V  (t,  ! t  ,x  )dt, 
1  o-o    +   r v  1  '    o  -o    1 


3.21 


The  displacement  and  velocity  deviation  in  the  angular 
direction  needs  additional  interpretation.   In  reference 
to  the  starting  point  or  origin,  x  ,  of  a  given  particle, 
the  angular  velocity  is  defined  as 


V  '  (t  |t  ,x 

0    o-o 


de(t|t0,x_0) 


dt 


(3.22 


and  the  angular  displacement  follows  from  Equation  (3.15)  as 


tit  ,x 
1  o-o 


/  V^(t1|t0,x0)dt1 


(3.23 


where  r   is  the  value  of  the  radius  at  time  t  =  t  .   Note 
o  o 

that  if  the  angular  velocity  is  defined  using  r(tft  ,x  ), 
then  when  Equation  (3.15)  is  used  to  compute  the  angular 
displacement,  the  displacement  would  not  necessarily  have 
the  proper  value  of  zero  for  a  closed  path.   Using  the 
above  definitions  for  displacements  and  deviations  in 
velocity,  the  elements  of  the  eddy  diffusivity  tensor,  E, 
are  consistent  with  Equation  (2.49)  where 
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vr 


3T  .     1  3£.     3T_, 

3r  -r   r  se^e   3z-z 


(3.24) 


for  cylindrical  coordinates,  and 


i_   is  the  unit  vector  in  the  radial  direction 
i_„   is  the  unit  vector  in  the  angular  direction 


i_z      is  the  unit  vector  in  the  axial  direction. 


CHAPTER  4 
RESULTS  AND  DISCUSSION 

Paths  of  suspended  particles  were  photographed  over 
large  distances,  and  the  turbulent  intensities  and  Lagrangian 
correlations  were  computed  from  the  paths  observed  in  a  0.0254 
meter  (1  inch)  I.D.  diameter  pipe.   The  Reynolds   number  at 
which  the  observations  were  made  was  110,000.   Particle  paths 
were  observed  for  elapsed  times  as  great  as  1.2  seconds,  cor- 
responding to  an  axial  distance  of  1.32  meters  (4.33  feet). 
Seven  particle  tracks  were  examined  with  a  total  of  1600 

time  intervals  recorded.   The  time  period  of  the  observa- 

_3 
tions  was  6.16  X  10    seconds.   The  direction  of  the  par- 
ticle paths  varied  very    little  over  this  interval  of  time, 
so  that  the  images  of  the  paths  were  nearly  straight  lines 
with  very  little  curvature. 

4.1   Results 

The  particle  track  data  were  used  to  compute  both 
turbulent  intensities  and  Lagrangian  correlations.   These 
variables  were  computed  for  five  radial  intervals.   The 
annular  zones  formed  by  the  radial  intervals  are  shown  in 
Figure  4.28.   In  view  of  the  frequent  use  of  homogeneous 
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models  for  turbulent  pipe  flow  and  the  relatively  small 
sample  size,  the  overall  averages  of  the  correlations  for 
the  five  zones  were  computed. 


4.1.1  Turbulent  intensities  and  shear  stresses 

The  computed  relative  turbulent  intensities  are 
shown  in  Table  4.1,  and  the  relative  shear  stresses  are 
given  in  Table  4.2.   The  relative  turbulent  intensities 
agree  with  those  given  by  Patterson  and  Zakin  (1967)  for 
organic  solvents  using  hot-film  anemometry.   However,  the 
radial  variations  are  not  in  agreement.   This  lack  of  agree- 
ment probably  stems  from  the  relatively  small  number  of  points 
available  to  form  ensemble  averages  for  the  five  radial  zones. 
The  relative  radial  turbulent  intensities  are  also  in  fair 
agreement  with  published  values  by  Laufer  (1954)  for  air  and 

by  Cermak  and  Baldwin  (1964)  for  water.   The  relative  tur- 

2 
bulent  shear  stress,  <V^,  ( t  )  \T  ( t  )  >/U_,  MAX  ,  is  in  good  agree- 
ment with  published  data.   The  relative  shear  stress, 

9 

<V'(t  )V'(t  )>/U   MAX,  appears  to  be  zero  as  expected.   How- 
ever,  the  nonzero  results  for  <vz(t0)VQ(t0)>/uz  MAX  cannot 
be  explained.   Both  shear  stresses  <V  '  ( t  ) V  '  ( t  ) >  and 

/    u    u    U 

<V'(t  )V'(t  )>  should  be  zero  since  there  is  symmetry  in 
r v  o   e   o 

the  tangential  direction. 

4.1.2  Lagrangian  correlations 

The  computed  correlations  for  fully  developed  turbu- 
lent pipe  flow  are  shown  in  Figures  4.1  through  4.27.  Figures 
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Table  4.1   Relative  Turbulent  Intensities 


Radial 

Interval  , 

2r/D 
o 

0.0  to  0.2 
0.2  to  0.4 
0.4  to  0.6 
0.6  to  0.8 
0.8  to  1.0 
Overal  1 


Axial 


1/2 


uz,MAX 


0.0927 

0.0799 
0.0695 
0.0869 
0.0879 
0.0819 


Radial 


1/2 


uz  ,MAX 

0,0295 
0.0315 
0,0351 
0.0397 
0.0344 
0.0359 


Tangential 

uz,MAX 

0.0367 
0.0378 
0.0400 
0.0470 
0.0451 
0.0429 


Table  4.2   Relative  Turbulent  Shear  Stresses 


Radial 

Interva  1  , 

2r/D 
o 

0.0  to  0.2 
0.2  to  0.4 
0.4  to  0.6 
0.6  to  0.8 
0.8  to  1.0 
Overal  1 


V'(t  )V!(t  ) 
z v  o   r   o 

uz  ,MAX 


0.000592 
0.000562 
0.000537 
0.001015 
0.000709 
0.000740 


^V^V 

V^oWV 

u2 

uz  ,MAX 

U2 

uz  ,MAX 

0.000543 

-0  .000160 

-0.000204 

0.000240 

0.000387 

-0.000013 

0.000851 

0.000044 

0.000581 

-0.000447 

0.000521 

-0.000035 
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4.1  to  4.3  give  a  representative  sample  of  the  results 
obtained  for  the  5  radial  zones.   However,  since  the  sample 
size  is  relatively  small,  only  the  overall  correlations 
given  in  Figures  4.4  through  4.27  will  be  discussed.   Speci- 
fic observations  about  the  correlations  are  made  below. 

The  time  period  between  observations  is  an  adjusted 
3.08  X  10   seconds  or  0.145  dimens i onl ess  time  units.   One 
dimensi onl ess  time  unit,  t*  =  1 ,  corresponds  to  the  time  re- 
quired to  traverse  the  pipe  diameter  at  a  velocity  equal  to 
the  maximum  center! ine  velocity. 

Displacement-displacement  correlations,  <X!  ( t)  X  '■  ( t )  >  . 
The  overall  displacement-displacement  correlations  are  given 
in  Figures  4.4  to  4.9.   These  correlations  are  normalized 
to  the  pipe  diameter,  D  .   The  correlations  form  a  symmetric 
second-order  tensor;  therefore,  only  the  six  elements  pre- 
sented need  to  be  considered.   The  diagonal  component  for 
the  axial  direction,  <z '  (t)>,  is  at  least  an  order  of  mag- 
nitude greater  than  the  other  components,  except  for  the 
diagonal  component  for  the  tangential  direction,  <  (  r  9  '  ( t )  )  > 

The  shape  of  the  correlations  of  the  diagonal  com- 
ponents of  the  above  tensor  are  the  form  expected.   The  cor- 
relation for  the  spread  in  the  axial  direction,  Figure  4.4, 
does  not  contain  the  time  linear  portion  predicted  by  Equa- 
tion (2.51);  however,  the  correlation  apparently  does  tend 
towards  an  inclined  straight  line  for  large  elapsed  times. 
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The  shape  of  the  curve  is  similar  to  that  reported  by 
Gielow  (1972)  for  shear  flow  between  parallel  plates.   For 
the  radial  diagonal  correlation,  <rl  (t)>  in  Figure  4.5, 
the  effect  of  the  confining  pipe  wall  is  evident  for  the 
larger  elapsed  times.   There  is  an  asymptotic  limit  for 
the  radial  displacement  correlation  because  of  the  confin- 
ing walls  of  the  system.   This  limit  is  found  by  assuming 
a  random  dispersion  of  the  particles  from  an  initial  radial 
position,  r  ,  i.e.  the  mean  square  displacement  from  r   at 
infinite  time  is 


lb. 

2tt  /   (r  -  r  )  rdr 


2 
2tt  /    rdr 

0 


2D  r     9 
3      o 


4.1) 


For  the  overall  correlation,  the  particle  may  start  at  any 
radius;  therefore,  the  asymptotic  limit  is 


\  2 

2  D  2D  r     , 

or  I        0  0   0   ,     C-  \            A 

2tt  /  [~  -  — ~ —  +  r^)r„dr 


o '    o   o    D 
=  36 


(4.2) 


2tt  /    r  dr 
0     °   ° 


Since  the  ordinate  of  Figure  4.5  is  normalized  relative  to 

D2,  the  asymptotic  limit  for  the  ordinate  of  the  figure  is 
o  K 

1/36  or  0.0278. 
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The  off-diagonal  components  of  the  tensor  are  given 
in  Figures  4.7  to  4.9.   These  correlations  are  of  the  same 
order  of  magnitude  as  the  radial  diagonal  correlation.   The 
positive  value  of  the  correlation  between  the  axial  and 
radial  components  of  the  fluid  particle  displacement  vector 
<z'(t)r'(t)>  in  Figure  4.7,  is  to  be  expected.   Due  to  the 
mean  velocity  profile,  if  a  fluid  particle  has  a  positive 
displacement  in  the  radial  direction,  it  has  moved  to  a 
flow  region  of  lower  mean  axial  velocity.   Consequently, 
the  average  effect  would  be  a  positive  axial  deviation  in 
the  particle's  displacement.   Similarly,  a  negative  radial 
displacement  leads  to  a  negative  axial  displacement,  on  the 
average.   The  product  of  the  two  displacement  combinations 
is  positive,  thereby  resulting  in  a  positive  correlation. 
The  two  remaining  correlations  are  expected  to  be  zero  be- 
cause of  angular  symmetry  in  the  pipe.   This  appears  to  be 
the  case  for  the  correlation  between  the  radial  and  tangen- 
tial components  of  the  fluid  particle  displacement  vector, 
<r'(t)r  6'(t)>  in  Figure  4.9.   However,  the  correlation  be- 
tween the  axial  and  tangential  components,  <z ' ( t)r  6 ' ( t)  > 
in  Figure  4.8,  has  a  significant  positive  nature.   This 
probably  occurs  because  of  the  relatively  small  number  of 
samples  present  in  the  correlation.   If  numerous  samples 
had  been  taken,  the  effect  of  multiple  crossings  of  the 
starting  angular  position  could  reduce  the  magnitude  of 
the  correl ati  on . 
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Displacement-velocity  correlations,  <X  !  ( t)V  '.  ( t)  >  . 
The  correlations  between  displacement  and  velocity  devia- 
tions are  presented  in  Figures  4.10  through  4.18.   These 
correlations  are  the  components  of  the  eddy  diffusivity 
tensor  for  various  mixing  length  theories.   In  addition, 
the  symmetric  off-diagonal  components  of  the  eddy  diffu- 
sivity tensor  suggested  by  Batchelor  (1949),  see  Equation 
(2.45),  are  presented  in  Figures  4.19  to  4.21.   Both  of 
these  types  of  correlations  are  normalized  relative  to  the 
product  of  the  pipe  diameter,  D  ,  and  the  maximum  mean 
axial  velocity  at  the  centerline,  U   .....  .   Again  the  diagonal 

component  of  the  tensor  for  the  axial  direction,  <z '  (t)V '  (t)> , 
is  an  order  of  magnitude  greater  than  the  other  components. 

The  diagonal  components  of  the  eddy  diffusivity 
tensor  behave  yery   much  as  expected.   In  Figure  4.10,  the 
axial  component,  <z  '  ( t )  V ' ( t ) >  ,  shows  a  gradual  increase  and 
appears  to  be  approaching  a  constant  value.   The  change  in 
direction  near  a  d imensi onl ess  time  of  9  is  probably  due  to 
decreasing  statistical  accuracy  since  the  sample  size  di- 
minished with  increases  in  the  time  coordinate.   After  a 
rapid  initial  increase,  the  radial  component  reaches  a 
plateau  and  begins  to  decline.   The  approach  of  the  axial 
component  to  its  asympototic  value  appears  to  be  more 
gradual.   The  assumption  of  a  constant  eddy  diffusivity 
coefficient  for  the  radial  direction,  often  used  in  modeling 


dispersion  in  pipe  flow,  seems  appropriate  for  a  portion  of 
the  dispersion  time.   As  seen  in  Figure  4.14,  the  radial 
diagonal  component  could  be  approximated  by  a  constant.   Of 
course,  for  large  times,  the  correlation  is  expected  to  go 
to  zero  since  the  radial  spread  is  bounded  by  the  pipe  walls. 
The  behavior  of  the  tangential  component,  <r  8'(t)Vg(t)>  in 
Figure  4.18,  seems  to  be  more  akin  to  the  behavior  of  the 
axial  component  than  the  radial  component.   The  rapid  de- 
crease in  the  tangential  correlation  after  a  dimensi onl ess 
time  of  5  is  probably  due  to  poor  statistics. 

Figures  4.11,  4.13,  and  4.19  give  the  correlations 
between  the  axial  and  radial  deviations.   The  correlations 
<z'(t)V'(t)>  and  <r  '  ( t ) V^ ( t)  >  ,  Figures  4.11  and  4.13  re- 
spectively, have  similar  shapes  which  suggest  that  the  eddy 
diffusivity  tensor  is  symmetric.   Using  the  eddy  diffusivity 
components  as  given  by  Batchelor  (1949)  symmetry  is  imposed 
by  definition,  see  Equation  (2.45).   Figure  4.19  gives  the 
axial-radial  component  of  Batchelor's  definition  of  the  eddy 
diffusivity  tensor. 

Because  of  symmetry  in  the  tangential  direction,  the 
remaining  off-diagonal  components  are  normally  assumed  to  be 
zero.   Since  the  data  were  readily  available,  these  correla- 
tions were  computed  also.   The  correlations  between  the 
radial  and  tangential  directions  are  apparently  zero  as 
indicated  in  Figures  4.15,  4.17,  and  4.21.   Figure  4.12 
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suggests  that  <z' (t)Vg(t) >  is  also  zero.   However,  the  cor- 
relation <r  9'(t)V'(t)>,  Figure  4.16,  does  not  show  a  ten- 
dency toward  zero.   It  appears  to  be  increasing.   The  axial 
tangential  component  as  defined  by  Batchelor   also  shows  a 
continuing  increase  with  time.   This  is  only  a  reflection 
of  <rQG  '  (t)V^(t)>  since 


-^  <z'(t)ro9'(t)>  =  <  z'(t)V^(t)>  +  <ro9'(t)V;_(t)> 


Considering  the  tangential  symmetry,  for  any  V'(t)  it  would 
be  equally  likely  that  the  tangential  displacement  could  be 
positive  or  negative.   This  would  result  in  a  zero  correla- 
tion.  The  problem  of  a  positive  correlation  appears  to  be 
in  the  data.   Perhaps,  since  the  sample  size  was  relatively 
small,  there  was  not  sufficient  crossing  of  the  initial 
angular  position.   The  portion  of  the  particle  path  remain- 
ing after  the  last  crossing  could  cause  an  excessive  positive 
contribution  to  <r  6'(t)VJ(t)>.   With  a  sufficiently  large  . 
sample  size,  these  correlations  should  approach  zero  for  all 
times . 

Velocity -velocity  correlations,  <VJ  (t„)V  '■  ( t)>.   The 

J. <i 1  — o — j 

correlations  between  the  velocity  deviations  at  the  initial 

time  t   and  elapsed  times  t-t   are  given  in  Figures  4.22  to 
o  o 

4.27.   These  correlations  are  the  Lagrangian  double  velocity 
correlation  coefficients  defined  by  Equation  (2.32).   Because 
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of  the  manner  by  which  the  correlations  were  generated  from 
the  particle  path  data,  the  Lagrangian  double  velocity  cor- 
relation tensor  is  symmetric.  In  view  of  this,  only  six  of 
the  nine  components  of  the  tensor  are  presented. 

The  diagonal  components  of  the  tensor,  Figures  4.22 
to  4,24,  start  at  one,  as  required,  and  decrease  toward  zero 
for  large  times.   All  three  correlations  are  identical  for 
elapsed  dimensionless  times  less  than  0.5;  after  which,  the 
radial  correlation,  <V  (t  )V  (t)>  in  Figure  4.23,  decreases 
most  rapidly.   Since  the  diagonal  components  of  the  eddy 
diffusivity  tensor  are  proportional  to  the  integrals  of  the 
respective  diagonal  components  of  the  Lagrangian  double 
velocity  correlations,  the  most  rapid  increase  in  the  radial 
eddy  diffusivity  coefficient  is  consistent  with  the  above 
observation.   Also  from  the  integral  relationship  between 
the  two  radial  components,  the  radial  Lagrangian  velocity 
correlation  would  be  expected  to  decrease  to  negative  values 
and  then  return  to  zero.   This  behavior  is  expected  since 
the  radial  eddy  diffusivity  coefficient  should  decrease  from 
its  relatively  constant  value  to  a  zero  value  at  large  times. 
Similarly,  since  the  axial  and  tangential  diagonal  components 
of  the  diffusivity  tensor,  <z'(t)V^(t)>  and  <rQe'(  t  )V^(t)>, 
tend  to  a  constant  value  for  large  times,  the  respective 
Lagrangian  double  velocity  correl a ti ens  shou 1 d  monotoni ca 1 1 y 
decrease  to  zero  for  large  times.   This  appears  to  be  the  case  for  the 
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axial  Lagrangian  correlation,  but  not  so  for  the  tangential 

correlation  <Vg(t  )V'(t)>  in  Figure  4.24.   Again  there  seems 

to  be  a  statistical  problem  related  to  the  tangential  element 

of  the  data . 

Off-diagonal  components  of  the  Lagrangian  double 

velocity  correlation  tensor  are  given  in  Figures  4.25  to  4.27 

The  correlation  between  the  radial  and  tangential  velocity 

deviations,  <V  (t  )V'(t)>,  is  zero  as  expected,  Figure  4.27. 

The  correlation  between  the  axial  and  tangential  velocity 

deviations,  <V'(t)V'(t)>  in  Figure  4.26,  is  not  zero  as 
z   o   8 

would  be  expected  from  symmetry.   Figure  4.25  gives  the  cor- 
relation between  the  axial  and  radial  velocity  deviations, 
<V'(t  )V  (t)>.   The  correlation  begins  at  a  positive  value 
and  decreases  to  zero  for  large  times.   This  behavior  can 
be  explained  as  follows.   If  the  radial  velocity  is  positive 
the  fluid  particle  would  enter  a  flow  region  of  lower  mean 
axial  velocity.   Consequently,  the  axial  velocity  deviations 
would  on  the  average  be  positive,  resulting  in  a  positive 
product  of  the  velocity  deviations.   A  similar  argument 
applies  for  a  negative  radial  velocity  deviation. 

The  interrelationships  between  the  three  different 
types  of  correlations  do  appear  to  be  satisfied  from  a  de- 
tailed study  of  Figures  4.4  through  4.27.   This  can  be  viewed 
as  a  confirmation  of  the  internal  consistency  of  the  computer 
program  used  to  generate  the  various  correlations. 
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4  .  2   Equipment  and  Procedure 

The  equipment  and  procedure  were  effective  in  ob- 
taining the  desired  data.   Turbulent  conditions  were  gen- 
erated in  the  flow  loop  without  abnormal  or  distorted  flow 
patterns.   The  suspended  particles  were  photographed  and 
their  positions  were  accurately  determined  from  the  photo- 
graphs.  Even  though  the  quality  of  the  photographs  was 
not  as  good  as  desired,  a  particle  could  be  followed  con- 
tinuously over  tens  of  diameters  of  the  flow  tube. 

Operation  of  the  flow  loop  was  satisfactory,  and 
it  is  felt  that  fully  developed  turbulence  was  established 
within  the  observation  section.   No  abnormal  bulk  flow 
characteristics  were  indicated  from  a  close  examination 
of  the  particle  path  data.   In  addition,  visual  observa- 
tions of  the  overall  flow  of  the  illuminated  particles 
revealed  no  abnormal  flow  patterns  within  the  observation 
section.   The  flow  rate  was  determined  using  a  rotameter 
calibrated  with  the  process  fluid,  tri  chloroethyl  ene  .   The 
temperature  was  also  measured  to  within  +  1°F.   The  mea- 
sured flow  rate  was  0.52+0.01  liter/sec  (8.3+0.2  GPM)  which 
corresponds  to  an  average  velocity  of  1.04+0.03  meters/sec 
(3.41+0.06  ft/sec).   Using  an  equation  form  developed  by  Pai 
(1953)  for  the  mean  velocity  profile,  a  least-squares  fit 
to  the  particle  velocity  data  yielded  an  average  velocity 
of  0.969  meter/sec  (3.18  ft/sec).   The  exact  equation  used 
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for  the  mean  velocity  profile  will  be  discussed  later.   The 
Reynolds  number  computed  from  the  above  data  is  110,000, 
a  flow  condition  well  within  the  turbulent  region  of  fluid 
f  1  ow . 

Steady  flow  conditions  were  maintained  during  the 
photographing  of  the  suspended  particle  motion.   Because 
the  flow  system  was  closed,  the  heat  generated  from  pump- 
ing caused  an  increase  in  fluid  temperature  relative  to 
the  room  temperature.   Sufficient  time  was  allowed  in  order 
to  reach  a  steady  fluid  temperature.   After  the  flow  rate 
was  adjusted,  the  rotameter  float  position  varied  only 
slightly  corresponding  to  a  flow  rate  difference  of  +0.01 
liters/sec  (+0.2  GPM).   The  particle  injection  rate  was 
adjusted  so  that  approximately  10  to  20  particles  visibly 
passed  through  the  observation  section  within  a  10-second 
period.   After  due  preparation,  the  film  was  loaded  into 
the  camera  and  the  shutter  opened  for  approximately  seven 
seconds  to  expose  the  film.   The  7-second  exposure  time 
was  sufficient  to  obtain  an  image  of  the  pipe  without  over- 
exposing the  film.   Images  of  the  pipe  walls  were  possible 
because  of  the  light  being  scattered  from  dirt  particles  on 
the  pipe  wall  and  slight  differences  between  the  refrac- 
tive indices  of  the  TCE,  glass  pipe,  and  the  glycerin  sur- 
rounding the  pipe  (see  Section  3.1.1  for  a  description  of 
the  pipe  enclosure). 
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The  combined  optical  geometry  and  the  photographing 
technique  provided  photographs  from  which  consistent  data 
were  obtained.   Care  was  taken  during  the  construction  of 
the  camera  and  camera  support  to  obtain  a  trueness  of  the 
optical  geometry.   The  angle  between  the  glass  plates  of 
the  pipe  enclosure  was  adjusted  so  that  the  glass  plates 
coincided  with  the  edges  of  a  squaring  tool.   Also  the 
right  angle  of  the  angle-iron  on  which  the  lense  assemblies 
were  mounted  was  checked  for  squareness.   In  support  of  the 
trueness  of  the  optical  geometry,  no  distortions  were  no- 
ticed in  the  images  registered  on  the  film. 

The  photographic  quality  of  the  particle  images  was 
not  as  good  as  was  hoped.   A  fully  automated  tracking  al- 
gorithm using  a  PDP-11  minicomputer  was  tried  but  was 
abandoned  because  the  approach  required  excessive  time. 
Extraneous  images  and  a  varying  background  light  intensity 
required  time-consuming  scans  by  the  algorithm  to  establish 
the  particle  track  position.   Consequently,  an  algorithm 
was  employed  using  the  minicomputer  to  data-log  coordinate 
positions  which  were  located  manually.   The  extraneous 
images  were  present  because  of  scratches  on  the  glass  pipe 
and  dirt  and  dust  particles  which  were  ever  present  in  the 
work  area  and  would  settle  on  the  glass  plates  of  the  en- 
closure.  The  light  intensity  available  for  forming  track 
images  on  the  film  varied  because  of  forward  scattering  of 
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light  from  the  hollow  glass  beads.   Within  any  given  frame 
of  film,  the  particle  images  were  the  darkest  at  positions 
nearest  the  light  source.   For  many  of  the  particle  tracks 
photographed,  the  images  of  the  particle  paths  were  not 
detectable  over  a  sufficient  length  of  the  frames  to  match 
up  tracks  on  the  adjacent  frames.   If  the  above  two  problems 
could  be  eliminated,  a  less  time-consuming  computer  data 
collection  algorithm  could  be  implemented.   The  extraneous 
images  could  be  reduced  by  polishing  the  surface  of  the 
glass  pipe  to  remove  scratches  and  by  moving  the  apparatus 
into  a  more  controlled  environment.   The  particle  images 
can  be  enhanced  in  several  ways;  increased  light  intensity, 
use  of  a  more  reflective  particle,  and  use  of  a  more  harsh 
developer.   The  latter  two  are  probably  the  more  promising. 
A  possible  reflective  particle  which  was  not  tested  is 
titanium  oxide.   Other  Kodak  developers  which  could  be 
tested  are  HC- 110  and  D-19 . 


4  .  3   Data  Collection  and  Analysis 

Reliable  coordinate  data  on  particle  positions  in 
time  were  extracted  from  two  photographed  views  of  the 
flow  loop  using  a  PDP-11  minicomputer.   The  coordinate 
data  collected  from  the  two  views  for  a  given  particle 
track  were  merged  together  to  yield  the  three-dimensional 
Cartesian  coordinate  points  of  the  particle  path.   These 
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coordinate  points  were  then  converted  to  cylindrical 
coordinates  from  which  the  Lagrangian  correlations  were 
generated . 

As  was  mentioned  earlier,  since  the  quality  of  the 
photographic  images  was  not  as  good  as  desired,  a  manual 
data-logging  algorithm  was  used.   The  projected  image  from 
the  photo-enl arger  yielded  a  total  magnification  varying 
from  2.9  to  3.7.   The  objects  which  were  farthest  from  the 
glass-air  interface  were  magnified  the  least.   The  voltage 
from  the  D/A  unit  of  the  PDP-11  minicomputer  ranged  from 
-10  volts  to  +10  volts  with  a  resolution  of  1/1024  volts. 
The  calibration  of  the  X-Y  recorder  was  set  at  2  volts/inch. 
This  combination  of  image  magnification,  volt  resolution, 
and  voltage  scale  provided  a  length  resolution  of  +0.0001 
meter  (+0.005  inch).   However,  due  to  hysteresis  in  the 
recorder  carriage  positioning  mechanism  and  errors  in  manual 
positioning,  the  accuracy  of  the  position  data  was  about 
+0.0002  meter  (+0.01  inch).   A  particle  path's  image  was  a 
dashed  line  in  which  the  dash  was  about  three  times  longer 
than  the  space  between  dashes,  i.e.  about  0.768  fraction 
of  the  track  was  visible.   This  fraction  was  computed  from 
the  blade  area  of  the  light  chopping  disk.   Since  correspond- 
ing views  of  a  given  segment  of  the  flow  tube  were  projected 
simultaneously,  corresponding  points  on  a  given  track  could 
be  measured  at  the  same  time.   All  measured  points  from 
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corresponding  frames,  along  with  datum  points,  were  stored 
on  punched  paper  tape,  and  then  the  data  were  transferred 
to  magnetic  tape  for  use  with  an  IBM/370  computer.   The 
larger  IBM/370  computer  was  used  to  transform  the  two 
2 -dimensional  views  of  the  flow  tube  into  a  3-dimensional 
cylindrical  coordinate  system.   The  coefficients  of  the 
transforming  equations  for  each  frame  were  computed  as 
discussed  in  Section  3.2.1.   These  equations  were  then 
used  to  compute  the  coordinates  of  the  desired  points  on 
the  particle  paths.   The  longitudinal  distance  was  taken 
as  the  average  obtained  from  the  two  views.   The  standard 
deviation  of  the  differences  between  the  corresponding 
longitudinal  distances  computed  from  the  two  views  is 
0.0002  meter  (0.01  inch).   Position  and  time  coordinates 
for  a  section  of  a  typical  particle  path  are  given  in 
Table  4.3.   Also  in  Figure  4.28,  r-6  data  for  all  the 
tracks  are  presented  graphically.   As  can  be  seen  from 
Figure  4.28  even  though  the  particles  are  less  dense  than 
the  tr  i  chloroethyl  ene ,  they  do  not  tend  to  collect  at  the 
top  of  the  pipe. 

The  data  were  smoothed  to  give  a  uniform  time 
period  and  to  generate  velocity  data.   Smoothing  was  ac- 
complished using  a  least-square  fit  of  eight  points  of  a 
particle  path  to  a  fourth-order  polynomial.   The  eight 
points  were  centered  around  the  point  to  be  smoothed. 
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Top  of  Pipe 


Figure  4.28    Projection  of  the  Particle  Paths  onto 
the  r-8  PI ane 
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Table  4.3   Typical  Particle  Path  Coordinates 


Coordinates  of  Particle  Position 


Time  , 
msec 


Radi  al  , 
i  nch 


Tangential  , 
radians 


Axial 
inch 


0.00 

0.407 

4.72 

0.406 

6.16 

0.395 

10.88 

0.390 

12.32 

0.377 

17.04 

0.374 

18.48 

0.363 

23.20 

0.359 

24.64 

0.342 

29.36 

0.337 

30.80 

0.319 

35.52 

0.312 

36.96 

0.290 

41  .68 

0.286 

43.12 

0.283 

47.84 

0.283 

49.28 

0.282 

1  .359 
1  .367 
1  .333 
1  .317 
1  .242 
1  .218 
1  .161 
1  .155 
1  .151 
1  .146 
1.111 
1  .099 
1  .036 
1  .019 
0.983 
0.977 
0.958 


551 
636 
790 
837 
992 
042 
218 
269 
430 
479 
642 
693 
870 
926 
117 
175 
365 
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Since  the  points  of  the  particle  paths  were  taken  as  the 
end  points  of  the  dashed  segments  of  the  particle's  image 
on  the  film,  these  points  were  not  distributed  uniformly 
in  time.   Therefore  during  the  smoothing  step,  the  particle 
positions  were  computed  for  a  uniform  time  period  of  3.08 
milliseconds,  i.e.  half  the  time  period  of  the  dashes.   The 
velocities  of  the  particles  were  also  computed  for  each 
time  using  the  derivative  of  the  fourth-order  polynomial 
generated  from  the  least-square  fitting. 

Using  a  discrete  fast  Fourier  transform  of  the  data, 
the  equivalent  convolution  function,  g(t),  for  the  least- 
square  smoothing  technique  was  examined.   The  effective  con- 
volution function  is  shown  in  Figure  4.29.   The  function 
shown  in  the  figure  was  computed  from  the  data  of  one  track, 
but  it  is  typical  of  the  functions  examined  for  each  of  the 
tracks.   As  would  be  expected,  the  shape  of  the  curve  indi- 
cates that  only  those  points  near  the  point  being  smoothed 
contribute  to  the  smoothing.   Since  eight  time  intervals 
were  used  in  the  smoothing  and  they  were  centered  around 
the  point  being  smoothed,  the  function  is  expected  to  have 
a  zero  value  for  times  greater  than  four  time  intervals. 
From  Figure  4.29  the  function  goes  to  zero  after  three  time 
intervals.   This  indicates  that  the  smoothing  technique 
placed  nearing  all  the  weight  on  only  five  time  intervals, 
the  point  itself  and  two  points  on  each  side  of  it.   Figure 
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Figure    4.29  Equivalent    Convolution    Function    for    the 

Least-Square    Smoothing 
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4.30  shows  the  discrete  Fourier  transform  of  the  function, 
G(f),  where  f  is  the  frequency,   Note  that  the  smoothing 
reduces  the  high  frequency  content  of  the  data.   This  was 
not  considered  detrimental  in  the  analysis  since  the  parti- 
cle path  iamges  were  slow  varying  curves  and  better  estimates 
of  the  particle  velocities  were  obtained  than  by  finite  dif- 
ferencing of  the  experimental  position  data.   Also,  the 
smoothing  enhanced  the  manipulation  of  the  data  by  making 
it  possible  to  compute  the  position  of  the  particles  with 
a  uniform  time  interval. 

The  deviation  of  the  axial  velocity  was  computed  as 
the  difference  between  the  actual  velocity  and  a  mean  velocity 
determined  from  a  mean  velocity  profile  equation.   The  mean 
axial  velocity  profile  was  established  using  an  equation  form 
suggested  by  Pa  i  (1953).   A  least-square  fit  of  the  instan- 
taneous axial  velocities  to  the  following  equation  form 


uz,MAX        '  ro         o 


(4.1) 


yielded  the  following  parameters 

-0.3626 
-0.6374 


1 


75 


z.MAX 


z.AVE 


=  1.20  meters/sec  (3.93  feet/sec 


0.969  meter/sec  (3.18  feet/sec) 


129 


o 
a 


U3  4- 


-0.80 


0.00 


0.80 


1.60 


2.40 


Log  (f) 


Figure  4.30 
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Averaged  velocities  for  50  equally  spaced  radial  intervals 
are  plotted  in  Figure  4,31  along  with  the  fitted  equation. 

From  a  statistical  viewpoint,  the  correlations  ob- 
tained in  this  study  should  be  viewed  cautiously  since  a 
minimal  of  data  were  used  in  generating  the  correlations. 
However,  even  though  only  seven  particle  tracks  were  exa- 
mined, approximately  1600  particle  positions  were  used  in 
the  correlations.   The  number  of  paths  used  in  the  correla- 
tions is  not  the  sole  consideration  in  examining  the  statis- 
tical accuracy  of  the  data.   Jenkins  and  Watts  (1969)  note 
that  the  covariance  function  of  a  stochastic  process  may 
be  estimated  with  arbitrarily  small  error  from  a  single, 
sufficiently  long  record.   For  the  covariance  function, 
the  average  over  one  infinite  record  is  equivalent  to  an 
ensemble  average,  and  hence  the  covariance  function  is  said 
to  be  ergod i  c . 

As  seen  in  Figures  4.22  through  4.24,  the  period 
over  which  the  velocity  deviations  are  correlated  is  rather 
small,  about  3  dimensionless  time  units.   This  time  cor- 
responds to  20  time  intervals  of  the  experiment.   In  a 
rough  statistical  interpretation,  the  1600  particle  posi- 
tions may  be  considered  equivalent  to  80   independent 
particle  paths  having  a  20-ti me- interval  length.   As  the 
path  length  is  increased,  the  number  of  equivalent,  inde- 
pendent paths  decreases.   In  view  of  this,  the  values  of 
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the  correlations  are  most  reliable  near  time  zero  and 
become  less  and  less  credible  as  the  value  of  the  time 
argument  increases.   This  is  also  indicated  by  the  number 
of  contributions  to  the  correlations  for  the  different 
times  ,  see  Tabl e  4.4, 

In  order  to  estimate  the  statistical  error  involved 
in  the  correlations,  consider  the  following  analysis  of  the 
axial  dispersion  correlation.   For  naturally  occurring  ran- 
dom variables,  the  normal  distribution  can  be  used  to  ap- 
proximate the  statistical  distribution  of  their  values 
(Jenkins  and  Watts,  1969).   This  fact  can  be  explained  by 
the  central  limit  theorem,  which  states  that  the  p.d.f.  of 
the  sum  of  N  random  variables  tends  very    quickly  to  the 
normal  form  as  N  increases.   This  theorem  applies  to  any 
form  of  the  p.d.f.  for  the  individual  random  variables. 

Therefore,  assume  the  random  variable  z'(t)  to  be  normally 

2 
distributed  with  zero  mean  and  a   variance.   Note  that 

<z'  (t)>  is  the  sample  variance  of  N  measurements  of  z'(t). 

2      2 
Under  these  conditions,  the  random  variable  N<z'  (t)>/o 

has  a  chi -square  distribution  with  N-l  degrees  of  freedom 


(Hogg  and  Craig,  1965).   For  large  N,  the  distribution  of 
(N<z   (t)>/o  -N)//2N  approaches  a  normal  distribution  wit! 
zero  mean  and  unity  variance.   Now  the  probability  that 


-1.96  <  (N<z,2(t)>/a2  -  N)//2N  <  1.96        (4.2) 
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Table  4.4   Number  of  Contributions  to  the  Correlations 


Time, 

Sampl e 

t* 

Size 

0.145 

3186 

0.726 

3126 

1  .452 

3056 

2.178 

2986 

2.904 

2916 

3.631 

3846 

4.357 

2776 

4.083 

2706 

Time  , 

Sampl e 

t* 

Si  ze 

5.809 

2636 

6.535 

2566 

7.261 

2500 

7.987 

2440 

8.713 

2380 

9.439 

2320 

10.166 

2260 

10.892 

2200 
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is  0.95.   Equation  (4.2)  implies  that 


a2(l  -  2.77//N)  <  <z,2(t)>  <  a2(l  +  2.77//N)   (4.3 


is  the  95%  confidence  interval  for  the  axial  variance.   A  p  - 

? 
proxitnating  o  with  the  variance  of  the  axial  spread,  i.e. 

2 
<z '  (t)>,  an  estimated  statistical  error  can  be  computed. 

The  upper  and  lower  bounds  as  estimated  from  the  above  analy- 

2 
sis  are  shown  in  Figure  4.32  for  < z '  (t)>. 

As  given  in  Section  2.1.1  for  a  time  stationary  pro- 
cess, the  Lagrangian  correlation  tensor  is  symmetric  about 
the  zero  elapsed  time.   Figure  4.33  through  4.36  show  the 
Lagrangian  double  velocity  correlations  for  selected  compo- 
nents of  the  tensor  for  the  fourth  radial  zone  from  the 
center  of  the  pipe.   This  radial  zone  was  selected  because 
it  contained  the  most  contributions  to  its  correlations. 
Symmetry  for  small  elapsed  times  is  expected  because  of  the 
manner  in  which  the  correlations  were  computed  from  the  data 
However,  symmetry  of  the  diagonal  components  is  quite  notice- 
able even  for  large  times.   The  off-diagonal  component 

<V'(t  )V  (t)>,  shown  in  Figure  4.36  does  not  demonstrate 
z  \  0  /  r  v  /  a 

the  expected  symmetry.  If  either  of  the  curves  were  shifted 
vertically,  the  symmetry  would  be  excellent.  The  reason  for 
the  offset  is  not  apparent  from  a  review  of  the  path  data. 

The  power  spectra  were  computed  for  the  three  prin- 
cipal directions  using  a  discrete  fast  Fourier  transform 
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Figure  4.32    An  Estimate  of  the  Statistical  Error 
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Figure    4.36 
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algorithm.   The  spectra  were  computed  using  the  velocity 
data  which  were  obtained  from  the  smoothed  particle  posi- 
tion data.   When  applying  the  Fourier  transform,  the  data 
are  treated  as  though  it  were  from  a  homogeneous  process. 
As  a  consequence,  detection  of  radial  variations  in  the 
spectra  is  not  possible.   The  energy  functions  for  the 
three  coordinate  directions,  F  (f),  F  (f),  and  FQ(f),  are 
shown  in  Figures  4.37  to  4.39  respectively,  where  "f"  sym- 
bolizes the  frequency.   There  is  an  appreciable  amount  of 
variation  in  the  spectra  which  is  probably  because  of  the 
relatively  small  sample  size.   However,  overall  character- 
istics can  be  noted.   From  Figure  4.30,  the  contributions 
from  frequencies  above  100  cps  would  be  drastically  reduced 
by  the  smoothing  technique.   The  Nyquist  frequency  for  this 
data  is  162  cps.   A  close  examination  of  the  power  spectra 
shows  that  at  100  cps,  the  power  function  has  decreased  to 
about  0.01  of  its  original  value.   The  power  is  distributed 
similarly  for  each  of  the  principal  directions.   Decreases 
in  the  power  of  4  orders  of  magnitude  occur  over  a  2  decade 
change  in  frequency. 

As  shown  in  Hinze  (1959),  the  velocity  autocorrela- 
tion functions,  RJF'(t),  Rj.F^(t),  and  Rg   (t)  are  related 
to  the  energy  spectra  functions  F  ,  F  ,  and  F   respectively. 
The  autocorrelation  functions  were  computed  by  taking  the 
inverse  Fourier  transforms  of  the  energy  spectra  functions. 
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Figure   4.37  Axial    Energy    Spectra    Function 
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Figure    4.38  Radial    Energy    Spectra    Function 
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Figure  4.39    Tangential  Energy  Spectra  Function 
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They  are  shown  in  Figures  4.40  through  4.42.   In  comparing 
these  correlations  with  those  computed  with  ensemble  aver- 
ages, Figures  4.22  through  4  .  24,  excel  1 ent  agreement  is  found, 
This  is  expected  since  in  both  methods  the  flow  was  consid- 
ered to  be  homogeneous. 

It  has  been  noted  in  the  results  that  most  of  the 
correlations  involving  the  tangential  direction  have  behaved 
unexpectedly.   These  correlations  are  the  <z  '  ( t )  r  0  '  ( t )  >  , 
<r  6'(t)V^(t)>,  <r  6'(t)V£(t)>,  and  <V^( tQ)V£( t)>  given  in 
Figures  4.8,  4.16,  4.18,  and  4.26  respectively.   The  common 
reason  given  for  this  behavior  was  that  the  sample  size  was 
relatively  small.   Additional  insight  may  be  found  for  this 
behavior  in  the  nature  of  the  tangential  displacement  for 
particle  paths  passing  near  the  center  of  the  pipe.   In 
particular,  consider  the  tangential  displacement  for  a  par- 
ticle path  which  coincides  with  a  diameter  of  the  pipe.   The 
tangential  displacement  will  be  zero  until  the  particle 
passes  over  the  center  of  the  pipe.   After  passing  through 
the  center,  the  displacement  becomes  2irr  ,  where  r   is  the 
starting  radius  for  the  path.   Similar  behavior  in  the  tan- 
gential displacement  would  occur  for  particles  which  pass 
near  the  center.   One  such  particle  track  was  obtained  from 
the  experiment.   If  the  positive  contributions  of  this  path 
to  the  above-mentioned  correlations  is  not  countered  by 
equal  negative  contributions,  the  observed  behavior  could 
occur . 
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Figure    4.40  Axial    Autocorrelation    Function 

Obtained    from    the      Axial    Energy 
Spectra    Function 
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Figure  4.41  Radial  Autocorrelation  Function 
Obtained  from  the  Radial  F.nergy 
Spectra  Fu  nc ti  on 
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Figure  4.42  Tangential  Autocorrelation  Function 
Obtained  from  the  Tangential  Energy 
Spec  tra  Function 
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It  should  also  be  noted  at  this  time  that  the 
smoothed  velocities  for  the  tangential  direction  differed 
significantly  from  the  tangential  velocities  computed  by 
differencing  consecutive  time  intervals  of  the  particle 
tracks.   The  magnitude  of  the  velocities  computed  from  the 
smoothed  data  was   several  times  larger  than  those  veloci- 
ties computed  by  differences.   The  correspondence  between 
the  smoothed  and  measured  angles  was  excellent.   Also 
velocities  computed  by  the  two  methods  for  the  other  two 
directions  were  in  good  agreement.   In  view  of  the  fair 
agreement  between  the  two  methods  for  the  axial  and  radial 
directions,  it  was  decided  to  use  the  smoothing  technique 
to  compute  all  the  velocity  components.   The  larger  magni- 
tude of  the  tangential  velocities  computed  from  the  smoothed 
data  would  be  expected  to  only  increase  the  positive  or  nega- 
tive values  of  the  non-zero  correlations. 

In  this  chapter  the  various  Lagrangian  character- 
istics of  turbulent  pipe  flow  have  been  presented  in  the 
context  of  the  results  of  the  experiment.   The  overall  re- 
sults support  and  enforce  certain  concepts  of  the  Lagrangian 
aspects  of  the  flow  field.   The  photographing  technique  de- 
veloped in  this  investigation  could  be  used  with  some  im- 
provements to  generate  additional  particle  path  data  with 
more  ease.   With  additional  studies  of  this  type,alarger 
sample  size  could  be  accumulated;  and  consequently,  the 
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statistical  errors  and  uncertainties  could  be  reduced. 
Lagrangian  characteristics  for  a  range  of  flow  conditions 
could  also  be  established.   With  minor  changes  in  the 
equipment,  drag  reduction  could  be  studied  too. 


CHAPTER  5 
SUMMARY  AND  CONCLUSIONS 


1.  The  motion  of  suspended  particles  in  fully- 
developed  turbulent  pipe  flow  was  experimentally  observed. 
From  a  close  examination  of  the  images  of  the  particle 
paths  and  visual  observations  of  the  flow,  fully-developed 
turbulence  appeared  to  be  present  in  the  flow  loop  at  a 
Reynolds'  number  of  110,000. 

2.  A  photographing  technique  was  developed  using 
black  background  illumination  which  enabled  the  photograph- 
ing of  the  motion  of  suspended  particles  over  large  dis- 
tances.  Photographic  images  of  particle  paths  were  obtained 
for  tens  of  pipe  diameters.   Even  though  photographic  quality 
was  not  as  good  as  desired,  particle  path  images  were  con- 
trasting enough  to  establish  the  particles  3-d i mens ional 
coordinate  positions  in  time.   The  interrupted  illumination 
of  the  particles  was  an  effective  method  for  embedding  a 
time  scale  into  the  particle  path  images. 

3.  Lagrangian  time  correlations  were  computed  from 
displacement  and  velocity  deviations  obtained  from  the 
particle  trajectory  data.   Displacement-displacement,  dis- 
placement-velocity, and  velocity-velocity  correlations  were 
obtained. 
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4,   The  time  dependence  of  the  various  correlations 
is  apparent  from  close  examination  of  the  plotted  correla- 
tions.  The  effect  of  a  confining  wall  on  the  radial  dis- 
persion was  visible  at  large  times.   Even  though  the  radial 
eddy  diffusivity  coefficient  is  time-dependent,  the  rate  of 
change  of  the  correlation  is  small  after  a  rapid  initial  in- 
crease from  zero.   Consequently,  it  appears  the  radial  dis- 
persion can  be  effectively  modeled  with  a  constant  radial 
diffusivity  coefficient. 


APPENDIX  A 
DATA  COLLECTION  ALGORITHM 

This  appendix  discusses  the  equipment  used  in  ob- 
taining data  from  the  photographs  and  the  algorithm  for 
collecting  the  data.   A  general  purpose  real  time  BASIC, 
CHE  BASIC,  was  developed  and  used  for  controlling  the  data 
collection  equipment.   A  brief  introduction  will  be  given 
later  to  the  computer  language  CHE  BASIC,   The  form  of  the 
experimental  data  was  photographs  of  paths  of  solid  par- 
ticles suspended  in  turbulent  pipe  flow.   The  form  of  the 
data  needed  for  computational  purposes  was  the  coordinate 
positions  of  the  suspended  particles  at  sequential  intervals 
of  time.   The  particle  paths  were  photographed  from  two  per- 
pendicular views  in  order  to  obtain  the  3-dimens ional  coordi- 
nate position,  see  Figure  3.2. 

Data  were  collected  through  the  use  of  three  pieces 
of  equipment.   A  Digital  Equipment  Corp.  PDP-11/20  mini- 
computer was  used  in  conjunction  with  a  Biesler  photo-en- 
larger  and  a  Moseley  model  2D-2  X-Y  recorder.   The  photo- 
enlarger  was  used  to  simultaneously  project  two  photographic 
images  of  the  perpendicular  views.   The  images  were  focused 
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onto  a  plane  surface  of  a  board  mounted  on  the  pen  carriage 
of  the  X-Y  recorder.   A  visible  point  marked  on  the  board 
will  be  called  the  congruent-point  hereafter.   The  congruent- 
point  could  be  moved  at  will  by  varying  the  input  voltages 
to  the  X-Y  recorder,  thereby  making  it  possible  to  match  any 
point  on  a  particle  path  segment. 

In  order  to  control  the  X-Y  recorder  and  to  measure 
the  position  of  the  congruent-point  on  the  pen  board  of  the 
recorder,  an  interface  panel  between  the  minicomputer  and 
recorder  was  constructed.   A  schematic  of  the  interface 
panel  is  givenin  Figure  A.l.   Both  analog  and  digital  sig- 
nals were  used  to  control  the  X-Y  recorder.   The  position  of 
the  congruent-point  could  be  controlled  by  analog  signals, 
either  manually  from  the  interface  panel  or  remotely  from 
the  minicomputer.   With  the  toggle  switch  on  the  panel  posi- 
tioned at  remote,  a  digital  signal  from  the  minicomputer  was 
used  to  select  the  mode  of  operation.   An  "OFF"  signal  estab- 
lished manual  mode  and  an  "ON"  signal  established  remote  mode 
If  the  toggle  switch  was  in  the  manual  position  only  manual 
operation  was  possible.   The  recorder  pen  position,  up  or 
down,  was  also  controlled  by  a  digital  signal  from  the  mini- 
computer, "OFF"  for  pen  up  and  "ON"  for  pen  down.   A  digital 
signal  compatible  with  the  digital  input  channels  of  the 
minicomputer  could  also  be  generated  from  the  interface  panel 

To  ease  the  task  of  controlling  the  X-Y  recorder  with 
the  analog  and  digital  signals,  a  general  purpose  real-time 
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BASIC  computer  language  was  developed.   Six  new  statements 
and  five  new  functions  were  added  to  the  PDP-11  1  of  8  USER 
BASIC  so  that  BASIC  programs  could  be  written  to  input  and 
output  analog  and  digital  signals  to  and  from  the  computer. 
The  additions  were  modeled  after  similar  statements  and 
functions  in  the  Fischer  Controls  pC   language.   The  language 
thus  created  is  referred  to  as  CHE  BASIC.   A  programmer's 
manual  has  been  prepared  for  easy  user  application  of  a  real- 
time computer  system. 

The  algorithm  used  for  collecting  the  data  not  only 
incorporated  the  use  of  the  above-mentioned  equipment  but 
also  a  CHE  BASIC  computer  program  listed  at  the  end  of  this 
appendix.   The  computer  program  used  two  analog  output  chan- 
nels, and  one  digital  input  channel.   This  program  used  only 
part  of  the  total  number  of  channels  available  on  the  PDP-11 
computer.   Initial  attempts  were  made  to  automate  the  data 
collection  algorithm  using  a  photocell  mounted  on  the  pen 
board  to  detect  the  light  intensity  of  the  image.   In  order 
to  establish  the  position  of  the  track,  excessive  scanning 
of  the  image  was  required.   In  view  of  the  relatively  slow 
response  of  the  equipment,  the  collection  of  data  was  very 
slow  using  a  fully  automated  algorithm.   Because  of  the 
slowness  of  the  automated  procedure,  a  more  operator-oriented 
algorithm  was  used  to  collect  the  data.   A  flow  diagram  of 
the  algorithm  used  is  given  in  Figure  A. 2.   A  description  of 
the  a  1 gor i  thm  f ol 1 ows : 
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(1)  The  following  variables  are  read  into  the  program; 

a)  Input/Output  channel  numbers  for  the  analog  and 
digital  signals. 

b)  Correlating  matrix  for  least-square  fitting  of 
data  points. 

c)  Number  sequence  for  order  in  which  the  datum 
points  are  to  be  read . 

(2)  The  program  requests  the  following  input  from  the  key- 
board of  the  teletype; 

a)  Run,  Frame,  and  Track  number 

b)  Status:   0  for  new  track,  1  for  restart  of  old 
trac  k 

(3)  If  starting  a  new  track  go  to  (6). 

(4)  Input  datum  points  and  track  points  from  previous  frame 
Use  the  high  speed  paper  tape  reader  for  input  of  data 
from  previous  frame. 

(5)  Go  to  (7) . 

(6)  Datum  points  are  input  for  current  frame  by  positioning 
the  congruent-point  manually  and  pressing  the  acknowl- 
edge  button.   This  is  done  for  all  14  datum  points 

s  hown  in  Fi  gu  re  3.7. 

(7)  From  the  datum  point  data  align  the  views  and  compute 
the  coefficients  of  Equation  (3.7). 

(8)  If  not  starting  a  new  track  go  to  (12). 
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(9)   Set  switch  register  switch  1  of  PDP-11  console  to  zero 
(down  position).   Position  the  congruent-point  manually 
to  correspond  with  a  point  of  the  new  track  in  top  view. 
When  the  acknowledge  button  is  pressed,  the  coordinates 
of  the  point  are  read  into  the  computer.   The  program 
then  requests  the  corresponding  point  from  the  bottom 
view.   At  this  point  the  program  reads  the  position  and 
computes  the  X  position  to  correspond  with  the  X  posi- 
tion read  for  the  top  view  and  moves  the  carriage  accord- 
ingly.  By  manually  adjusting  the  vertical  position,  all 
possible  vertical  positions  in  the  bottom  view  can  be 
scanned.   The  X  position  can  also  be  adjusted  within  a 
2  inch  span.   When  the  appropriate  point  is  found,  press 
the  acknowledge  button.   This  point  is  read  in  and  the 
program  checks  for  agreement  in  the  X  coordinates  read. 
If  the  X  positions  for  the  points  of  the  two  views  do 
not  agree  within  0.015  inch,  another  point  from  the 
bottom  view  is  requested.   The  first  4  track  points  are 
sequentially  read  into  the  computer  in  this  manner.   The 
remaining  track  points  for  the  frame  are  predicted  from 
the  previous  2  points,  displayed  on  the  X-Y  recorder, 
and  manually  adjusted  to  correspond  to  the  projected 
points  in  the  images.   Pressing  the  acknowledge  button 
initiates  the  reading  of  the  adjusted  point  and  the 
display  of  the  predicted  point  for  the  other  view.   The 
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program  will  display  the  predicted  point  in  alternate 
views  for  a  given  point  until  a  change  in  the  position 
of  switch  register  switch  0  on  the  console  is  observed. 
At  this  time  the  program  predicts  and  displays  the  next 
point  unless  switch  register  switch  1  is  set  to  1  (up 
position)  in  which  case  the  reading  in  of  points  for 
the  current  frame  is  discontinued. 

(10)   Compute  the  coordinates  in  inches  of  the  last  8  points 
of  the  current  frame  to  be  used  to  predict  the  first 
point  of  the  track  in  the  next  frame. 

(11  )   Go  to  (18) 

(12)  If  a  restart  is  in  progress  go  to  (10). 

(13)  Predict  next  point  using  correlating  matrix  and  last 
7  points. 

(14)  Display  point  on  X-Y  recorder  in  top  view.  Position 
switch  register  switch  1  to  down  position.  Manually 
adjust   position  of  congruent-point  to  match  image. 

(15)  Press  acknowledge  button  to  initiate  reading  of  point. 
If  position  of  switch  register  switch  0  has  changed 

go  to  (17) . 

(16)  Display  point  on  X-Y  recorder  in  other  view.   Manually 
adjust   position  of  congruent-point  to  match  image. 

Change  position  of  switch  register  switch  0  if  points 
of  the  views  agree.   If  it  is  the  last  point  of  frame 
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position  switch  register  switch  1  to  1  (up  position). 
Go  to  (15) . 

(17)  If  position  of  switch  register  switch  1  is  down,  go 
to  (13)  . 

(18)  Output  datum  points  and  track  point  data  for  current 
frame  on  the  low  speed  punch.   Go  to  (6). 

The  BASIC  computer  program  of  the  data  collection  algorithm 
is  listed  on  the  following  pages. 
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THIS  P."!a  GRAM  WAS 
DI  SPLAYED  ON  A  X- 
POSITION  OF  THE  : 
SIGNALS  FROM  THE 


written   to    collect   data 
y   recorder:,      the   pen    carriage 
ecorder    is    control   hy    voltage 
pdp-11    minicomputer, 


10  1 
102 
105 
1  06 
107 
1  08 
109 
1  1  5 
1  16 
1  50 
151 
152 
1  53 
154 
155 
1  58 
168 
169 
170 
200 

20  1 
209 

21  0 
21  1 
212 

21  4 
219 
220 
222 
224 
225 
226 
227 

22  3 
230 
231 
232 

23  J 
23  4 
23  6 
236 
240 
25  0 
252 
254 
256 
257 
258 


O/A  CHANNt 
:  *  A  1     A  2 
1  t     2 


L! 


,  Yl 
A  i 
0  , 


AND 
D  I 

0 


21 


u2 
1 


U.J 
0, 


Y  1 

1 


-.01 015699 , 
.59 10  213 
. 01 0  559  05 , 
.  0358977o  . 


0 

.  1  07341  J 
. 03583145 
.0 1 6850  78 
.9 9916  04 


REMARK  * 

REM 

REM 

RE  M 

REM 

REMARK 
DATA 

REMARK  *.**  CORRELATING  MATRIX 
DATA  -.01458767,  .06993961. 
DATA  -.02588002,  .4729543, 
DATA  .0  02071306,  -.01449972, 
OAT  A  -.007939577,  -.0357  0045, 
REMARK  «**  ORDER  OF  DATUM  POINTS 
DATA  2,8,6,5,4,0,1 ,0,1,4,5,0,2,3 

REMARK  ***  DEFINITION  OF  FUNCTIONS  USED   IN  PROGRAM 
DEF  FMa(J)  =  D(J,K  +  I-2)+-D(  J,K+I-l)-D(J.K+I-3) 
FNQ{ K  )=2*K-4* INT  (K/2  )- 1 
PNC ( I ) =F ( L , I ) /h ND ( FNH ( I ) ) 
FND( X)=SQR( 1  . 1 9*X  +  2  0  ) 
FNE( X)=X+2-3* INT ( (X+2J/3) 
FN  I (  I  }=2*INT(  (  1-1 ) /2 )  +  1 

FNP.  (  I  )=(F(K-l  ,  I  )*F(K-l,  I  )+F{K,  I  )  **F(K,  I)  J  *C  (  J,  6) 
FNi(K)=(A(0,K)'A(0.K)  +  A(l,K)»A(l,K))*C(0,6) 
FNT(K)=(  A<2 , K) *M2,K  )  +A(  3,K)  *A(3,K)  )  *Cl  1  ,6) 
A{  3,  1  )  ,d(  1  ,5)  ,C(  1  ,9  )  ,0(3  ,80)  ,00(3  ,80J 
F  (  7 , 6 )  , T (  1,7)  , X (  2  ,  9  )  ,  Y  0  (  J  )  , N (  1  >  ,  E  <  1  , 5  } 
READ  I/O  CHANNEL  NUMBERS 
.  A  3  ,  D  1  ,  D  2  ,  D  3  ,  Y  1  ,  L  1 


DEF 

DE  F 
DEF 
DEF 

DEF 

OFF 

DEF 

DEF 

DI  M 

DIM 

REMARK 

f  '■'  E  A  0     A  1 


REMARK  *  *  * 

FOR  1-0  TO 

FOR  J=0  TO 

.-.  E  M  A  3  K  *  •*  * 

FOR  1=0  TO 

F  0  R  J  -  o  r  o 

REMARK  *** 

REM 

PRINT 

INPUT 


CORRELATING     MATRIX 


:AD     TCI, J):     NEXT     JC      NEXT 
O^DEfi     OF     DATJM    POINTS 


R  E  AD 

1 

7  :  R  E 

READ 

1 

6:  READ  F(I,J):  NEXT  J:   NEXT  I 

REQUEST  RUN,  FRAME, 

AND  TRACK  NUMBER,   AND  STATUS 
•INPUT  FRAME  DATA'.  R  0  ,  F  0  ,  J  0  ,  I  0  "  ; 
<  0  ,  F  0  ,  J  0  ,  I  0 


FROM  X-Y  RECORDED 


print  :  if  10  go  to  700 
remark  ***  input  frame  datum  points 
print  "datum  pc  i  nts  for  frame'sfoj 
print  ...  run'»;ro;m.  track";  jo;  ".•• 
contact  di  to  0:  contact  c2  to  1 
for  j=2  to  0  step  -2 

FOR  K=0  TO  t>:      GOSUB  1100 

let  f(  j,e<  j/2  ,k)  )-z".    let  f  (  j  +  1  ,  e  <  j/2  ,  k  )  )  =  y  :  g0su8  1005 
next  k:  print:  print:  next  j 

remark  ***  align  the  views,  i.e.  rotate  axis  and  center 
for  j=0  to  2  step  2:  let  yo{j)=0:  let  y0( j+1)=0 

LET  70  =  ATN((F(Jt-l,3)-F(J  +  l,0))/(F(J,J)-F(J,0)  )) 

LE  T  C( J/2,7) =ATN(  (  r<  J  +  1  ,2  )-F (  J +  1  ,  1  )  )  /  (F  {  J  ,  2  J-F  (  J  ,  1  )  )  ) 

LET  C( J/2 ,7)=< 20+C< J/2 ,7 ) )/2 

FOR  K=0  TO  6:  LET  Sl=-C(J/2,7) 
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260 
262 
268 
270 
272 
274 

27  a 

273 
2  79 
280 
232 

28  3 
28  4 
286 
28  3 
290 
291 
29^ 
294 
300 

30  2 
304 
305 
306 
308 

31  0 
312 
314 
31  5 
316 

31  7 
318 
320 
321 
322 

32  4 
326 
328 
330 
332 
334 
336 
337 
336 

33  ) 
340 

34  6 
348 

35  0 
352 
354 
356 
358 


LET  Z=F(J,K>:  i_ET  Y=F(J+1,K):  GOSUB  1200 

LET  F(J+4.K)=Z:  LET  F(J+5.Ki=Y:   IF  K>3     G 

LET  Y0(  J)  =Y0(  J)  +  F(  J  +  4  ,K)  :  LET  Y0(J*1)=Y0 

NEXT  K 

LET  Y0(  J  ) -Y0(  J)/4  :  LET  Y  0  (  J*  1  )  -Y  0  (  J  +  1  )  /  4 

FOR  K=0  TO  6:  LET  F  (  J+4  »  K  )  =F  (  J+4  *  K  )  -  Y  0  (  J 

LET  F( J  +  5  ,K)=F( J+5,K)-Y0( J+l  ):  NEXT  k:  N 

FEMAiK     ***     COMPUTE     THE     COEFFICIENTS     OF     T 

REM  TRANSFORMATION    EQUATIONS* 

FOR  J=0     TO     i:     LET     K=2*3+5:     LET     L-K 

LET  C(J.0)=(F(K.0)+F(K,3)-F(K,l>-F(K,2)   ) 

LET  C(J,6)=l/C(JtO) | 2 

LET  C(  J  ,  1  >  =  ( F N  C (  0 ) +  FNC (  3 } -FN  C (  1 ) - F NC ( 2)  ) 

LET  C(J,4)=(2«J-1)+F<K,6):     LhT    C<J,5)-(2 

LET  L  =  K-1  :     LET     C ( J , 2  )  =F ( L , 6 ) -F ( L , 4  J 

LET  C( J ,3)=FNC(6)-FNC(4 ) :     NEXT     J 

F  0  R  J  =  0     TO     1  :     LET     B<J,2)=l/C(J»l)-l/C<J» 

LET  d( J ,0) =C( J,0)/C( J, 1  )-C( J  ,2)/C( J , j)  +  C 

LET  S(J,1)=C(J,4)-C(J,2)*C(J.5)/C(J,J): 

FOR  J=0     TO     i:     LET     K=2*J+4 

LET  3  (  J  ,  3  )  =  d  (  0  ,  0  )  *  L  (  1  »Q)-B<J.l)*B(J.2)  +  0 

LET  B(J,4)=B(J,3)*B(J,3)+4*B(Q,0)*D( 1 ,0) 

LET  8{ J ,4 >  =  SQ^ ( 8( J ,4  )  ) 

LET  C  (  J  ,  2  )  =  <  B  (  J  ,  4  )  -  0  (  J  ,  3  )  )  /  (  2*  0  (  J  ,  0  )  *  S  (  J 

L^  T  C(J»4)=B(J,1)*C<J.2)+C{J.5)/C(J,3) 

LET  C(J.3)=(1/C<J.2)-C(J.0))/C(J,1 ) 

LET  C( J , 0 )  =  ( 2  +  J- I) *(C( J ,4)*F (K,6)/F(K+1  , 

LCT  C(  J ,6  )= 1 0*C( J .2  )/C( 0  ,2) :      LET     K=2*J 

LET  C(J,3)  =  (F(K*-1  , 0)-F( K+l  il  )-F(K»6)+F(K 

LET  C(Ji9)=C(J.8)/(F(3-K»0)tF( 3-K, 1 )  -2  +  F 

LpT  C(J,3)=F(K+1  ,  0)-F(K-H  » 1 i -{ F ( 3-K . 0 ) +F 

LET  C( J. 8  1  =  2*C( J .3)  :     NEXT     J 

LET  Z2  =  INT( 1  .5*F0  +  1  .25) :     LET     00=10:      IF     I 

REMARK     **«      INPLT     T"ACK     POINTS     FOR     1ST     FR 

PRINT     "INPUT      1ST     4     POINTS":      LET     Z3-F:     LE 

FOR  K=0     TO     2     STEP     2".     LET     J2  =  0 

FOR  J=0     TO     i:     LET     1=131     LET     Y=Y2:     GOSUB 

LET  M=SWH(1):     LET     L=2*M:     LET     Z3=Z:     LET     Y 

LET  D(L«K+J)=Z:     LET     D { L+ 1  ,  K+ J ) = Y 

LET  Y=Y3:     GOSUB     1210 

LET  D( 2-L  .K  +  J )=Z :     LET     0 ( 3-L .  K  +  J  )  =  Y 

LET  X0  =  A(L,J):     LET     X1=A(3-L.J) 

LET  X2  =  A(  L+- 1  .  J  J  :      GOSUB     1150 

LET  X0=A(L,J)-A(2-L,J):      PRINT     XO, 

IF     ABS( XO ) >0. 015    GOTO     3  32 

PRINT:      NEXT     J:      NEXT     K.  :      PRINT 

P=>INT     "TRACK     END     POINTS":     CEMRK     ***     NOW 


D  TO  2  7  0 

(  J  +  l  )+F(  J+5 , K) 


) 

EXT  J 

HE  COORDINATE 


/  2 


/  2 

*  J- I ) *FNC(6) 


3) 

( 1- J ,5) /C( 1- J ,3) 

N  EXT  J 

( 1-J , 1) *B( 1-J ,2) 
*B(J,1)*S(J,2) 

.  1  )  ) 


b  )  ) 

,4)  ) 

( 3-K, 4) ) 

( 3-K, 1 } )*C{ J, 3) 

0  GO  TO  390 

A  ME 

T  Y2  =  0:  LET  Y3=0 

1010 
2  =  Y 


LET  Y3=Y:  GOSUB  1060 


GET 'EM  ALL 


POINT:   PRINT  " 
FOR  K-0  TO  30 
FO^  1  =  0  TO  1  : 
GOSUB  12tj2 
GO  TO  36  0 
FOP  J-0     TO  3: 


UPPEF  FRAME 
5TE=>  2 
IF  K  GE 


•LOWER  FRAME' 


GO  TO  35  8 


LET  D0(  J  .  0  )=FNA(  J)  :  NEXT  j:  GOSUB  I 
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360 
362 
364 
365 

36  5 
370 
371 
372 

37  4 
376 
378 

38  0 
390 
400 
402 
404 
41  0 
412 

41  4 
420 
422 
426 
427 

42  9 
42  9 
430 
431 
432 
433 
434 
440 
442 
444 
466 
487 
483 
490 
5  0  0 
501 
502 
5  03 
505 
51  0 
51  4 

51  3 
516 

52  0 
523 
524 
525 
52  9 
530 
535 


NEXT  I  :  GOSUG 
FO°  1  =  0  TO  i : 
IF  NOT ( SWR( 2) 


1  050 

PRINT  A(  0 ,  I  )  ; A ( 2,  I  )  .A(  1  .  I  ) 

AND  (K  LT  6))  GOTO  366 


A  (  3  ,  I  ) 


NEXT  I 


print  CH:n(7);:  wait  1000 

IF  SWR(2)=0  THEN  NEXT  K 
REMARK  ***  COMPUTE  COORDS 


RE  M 
LET 
FOR 
LET 

LET 


OF  FRAME 
K0=FNI(K+1):   LET  J2  =  0 
K=K0-7  TO  KO  STEP  2*. 


:    GOTO    364 

IN     INCHES 

JF     LAST     8 

:    let    J l  =  4  : 

LET     J3  =  3 

GQSUB     10  50 

POINT: 


J2  + 


I  F 


J 
J2=6 : 

L'T  K=KO: 
LET  K = K 0 : 
P  E  M  A  <  K  *■  *  *  N  C 
PRINT  "TRACK 
PR  I  NT :  PRINT 
FOR  K-0     TO     JO 


NEXT  K 

10  GO  TO  57  0 
GO  TO  5  00 
IF  I0>0  bO  TO  372 

A  FIND  POINTS  ON  FRAf 

END  POINTS" 

"UPPER  FRAME"  ,  ,  "LOWE." 
STEP  2 


REMARK  ***  PREDICT  NEXT  POINT 

LET  J1=FNF(J1):  LET  J2=FNE(J2):  LET  J3=FNE(J3) 
1  =  0  TO  2:  LET  XO  =  X(  I  .  J  1+  1  )  -X(  I  » J3  ) 

LET  X(  I  ,  J2  *-J  )  =  X(  I  ,  Jl  +  J  )+X0:  NEXT  J: 
LET  X0=0:  LET  Xl=0 
X  3  =  I  N  T  (  X  ) 
LET  J  4  =f:HE  I   J2+-  J) 
(  X(  0 , J4  )-X3  J 
X(1.J4):  LET  X2-X^fT(I,J) 


j=o  to  l : 
i=o  to  l : 

X2=0:  LET 
J  =  0  TO  7 : 
X  0  =  X  0  *■  T  {  I  ,  J  ) 
XI =X1 -T ( I , J ) 


X(2, J4) 


COORDS  AND  GET  POINT 


FO- 

FO3  J=0  TO  i:  LET  X(  I  , J2 *J ) =  X(  I  , Jl + J  )+X0:  NEXT  J:  NEXT  3 
FOR 
LET 
FOR 
LET 
LET 

N  E  X  T  J 

LFT  X0=X0+X3 

REMARK  *+•*  CONVERT  TO  RECORDER 
FOR  J  =  0  TO  i:  GO SOD  115  3 

LET  A4=M0:  LET  Z0=X2:  LET  X2=-Xi:  LET  X1=Z0 
LET  D0(2*J»0)=Z:  LET  D  0  (  2*  Jl-  1  •  0  )  =Y  :  NEXT  J:  G05UO  1250 
NEXT   I  I  GOSUB  1050 

REMARK  ***  LIST  THE  PAIR  OF  POINTS 

FOR  1  =  0  TO  i:  PRINT  A ( 0  i  I )  *A ( 2 , I )  . A ( 1  , I ) , A (3 , I  J :  NEXT  I 
IF  SWR(2)=0  THEN  NEXT  K 
LET  K=K+1 
REMARK  **••*  OUTPUT  DATA:   FIRST  LIST  THE  "C'«  MATRIX 

print:  print:  print  "C  matrix" 

for  1=0  tc  i:  print:  print  i; 

for  j=0  to  7:  print  c  (  i  ♦  j  )  ,  r  next  j:  next  i 

print:  print:  print  "prepare  for  output.";:  gosub  iiro 

LET  C 0-255:  G0SU3  1000 

rfmark    +**    punch   title    of    data    block 

print   "run   nom;ro;"«    frame   no";fo;".    track   no"  ;  jo  ;  •• .  " 

GOSUB  100  0 

LET  C0=0:  GO SUE  1000 

REMARK  ***  PUNCH  RUN,  F^AME,  AND  TRACK  NUMBER  AND 

REM         NUMBER  OF  POINTS  TAKEN  FROM  FRAME 

print  -  o ;  •• ,  "  ;  f  o ;  " .  ■■ ;  J o ;  *• , " ;  k  ;  chr  s  (  i  j  i  ; 

REMARK  **«  PUNCH  DATUM  POINTS 

FOR  J=0  TO  3 

FOR  1=0  TO  6:  PRINT  F ( J , I ) ; C HRi ( 1 3 ) I :  NEXT  I:   NEXT  J 
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540  LET     I l=SGN(D{ 0,0)-C( 0*KJ ) :     LET     K0=0:     LET     K1=K 

54  2  IF     I  KO     THEN     LET     K0=K:     LET    Kl-0 

54'4  REMARK     **«     PUNCH    TRACK    POINTS 

545  FOR     J=0     TO      i 

550  FOR  I  =  K0  TO  Kl  STE=  II  :  PPINT  I  NT ( D ( J , I  )  +  0 . 5 )  ; CHRS(  1 J )  ; 

560  NEXT  i:  NEXT  j:  GOSUB  10  00 

56  5  GOSU-i  1190 

570  print:  p-^int:  print:  let  fo  =  eo+ii:  let  io=-i:  go  to  ejo 

700  remark  *  *  *  restart  with  last  frame 

701  LET  ri=SGN(I0J:  LET  K0  =  0 

702  OPEN  "PR:"  AS  FILE  2".      INPUT  #2  .R  1  .  F  1  •  Jl  ,  K  1 

703  IF  (  <1  NE  R0)  GOTO  706 

7  0  '-*■  IF  IFO  EG  (Fl+Il)  AND  (Jl  EQ  JO)  GO  TO  710 

706  PPINT  :  PRINT  "E..AME  DATA  DO  NOT  MATCH."*.  CLOSE  2:   STOP 

710  FOR  J-0  Tn  3 

712  ro;;  <=o  to  c:  input  #2,fcj.kj:  next  k:  NexT  j 

'14  IF      IKO     THEN     LET     K0=Kl:      LET      K1=0 

72  0  FO-     J=0     TO     J 

722  FO"  K  =  K0  TO  Kl  STEP  I  1  :  I  NPUT  #2»D{J,K):   MEXT  K:  NEXT  J 

72  4  IF  K0=0  THEN  LET  K0=K1 

726  CLOSE  2:  LET  F0=Fl:  LET   10=1 :   GO  TO  250 
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1000 

REMARK  *** 

PR  INT 

1001 

FOR  1  =  1   TO 

20 :  p 

2  0     CHRS(  CO  )  •  S 

INT    CHRKCO)  ;; 


NEXT     I 


:tupn 


1004  REMARK     ***     LIST     AND     MARK    RECORDER     COORD    POINT 

10  0=  PTINT    "     (";  INK  /.);",»;  I  mt(Y  )  ',")"  , 

1006  CONTACT  Dl   TO  l:  WAIT  200:   CONTACT  Dl  TO  0:  WAIT  200 

1007  IF  2*INT( (K+l )/2) =  K+1  THEN  PRINT 
10  09  RETURN 


1010  REM^K  ***  -FAD  RECORDER  COCR0  POINTS  -  FINE 

1012    5lt    z\    to    z+ain(  ai  ) /do:     SET    Yl     TO    Y  +  AIN(A2)/B0 

10  14  IF  ABS  (  AIN(  Al  )  )>1  00  0  THEN  LET  Z-Z  +  5  *SGN  (  A  I  N  (  A  1  )  J 

1016  IF  AdS< AIN( A2) ) >1000  THEN  LET  Y=Y+5* SG N ( AI N ( A2 J ) 

10  16   IF  D I N 10  3) -0     GO  TO  1012 

1020  IF  DIN (03) =1  GO  TO  1020 

10  22  LFT  Z  =  Z+AIN(A1)/B0:  LET  Y  =Y  +  A  I  N(  A2  )  /GO  :   RETURN 


1050  REMARK***  CONVERT  2  POINTS  AT  A  TIME 
10:52  FOR  J=0  TO  i:  GOSUb  1060 
1 054  NEXT  J :   PE TU-N 


10(50 

REMARK  ♦ 

1  062 

FOP 

1=0 

1  0t>4 

LET 

Z  =  D( 

1066 

LET 

A(  I  . 

10  7  0 

LET 

XI  =F 

1072 

LET 

ZG  =  1 

1  074 

LET 

A4  =  ( 

1076 

LET 

A(li 

1078 

FOR 

I  =2 

1  080 

LET 

A<  I  . 

1082 

LET 

A  (  0  , 

10  34 

LFT 

A(  0. 

10  J  5 

LF.T 

A(  2. 

10  36 

LET 

X(  1  , 

1  088 

LET 

X(  Of 

*■*  CONVERT  COORDS  TO  INCHES 

TO  2  STE"3  2:  LET  SI  =-C(  1/2,7) 

Y-O (  I +1  ,K  +  J)  :     GOSUB     120  J 
LET     A{  I  +1  ,  J  >=Y-Y0<   1+  1  )  :     NEXT 
LET     X2=FND(FNT { J ) ) 
2)*A(1,J)*A(3,J)/(X1«X2) 


I  ,  K  +  J )  :     LET 

j) =  z-yo (  I  )  : 

ND(FNS<  J)  )  : 
+  C ( 0,2  )  *C(  1 


C(  0,3  >-C(  1  ,2)*A(3,J)*<1+C(1,3)/X2)  )/Xl 
J)=-A(1.J)*C(0,2)*(  1  +  A4J/Z0 

TO     3 
J>=C(l»2)*A(I,J)*(l+(C(l,3)-A(l,J))/X2) 

J)=C(0,2)*A(0.J)*(1+<C<0,3)-A<3,J) )/Xl ) 

J )  =  Z2-A(0 , J )+C( 0,0) 

J) =Z2-A<2 , J)+C( 1,0) 

J+J2)=A(1,J):     LifT     X(2,J  +  J2)=A(3,J) 

J+J2)  =  (A(0,J)*-A(2,J))/2:     RETURN 


NEXT     I 
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1100  REMARK  * 

110  2  P 3  INT  K ; 

1110  LET  I-   AIN<  Al J  :  LET 

1115  LET  Y=AIN(A2):  LET 

1120  IF  9IN<03)=0  GO  TU 

1  1  JO  IF  DI  N  (  03)  =1  GO  TU 

1140  RETURN 


READ  RECORDER  COORD  POINTS  -  COAP.Si 


Z=AIN( Al ) : 
¥=  A  I  N(  A2  )  : 
1110 
1  1  3  .) 


SET 
SET 


Zl 
Yl 


TO 

TO 


I  ISO  REMARK  ***  CONVERT  FOR  RECORDER  COORDS 

UE2  LET  Z0=Z2-X0+C( J, 0 ) :  _ZT  S1=C{J,71 

1  1  54  LET  MO^C(  J  ,2 >*(  1 +(  C ( J  , J  )-X2  )/END(  ZO* ZO  +  X  i*Xl  )  ) 

1156  LET  Z=Z 0/MO+YOC 2 *J J :  LET  Y  = X 1 / M 0 + Y 0 ( 2 * J+ 1 ) :  GO  SO 5  1200 

1158  RETURN 


11^0  REMARK     ***     WAIT  FOP     ACKNOWLEDGEMENT 

1192  IE     OIN<D3)=0     GO  TO     1192 

1194  IF     DIN(D3)=1      GO  TG      1194 

1196  RETJRN 


1200 
1201 
12  02 
1203 
1204 
1206 


REMARK     **     ROTATE     THE     COORD.     AXIS 

LET     K 1=SQR ( Z*Z+Y"Y ) 

IF  Z-0  THEN  LET  A4 =3 • 1 4 1S9* SGN ( Y J +S I :  GJ  TO  12  06 

LET  A4=ATN(Y/Z) 

LET  A4  =  A4-H  .57079*  SGN  {  A4)*(  SGN( Z )-l ) +S  1 

LET  Z=R l*COS( A4)  :  LET  Y =P 1  * S I N ( A4 )  :   RETJRN 


1210  REMARK  ***  RF.AD  ONE  OF  FIRST  FOUR  FOR  1ST  FRAME 

1212  LET  X2  =  C(M,3)*D(2¥M  +  1  .  K+J  )  «-C(  M  •  9) 

1214  LET  X0  =  Y  +  AIN(  A2  )/B0:  SET  Yl   TO  XO 

1216  LET  X1=C( 1-M,ai *XQ+C<1-M»9) 

1218  LET  Z=(D(2*M, K+ J)-T C2*M ,6) ) *X2/X1+F( 2-2+M, 6 ) 

1220  Ip  ArfS< AIN( A2) ) >1000  THEN  LtT  Y  =Y+5*SG  N<  Al  N(  A2  )  J 

1222  SET  Zl  TO  Z+ A  I N <  A 1  ) /0 0 :   IF  DIN(D3)=0  GO  TO  1214 

1224  GO  TO  10  2  0 

1250  REMARK     ***     GET     NEXT    PJINT 

1252  LET     J  =  0:     LET     Z0=SWP(0) 

1254  LET     Z=00(J.O):     LET     Y=D0(J+1,0) 

1256  SET     Zl      TO     Z:      SET     Yl     TO     Y:     LET     i30  =  C  {  J/2  ,  5  )  :      GOSUO     1010 

12  58  LET     D(J,K+I)=Z:     LET     D(J+1,K  +  I)  =  Y 

1259  IF  SW--(O)  NE  ZO  GO  TO  1262 

1260  LET  J=2-J:   GO  TO  1254 

1262  FOR  J  =  0  TO  2  STEP  2:  PRINT  K+  I  J 

1264  LET  Z=D(J.K+I>:  LET  Y=D(J+1.K+I) 

1266  SET     Zl      TO     Z:      SET     Yl     TO     Y."     WAIT     500:      GOSJO     1005 

12b8  NEXT      J:      PRINT:     'ETU^N 


APPENDIX  B 

SOLUTION  FOR  COEFFICIENT  TO  X-Y  RECORDER 
COORDINATE  TRANSFORMATION  EQUATION 


The  coordinates  as  read  from  the  plane  of  the  X-Y 
recorder  must  be  transformed  to  the  Cartesian  coord i nates  of 
the  flow  tube  in  order  to  fix  the  position  of  a  particle. 
The  equation  for  such  transformation  is  by  Equation  (3.8) 
of  Chapter  3  as 


C  .    .         M  •  v  ■ 
J  m  k  j    j 


d  . 
J 


J  2k 


ii    2    "     i \p2       ,       2M2    2*1/2 
((n      -    l)Cjk    +   n   M.v.) 


(B.l) 


where  the  nomenclature  is  given  in  said  chapter.   Equation 
(B.l)  is  rearranged  to  give 


D  .  .     u  • 

jmk  _   Hj 

C  .  ,    M  ■  v  • 
jmk    j  j 


(1  + 


d  .  -  D^0I 
J     J2k 


M,v-'   jk     Hj' 


).  (B.2) 


((n1 


yi 
The  values  of  ..  J   and  d.  for  each  view  of  a  tube  section 

MJVJ      J 
are  needed  to  compute  Djik  and  D j 2 k '  Cartesian  coordinates 

within  the  flow  tube.   To  accomplish  this,  first  express 

the  known  distances  between  certain  datum  points  in  terms 
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of  these  unknown  parameters,   Since  the  inside  diameter  of 
the  tube  is  one  inch,  then  for  points  (1)  and  (2)  on  Figure 
3.7, 


CT21mi(1  +  X^  "  Ci22'ni(1  +    l-7-z]    =    ' 


(B.3) 


and    for    points    (3)    and    (4) 


ci24"i«1  +  iqj'  "  ci23mi<1  +  rrj»  =  1  (B-4 


where 


(B.5) 


A.  . 
1  J 


(n1 


l)m2C2.    +   n2,2)1/2 


(B.6) 


Note    that 


Di21 


D  i  22         Di  23 


Di24    =    ° 


for    i    =    1,2.      Averaging    these    two    equations    and    rearranging 
the    resulting    equation    yields 


,r  r  c  r       >mi        ,Ci21       Ci24       Ci23       Ci22^midi 


i 21  i24         i23         i22'  2        U^         Ai4         Ai2  12 


rr'-?" 


=i. 


B.8) 
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Similarly,  the  distance  between  inch  markers  on  the  tape 
is  one  inch,  so  that  for  points  (5)  and  (6) 


Cil6    "    Cil5)mi    +    (   A, 


C  .  C  . 

i  1  6  i  1  5  \       /  ,  n         ^    _    i 

—7   "    TT^K^i    -    Df25)    "    ] 
i  6  1  5 


;b.9) 


A  third  equation  can  be  written  to  relate  D-^c  to  d.  and 


m  .  as 
i 


C 
(-l)1C125mi  +  (-i)l-J^m1(di   D,0,) 


i5 


i25 


D  . 


i  2  5 


(B.10 


The  factor  (-1)   is  required  because  view  1  is  inverted. 
Each  of  these  equations  apply  to  both  views,  so  there  are 
six  equations  and  six  unknowns.   These  equations  are  non- 
linear since  A.,  is  a  function  of  m  •  .   Before  attempting 
a  solution,  write  the  equations  in  a  more  compact  form, 
for  i  =  1,2 


where  A 


A . ,m .  +  A.0m.d.  =  1 
ill     i 2  i  i 


Ai3mi  +  Ai4mi(di  "  Di25}  =  ] 


A-r-m.  +  A.cm.(d.  -  Do,c)  =  D.QI- 
i  5  l     n  6  i  v  i     l  25  '     i 25 


^  C  i  3 1  +  C  i  24  "  C  i  23   C  i  2 2  ^ 

n  "         2 

A    =  1 ( C  i  2 1  +  Ci24  _  C  i  23  _  C  i  2  2 

i  2    2  A  .  1     A  .  4     A  .  3     A  .  2 


B.ll  ) 
B.12) 
B.13) 

B.14) 
B.15) 
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Ai3   Cil6  "  Cil5 


(B.16) 


r     c 
Lil6    i!5 


'i4    Ai6     Ai5 


B.17 


Ai5  =  M)1^ 


B.18) 


Ai6    t  D   Ai5 


(B.19) 


Because  the  equations  are  nonlinear,  an  iterative  method  is 
employed.   By  assuming  a  value  for  m-  the  function  X..  can 
be  evaluated  and  by  assuming  a  value  for  D.pc»  Equations 
(B.ll)  and  (B.12)  can  be  solved  for  m.  and  d..  .   Equation 
(B.13)  yields  a  new  approximation  for  D^c  and  ^22  5"   This 
procedure  can  be  repeated  until  values  of  m.  ,  d.,  and  D^5 
have  converged.   Once  m .  and  d.  are  known  for  a  given  tube 
section,  all  points  photographed  within  that  section  of  tube 
can  be  transformed  to  their  corresponding  Cartesian  coordi- 
nates of  the  flow  tube. 


APPENDIX  C 

COMPUTER  PROGRAM  FOR  COMPUTING 
THE  CYLINDRICAL  COORDINATES  OF  PARTICLE  PATHS 


The  theory  and  use  of  this  program  is  found  in 
Section  3.2.1  and  Appendix  B.   Documentation  is  also 
included  within  the  program. 
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THIS  PROGRAM  COMPUTES 
PARTICLE  POSITIONS  FRO 
DEFINITION  OF  MAJOR  VA 


T  J 



T  I  ME     P 

NUMT  3K 



EXP  EH  I 

NUMKri't 



NUMBER 

INITF 



F  I  3  S  T 
VIEWED 

NUMRUN 



LXPE-(  I 

NCF^M 



FF  AME 

N  0  T  F  K 



fcXPER I 

NUMSG1 



NUMBER 
V  I  E  W  E  D 

F 



ARRAY 

C 



AFFAY 

T  ION     E 

0 



A  f  R  A  Y 
TWO     2  - 

X 



ARRAY 
D I  NATE 

V 



AFFAY 

V  Z  A  V  G 



AhRAY 
FIVE     3 

N 

—  ~  — 

AFF  AY 
IN    THE 

THE  CYLINDRICAL  COORDINATES  OF 
M  THE  Tw/G  2-D  JMENSI  ONAL  VIEWS. 
R1ADLES   IN  THE  PROGRAM  FOLLOW: 
EFIOO  OF  TACK,  SECOND. 
MENTAL  TRACK  NJMJER. 

OF  PHOTO  FRAMES  WITH  DATA. 
FPAME  UN  WHICH  PARTICLE  PATH  IS 

WENT  RUN  NUMBER. 

NUMBER  . 

MENT  TRACK  NUM3E3 . 

DF  SEGMENTS  OF  =>  ARTICLE  PATH 

IN  FRAME. 
OF  DATUM  POINTS  FOR  FRAME. 
OF  COEFFICIENTS  = OR  TRANSFJRMA- 
QUAT ION. 

OF  PARTICLE  PATH  DATA  FOR  THE 
DIMENSIONAL  VIEWS  OF  THE  FRAME. 
OF  COMPUTED  CYLINDRICAL  COOK- 
S. 

OF  COMPUTED  VELOCITIES. 
OF  AVERAGE  VELOCITIES  FOR  THE 
ADIAL  ZONES. 
OF  THE  TOTAL  NJMBER  OF  SEGMENTS 

PARTICLE  TRACKS. 


DI  MENS  I  ON  A{ 4  ,2)  ,  C ( 2 .  I  0 )  , D ( 4  ,  2  6 )  , F ( 8 , 7 )  ,ri(  2, 6) 

,X(7,2,3.210)   ,V(  7, 2,  J,?.  10)  , XAVG(4 )  ,XNF(7  ) 
.  V<£AVG<2  .5)  »N  (  10) 

DATA     TO/0. 0  061 6/, VZAVG/1 0*0.0/ 

M  1  =  1 

READ  IN  TRACK  DATA 


R  E  4  D  {  1  )  NU  MTK  K  ,  NU  MF  R M  ,  INITF 

M  =  0 

PEAD(  1  )  NUM  HUN  ,  M OFF M  ,  NOTPK  .NUM3G1  ,F  ,( 

READ( 1 >  D 

ALIGN  FPAMES 


DO  11   1=1,4,2 
XAVG( I )=0 
XAVG( I +1 ) =0 

C(  (  1  +  1  )/2,7)  =  (  <\TAN2  (  F(  1+1 , 4)-F(  1+1  ,1  )  .  F  (  I  .4  )-F(  I,  1  )) 
+  ATAN2( F(  1  +  1  »3)-F ( 1  +  1 ,2) «r(  t  »  3)-F(  I  ,2) ) )/2 
DO  10  J  =  1  ,  7 
CALL  RQTAXIFII  ,  J)  ,F(  1  +  1  ,  J)  ,-C(  { 1  +  1  ) /2 .7 J  , 

F(  1+4 , J  )  ,F(  1+5, J  )  ) 
IF  (J  . GT.  4)  GOTO  10 
XAVG(  I  )=XAVG(  I  )  +F(  1+4, J )/4 
XAVG( I  +  l ) =XAVG< l  +  l )  +F ( 1+5, J  )  /  4 
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10  CONTINUE 

DO  11  J  =  l  ,  7 


DO  11  J  =  l  ,  7 
F{  1+4 , J )=F(  1+4, J)-XAVGt I  ) 
11  F(  1+5,  J  )  =  F(  1+5,  JJ-XAVGU  + 


1  ) 


DO     1 

K  =  2  - 

C(  I  , 

XNr  ( 

CO  , 

KI, 

1  2 

CONT 

3  1 

00     1 

K  =  2* 

:(  i  , 

l  3 


30 
20 


C 


COMPUTE  COFFF  FOR  EQUATION  TO  ADJUST  CHORDS 

2  1=1,2 
(  I  -  1  )  +  6 

l)=(F(K,t)+E(K,4)-F(K,3)-F(K,2))/2 

I  )  =  1  /  C  (  I  ,  1  ) 

5  )= (2*1-3 )*F(K, 5) 

3)=  (2*1-3)  *F  [  K,5)/(F(Kt-l  ,7)-F  (K-l  ,  5)  } 

INUE 

J  1=1,2 

(  I  -  1  )  +  6 

2 ) = ( FUNC ( F ( K , i ) , F ( K-l , 1 ) , XNF ( I ) ) 
+FUNC(F  (K, 4)  ,F(  K-l  , 4  )  ,  XNF (  I  )  ) 
-FUNC( F(K , J)  ,F (  K-l  , 3)  , XNF {  I )  ) 
-FUNC(F (K, 2)  ,F(  K-l ,2)  , XNF (  I  )  )  )/2 
C( I ,4} =FUNC(F( K-l ,7) ,F (K.7 ) ,  XNF( I) ) 
-FUNC(F( K- 1  ,5 )  ,F  ( K,5  )  , XNF (  I  )  ) 
C(  I  ,6  )  =  (2M-3  )  *FUNC(  F(  K.5)  ,F(  K-  1  ,5)  ,XN=  (  I  )  ) 
C(  I  ,.5)  =F(  K-l  ,7  )-F  (K-l  ,  5)-C  (  I  ,4)  *B(3-I  ,  J  ) 
CONT  INUi: 
DO     14      1=1,2 

A  1  =C  (  I  ,  1  >  *C(  I  ,4) -C(  I  ,2)*C(   1,3) 
0 (  I  ,1  )  = ( C (  I  , 4 ) -C {  1,2)  ) / A  1 
B(I,2)=(C(I,l)-C(I.3))/(Al*B{I,l)l 
K  =  2* 1+3 
C(I ,  I  0 )=b( I,1)*(F(K,7)+(B(I,2)-B(3-1,3) ) 

*FUNC(F (K, 7) ,F( K+l , 7) ,XNF( I ) ) ) 
CONT INUF 

IF{ ABS( ( XNF( 1 ) -0( 1 , 1 )) /XNF { 1 )) .LTi 0.0001 )     GOTO    20 
00    3  0     1=1,2 
X\IF{  I  )  =b(  I  ,  1  ) 

0(  I  ,  3  )  =0(  I  ,  1  >  *  (C(   I,5)+C(  I  ,  5)  *(B(  I  ,  2)-B(  3-1,3))) 
CONT I NUE 
GOTO     3  1 
IP0S=1 . 5*N0FRM+1 .25 

CONVFRT     DATA     TO     CYLINDRICAL     COORDS. 


1  7 


DO     15     I  =1  ,  NUMSG1  ,2 

DO      16     J  =  1  , 2 

DO     17     K=l ,4 ,2 

CALL     30TAX(D(K , l +J/2)   ,D(K+l.I+J/2), 

-CKKfl  )/2,  7  )  ,A(  K,J)  ,A(K+1,J)  ) 
A(K,J)=A(K,J)-XAVG(K) 
A(K  +  1  , J)=A( K+l  ,J)-XAVG(K+1  ) 
DUM1 -XNF ( 1 ) *XNF( 1 ) 
DUM2=XNF( 2) *XNF( 2 ) 
0UM1  =  SQST(1  .  19* (  \(  I  » J  )*A(1  ,J)+A(2.J)+A(2,J))*DUM1+16) 
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DUM2 

DUM3 

A(  2, 

no  l 

i  e 

A  (  i<  , 

A(  1  , 

A(  1  , 

1  6 

A(  3, 

M  =  ^  + 

no  l 

X(M1 

X(M1 

X(  MI 

I  F  (  M 

O.J  2 

^  6 

V  (  M  1 

GOTO 

2  7 

00  2 

V(M1 

IF(  L 

IF(A 

9     V 

V  (  M  1 

23 

IF(L 

10  AD 

IF(  I 

24 

VZAV 

VZ4.V 

2  2 

CONT 

1  5 

CO  NT 

=SQRT(  1  .  19*(  A(3  ♦  J) *A( 3  ,  J)+A(4 , J )* A(4  f  J)  )*DUM2+16) 

=A{4»J)*B{2, 1 )  +  ( 1+B{2, 2 )/0UM2 } 

J )  =- A (  2  ,  J ) *  3 (  1  ,  1 ) r { 1  +  (  0(1  , 2 ) -DU  M3 ) / DU  M 1  ) 

/  (  1«-A(2,  J  )  *A<  4,  J  )*3  (  1  ,  I  )  *D(  2  ,  1  )  /  (DUM1*DUM2)  ) 
8  K  =  3  i  4 

J)=A(K,  J)*B< 2, 1) M  1+(D(2»2»-A(2 .J ) ) /DUM2 ) 
J)=A(1,J)*0(1,1)*(1+(0(1,2)-A(+,J))  /DUM  1  ) 
J)=IPOS-A(  I  ,  J)+C(  1  dOJ 
J  )  =  I  PO  S-  A  (  J  ,  J  )  +•  C  (  2  ,  1  0  ) 
1 

5  J  =  l  ,2 

, J ,  1  ,M)  =  (  A(  1  , J )+A(  3,  J )  J/2 

»  J  .  2  ,  M  )  =S  Q  IT  (  A(  2  ,  J  )  *  A  (  2  ,  J  )  +A  (  4  ,  J  )  *  A  {  4  ,  J  )  ) 
,  J  ,  3  ,  M ) =  A  T  AN2 ( A ( 4  .  J )  ,  A ( 2  ,  J  )  } 

. N  E .   1)  GOTO  27 
5  L=l  ,3 
•  J  .  L  .  M  )  =  0 

15 
2  L=l ,3 
. J . L . M ) = ( X ( M 1 , J . L , M ) -X ( M 1 ,  J  ,  L , M- 1 ) ) /TO 

.  NE.  3)  GOTO  2.3 
OS(TO  +  V(M1,J,3»M))   . G  T .  5.5)  V(M1,J,3,M)  = 
(  Ml  ,  J  ,3,M)  -S  I  GN(6  .2331  32 /TO  ,  V  (M  I  ,  J  ,  3,  M  )  ) 
t  J  «  L  .  M  )  =  V(  Ml  f  J  .  L  ♦  M)  *  X  (  Ml  ,  J  ,  2  ,  M) 

.Nt.   2)  GOTO  2  2 
=  X  (  wi  ,J,2«M)/0.1+1 
"-.  AO  .GT.  6)   I  PAD  =5 

G  (  1  ,  I  R  A  D  )  =  V  Z  A  VG  {  1  .  IP  AD  )  +  V  {  M  1  ,  J  ,  1  ,  M  ) 
G(2,IRAD)=VZAVG(2tIRAD)+l 
I  NOP 
I  NUE 


c 
c 

•*  *  * 

LIST  DATUM  POINT  DATA 

WR  I  TE  (  6  ,  1  20  )  ND=«M  .NOTRK  .NUMRUN 

120 

FC ?M AT (  •  1  '  ,  ' FPAMc  •  ,  I  2  « 1  OX « •  TRACK »  , 13  » 1  OX  ,  'P UN •  *  I  3 , 
•          ///.  1  X,  '  F,~  AMR  DATUM  POINTS   :',/,  1  K  ,  5  {»=  =  ==  i  )  , 

///,  1  1  X  ,  'UPPER  VIEW           L0.VE3  VIE*'./. 

2X,  'PCS'   ,  2X  ,2(2X,  •=  =  =  =  =  ==  =  ==•==  =  =  =  =•  ),/, 

2X ,  • NUM          / 1          X           Z2          Y  •  ,/ , 

IX ,  •         «  .  4(  2X,  •          •  )  ) 
Wf>  I  TE  (  6  ,  1  02  )   (I,{F{J,IJ,Js=l,4),I=l,7) 

c 

c 

102 

F  0  3  M  A  T  (  1  X  ,  I  3  ,  2  X  ,  4  P  9  .  0  ) 

**  * 

LIST  SEGMENT  END  POINT  DATA  FHOV     DATA  COLLECTION, 

r 

CYLINDRICAL  COORDINATE  DATA,  AND  VELOCITIES  COMPUTED 

c 

c 

BY  DIFFERENCE. 

WRITE ( 6.121) 

12  1 

FO  .:MAT(////  ,  1  X  ,'  SEGMENT  END  POINTS    :•,/,  1  X  ,  5  (•==  =  =•)  , 

///,  1  X  ,  •  TI  M»  ,6X  , 

•TIM  ',  CX,  'UPPER  VIE. 7           LOWER  VIEW. 

13X. • CYLINDRICAL    COORD  I  MATES '  ,  20 X, 
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•VELOCITIES.  DI A/SEC »,/, IX , 
•INT        Ol  D2         Do 

1 3X  ,  •  Z . DI A  * i 5X  » • R. D I  A • . 5X . • THETA • 
•  Ml*   ,  12X,  • VR •  ,  1  OX,  ■ VT  HETA ■  . 

/,1X,« ,.4(« '),5X 

j(t '),4X,3(*     

L=M-< NUMSG1-1 ) /2 
DO     19     J  =  l  .NUMSG1  ,  2 
0  0     i  9     K  =  1  , 2 
Kl=J+K/2+2ML-l) 
19    i«oiTE(6.101  J     Kl . ( D( I. J*K/2) , 1=1 .4) , 
(  X(  Ml  fKf  I  «L+J/2)  .1=1  .3) 
,  ( V( Ml  , K  ,  I , L+J/2  )  ,  I  =1  .  J) 
.01     FORMAT C IX  i  I3,lX,AF9.0.dX»3F10.3.6X.lP3E14.3) 
IF     (NOFPM     . LT.     NUMFPM)     GOTC     2 
N(  Ml  )=M 
Ml  ="11  +  1 

I1-"     (Ml      .LE.     7)     GU     TO     1 
M  1  =  M  1  -  1 
0  0     2  5     1=1,5 

IF(VZAVG(2.I)      .EQ.     0)     GOTO     25 
VZAVC  (1,1  )=VZAVG(  1  ,  I  )/VZAVG(  2,1) 
25    CONTINUE 


D4 
•  14X. 


)  ) 


SAVE     CYLINOPI 
TUN  5 


AL     COORDINATE     DATA     F0,~<    PARTICLE    P05I- 


WR I TE ( 2  )     Ml  ,N  ,  VZAVG 

DO    6  0     I = 1 . M 1 

DO     6  0     K=l ,3 

WPITE(2)      ( ( X( I ,J,K,L) ,J=1 ,2) ,L=1 ,21 0) 

siK^ITE  (2)      (  (  V(  I  ,  J.K.L)  ,  J=l  ,  2)  ,L  =  1  .21  0) 

CUNT  I NUE 

END     KILE     2 

ST  T^> 

END 
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c_ 
c] 
cl 
c| 


THIS    FUNCTION     ROUTINE     COMPUTES     THE    LAMODA     FUNCTION     !JF  | 
TRXNSFOHMATICN     EQUATIONS  | 

I 

FUNCTION     FUNC(X,Y,D) 

FUNC=X/SQ ?T  (  1  .  1 9  *  (  X*X+ Y+Y) *D *D+ 1 6 J 

RETU.'  N 

END 


C_ 

cl 
c  I 


THIS     FUNCTION     ADJUSTS     THE     ANGLE 

FUNCT ION     ANGLE ( A ) 
ANGLE=A-S [GNC6.2B3186. A) 

RET  J- N 

END 


C_ 

cl 

cl 


SUdROUTINE     TO     ROTATE    COORD     AXIS 


SUJr-OUT  INE     ROT  AX  <  X,Y  ,  A  ,XNEW,  YNErf  ) 

A1=ATAN2( Y.X)  4- A 

R1=SQRT(X*X+Y*YJ 

XNEW=RI *COS ( A  I  ) 

YNEW=R 1*SIN(  Al  ) 

RETURN 

END 


APPENDIX  D 

COMPUTER  PROGRAM  FOR  THE 
COMPUTATION  OF  THE  LAGRANGIAN  TIME  CORRELATIOI 


The  theory  of  this  program  is  found  in  Section 
3.2.2.   Documentation  is  also  included  within  the  program 
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THIS  PROGRAM  COMMUTES  THE  LAGRANGIAN  CORRELATIONS  ANJ 
TU^ajLENT  INTENSITIES  FROM  PARTICLE  DATA  KtLATING 
POSITION  OF  PARTICLE  IN  TIME.   CORRELATIONS  FOR  DOTH 
POSITIVE  AND  NEGATIVE  TIMES  ARE  COMPUTED  FO»  F I Vfc 
ANNULAH  ZONES  OF  A  ONE  INCH  DIAMETER  PIPE.    DEFINI- 
TIONS OF  THE  MAJOR  VARIABLES  IN  THE  PROGRAM  FOLLOrf: 


COR"EL (I  .  J  .K,L  ) 


NCOP( I ,  J  ,K) 


X(  I  , J ,K) 


V( I , J ,K> 

Tl 

T2 

TO 

Ml 
NCI) 


Z  ST  A  R 


*  *  * 
■**  * 


*  *  * 


THE  J-TH  CORRELATION  FD?  THE 
L-TH  TIME   INTERVAL  OF  THE  K-TH 
"ADIAL  INTERVAL.    1  =  1  FOR 
NEGATIVE  TIMES,  AND  1=2  FOR 
POSITIVE  TIMES. 

NUMQER  OF  CONTIBUTIONS  TO  J-TH 
COORELATICN  AT  TIME  INTERVAL. 
K.    1=1  FOR  NEGATIVE  TIME,   AND 
1=2  FOR  POSITIVE  TIME. 
THE'  J-TH  COORDINATE  OF  THE 
PARTICLE  POSITION  FOR  TIME 
INTERVAL  K.    1=2  FOR  RAW  DATA, 
AND  1=1  FOR  SMOOTHED  DATA. 
J=l  FOR  Z  COORDINATE  J=2  F3^ 
R  COORDINATE.  J  =  3  FOR  THETA 
COORDINATE- 

THE  J-TH  VELOCITY  COMPONENT  OF 
PARTICLE  AT  TIME  INTERVAL  K. 
MATRIX  FOR  SMOOTHING  PARTICLE 
POSITION  COORDINATES. 
MATTIX  FOR  SMOOTHING  VELOCITY 
COMPONENTS . 

TIME  INTERVAL  OF  PARTICLE 
POSITIONS,  SECONDS. 
NUM3ER  OF  PARTICLE  TRACKS. 
NUMBER  OF  INTERVALS  FOR  PARTI- 
CLE TRACK  I  . 
MEAN  Z  POSITION  OF  PARTICLE. 


THE  CORRELATIONS  ARE  NUMBERED  AS  FOLLO*: 


i:  <Z»(T)*Z,CT)>  2 ; 

3  :  <Z«  l   T)  «THETA  •  (  T)  >  41 

5:  <F<«  (T)*THETA«  (T}>  6: 

7  :  <^(T)*v^(T)>  a: 

9:   <Z« ( T}*VTHETA • ( T ) >      10: 

l  l  :  <  '<  •  (  t  )  *  v  R  (  T  )  >  12: 

13:  < THETA' { T)*VZ»<T>>      14 

IS:  <THETA» ( T) *VTHETA« ( T )>16 

17:  <VZ* (TO) *VR(T)>        la: 

19:  <VR(T0)*VZ« (T) >         20 

Zi:      <VP (TO )* VTHCTA' ( T )>     22: 

23:     <VTHETA« { T0)*VR< T)> 

24:  <VTHETA' (T0)>*<VTHETA« (T)> 
«ERE  Z'       =  Z  DISPLACEMENT  FROM  ML  AN  PATH 

R«       =  3  DISPLACEMENT  FROM  INITIAL  ORIGIN 
THETA*   =  (INITIAL  kAUIUS)  MTHETA  DISPLACEMENT) 


<Z' IT ) *R« (T) > 

< R •  (  T  )  *  R  •  (T)> 

<THETA« (T)*THETA« (T)> 

<Z'  (T) *VR (T) > 

<R« (T) * VZ« ( T) > 

<R'  (T  )*VTHETA ' (  T)  > 

<THETA»  (T  )*Vt.(T  )> 

<VZ • ( TO) *VZ« ( T)> 

<VZ»(TO)*VTHETA'(T)> 

<  V:.  (TO)  *  VTHir  TA'(T)> 

<  VT  HET A ■ ( T  0 ) *  VZ • ( T ) > 
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c| 
c| 

cl 

cl 


VZ1      =  DEVIATION  OF  Z  VELOCITY  COMPONENT 
VR       =  R  VELOCITY  COMPONENT 

VTHETA»  =  (INITIAL  RADIUS) *{ T I  ME  KATE  OF  CHANGE 

THETA  DISPLACEMENT) 


DIMENSION  Tl  {  2  1  ,21  )  , T2{21,21  )  .  :?INT(  10)  ,N(  1  J) 
COMMON  X(2,3,420) ,  V  (  2  ,  3  ,  42  0  )  ,  CORREL  (  2  ,  2  4  ,  5  ,2  1  0  ) 

,ARG<9)  ,  NCO"^  (  2  .5  .21  0)  .  IP.AD.ZSTAR 
DATA  FINT/0.0.2fc0.2,2* 0.4 t 2*0. 6 .2*0 .8 i 1 . 0/ 

DO  10  I  =1  ,2 
HO  10  J=l ,5 
DO  10  K  =  1  , 2  1  0 
NC  3R  (  I  .  J  .  K  )  =0 
DC  10  L=l .24 
CORREL{  I  iLiJ,K)-0 
CONTINUE 

n  p  t = a 

T0  =  0. OOol 6 


1  0 


c 

c 


1  01 

c 

c  **•* 

c 

c 


16 

1  7 


READ  IN  SMOOTHING  MATRIX 

READ (5» 101)  (( Tl( I,J>  ,1=1 , NPT) ,J=1 , NPT) 
READ (5. 101  J  (( T2(  I , J) » 1=1 . NPT) »J=1 ,N3T) 
FORMAT I 5F15.9) 

READ  IN  NUMBER  OF  TRACKS.  Ml.  AND  NUMBER  OF  TIME 
INTERVALS  IN  EACH  PARTICLE  TRACK 

r-"E  *D  (  3  )  M  1  ,  N  .  (  DUMMY  ,1  =  1,10) 

PI2=6. 283186 

DO  15  I  5=1, Ml 

N1=2*N{ 15) 

DO  17  K=l ,3 

^FAD(3)   ( X(2,K,L) ,L=l . 420 ) 

»EAD(3)   (  V(2,K,l_)  ,L=1  .420) 

IF  (  K  .NE.  3)  GOTO  17 
DO  16  L  = 1 . N I 

IF  (X{2,K,L)   .LT.  -2.7)  X  (  2  .  K  ,  L  )  =P  I  24-  X  {  2  ,  K  ,  L  ) 
CQNT INUE 
CONT INUE 

SMOOTH  PARTICLE  POSITION  DATA  AND  COMPUTE  VELOCITIES 
FROM  SMOOTHED  DATA 

DO  3  1  K  =  l  ,3 

DO  31  L=l ,N1 

INTX-XI 2.K ,L) 

X  0  =  I  N  T  X 

V0  =  0 

M2-2«(L/2)-A 

M3=L-2 *(L/2)+4 

IF((L  .GT.  3)   .AND.  (L  .LT.  Nl-2))  ^DTD  22 

M2  =  0 
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?  o 


:n 


M3-L 

I  F ( L   .  L  T .  4 )  GOTO  2  2 

M2  -Nl  -'3 

M3-L- M2 

DO  20  M=l . NPT 

X0  =  X0-K  X  (  2  ,  K  ,  M  2+  M  >  -  I  NT  X  >  -  T  1  (  M  J 

M) 

V  0  -  V  J  «-  X  (  2  ,  K  ,  M  2  +  M  )  *  T  2  (  -1  J  i  M  ) 

CO NT INUE 

X(  1  ,K »L)=X0 

V(  1  t KiL) -VO/TO 

CONTINUE 

COMPUTE  THE  CORRELATIONS 

DO     3  5     K,  =  1,N1 
N2  =  K  + 1 
N3-N1 

I  •?  A  D  =  X  (   1, 2, K)/0. 1+1 
IF     ( IRAD     .GT.     5)      IhAD=5 
ARG  (  7)  =V(   1  ,  1  ,K)  -  VEL(  X(  1  .2  ,  K)  ) 
ARG(8) =V( 1 .2 ,  K ) 
\R  G  (  9  )  =  V  {  1  ,  3  ,  K  )  «  X  (  1  ,2  ,  K  ) 
ZSTA-  =X(  1  .1  ,  K) 


FO  ?     NEGATIVE 


I  Me 


2  5 


OH     2-3     M-l  t  K 

M  2  =  K  -  M  +  1 

IF(K-M2     .GT.     2  09)      GOTO     2d 

CALL     COPLAT (K ,M2 ) 

ZSTA-  =  ZSTAR~TO*VEL(X{  1   ,  2  ,  M  2  )  )/2 

CONT INUE 


C  **  * 
C 


2  9 

j{: 

1  5 


37 


FO?  POSITIVE  TIMES 

ZSTAR=X ( 1 ,1 ,K) 

IF  (N2  .GT.  N3)  GOTO  3  6 

DO  2-i     M  =  N 2  ♦  N 3 

IF(M-K   .GT.  2  09)  GOTO  36 

ZSTA-  -I  STAR  +T0*VEL(  X(  1  ,2.MJ )/2 

CALL  CQPLAT(K.M) 

CONTINUE 

CONT I NUE 

CENT  I NUE 

DO  3  9  L  = 1  t  2 

DO  3  B  J  =  1  ,  5 

DO  3  3  1=1,210 

IF  ( NC OH ( L.J.I ) 

DO  3  7  K=l ,24 

C03RFLI L,K. J , I ) 

CONT INUE 


•EO.  0)  GOTO  38 

CORBEL  (L.K  »  J,  I  )/NCORIL.  J.  I  ) 
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***  SAVE  TURBULENT  intensities 

WRITE  (  1  )  (  (  CURREL  (l.IiJil)  iJ  =  l  .C)),I-16,«) 


LI5T  TURBULENT  I  NT  ENS  I  T  I  ES 


'.*"ITE  (6,  142) 
142  FORMAT ( // , 1 X, • TU3DULE JT  I MTENS I T I ES • ) 

A)RITE(6  ,1  35)  (  :UNT(  I  )  ,  1  =  1  ,  10) 
135  FORMAT (//////  ,4BX ,' RAO  IAL   I  NT El V AL •  , / .  3 X  ,  •       ',OX, 
39  {  •  ='  )  ,  /  ,  3X  ,  '   NUM  •  »2X,5(F9.  1  ,  •<R/RO<«  ,  F3.1  ) 
/,  Uf9(  •  =  •)  .2X,  5(  IX  ,'  ==  =  =  =  =  =  =  =  =  =  =====  •)) 

00  3  a  1=10,24 
3  4  WRITE  (6  il  32)  I  .  (  :  Of-  RE  L  (  1  ,  I  ♦  J  ,  I  )  ,  N  COP  (  1  ,  J  ,  1  )  ,  J=  1  , 5  ) 
13  2  F  0  EM  A  T  (  I  X  ,  I  S  ,  1  X  ,  5  (  1  P  E  1  3  . 4  ,  •  /  •  ,  I  4  )  ) 

**♦  NORMALIZE:  DOUBLE  VELOCITY  CORRELATIONS 


OH  4  0  L=l , 2 
DO  4  0  J=l , 5 
DO  4  0  1=16,24 
11=1 6+4*1 ( 1-16  1/3) 
Cl=CORREL  (1  ,  I i  ,  J,  1  ) 
C2=C0RREL (  1  ,  I  1 +4v (  1-1  1  )  , j,  1 ) 
IF  ( ( CI   .EG.  0)  .CR.  ( C2  , bu. 
CJ=SCST(Cl)t50"T( C2 ) 
DP  4  1  K  =  1  , 2 1 0 

IF  (CORR£L(L , I  ,  J ,K)   .EQ.  0)  GOTO  4  1 
I  F  (  (  K   .  E  Q  .   1  )   .AND.   (  4  *  {  I  /  4  )   .  E'O  .   I  )  ) 
C  OR  PEL  (  L  .  I  ,  J  »  K  )  =C  0  :~  RE  i_  (  L  ,  I  ,  J  ,  K  )  /  C  3 
4  1  CONTINUE 

4  0  CONTINUE 

00  5  0  L  =  l  ,2 

DO  5  1   J  =  1  ,  5 

DO  5  1  M  =  1  6  ,  2  4  ,  '+ 

C0RREL( L,  4. J, I ) =1  .0 

5  1  CONTINUE 


0)  )  GOTO  40 


GOTO  41 


■*    LIST  CORRELATION! 


DO 
IC  = 

00 
c  1  = 
DO 
CI  = 
I  F 
IF 
WRI 
FO  'i 


>  0     J  =  1  ,  2  4 

■2      M  =  1  ,  2  1  0 

I 

K  =  l 


)  )      GOTO     53 


-i     K  =  I  ,  5 

1 +AOS  f CORREL (L , J , K  ,M )  ) 

(CI      .EQ.     0)      .AND.      (M     » NE.      1 

I  C     .  L E .     7  0)     GOTO     53 

E  (  6  ,  1  3  0  )     J  ,  (  R  I  N  T  (  I  )  ,  I  =  1  ,  I  0  ) 

AT(  'lS'CO^fcLATIUN'  ,1  J,'      /     NUM3tiR'    ,///, 

4faX,  •  RADIAL      INTERVAL'  ,/  ,  3X,  •  T I  ME •  , 6X , 09 (  •  =  ♦  ) . /  , 
3X,  • MSEC'   , 2X  ,5( F9. 1   ,  ■ <R/RO<«  ,F3.  1  )  ,/ ,  1 X  , 
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5  3 

52 
131 

c.  0 

r 


9(  •=•)  »2X,5(  IX,  •==========  ===  =  =     •)) 

IC  =  1 

A2 -5  00  ~  TO  *(  M-l  ) 

IC=IC+  1 

WP  ITt.  (  5,  1  3  1  )     A2.  (CORREHL,  J»K,M)  ,NCOP .  (L  .  K  ,  M)  , K=  1  ,  5  ) 

F  n  i  M  ,\  T  (  lX,F8.2,lX,b(lPE13.4,,/«  ,14)  ) 

CO NT  I NUE 

SAVE     COF DELATIONS     AND     NUMBER    OF     CONTRIBUTIONS     TO     EACH 
COR^ELATI ON 


DO     5 1      I  -1  ,24 
DC     6  1      J  =  l  ,  5 

61  wRITE(l)      ( (CORREL <L, I  » J,K)  ,K  =  1  ,210)  ,u  =  l  ,2) 

DO     62     J-l , 5 

62  iVRITE(l)      (  (NCO  P(  L  ,  J  ,K  )  ,K-1   ,2  10  )  .L  =  l  ,2  ) 

EN  J     FILE1     1 

STOP 

END 
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C|  SUBROUTINE  FOR  COMPUTING  C ONTR I  OUT  I  ON S  TO  CU3RELA-  | 
C|  TUNS.  PARAMETER  K  IS  THE  NUMBER  OF  THE  STARTING  TIME| 
C|  INTERVAL  AND  PARAMETER  M  IS  THE  NUMBER  OF  THE-.  ELAPSEoj 
C\                TIME  INTERVAL.  | 

C  1 I 

SUBROUTINE  CORLAT (K ,M ) 

COMMON  X(  2,3.42  0)  ,V{  2  ,  J  .420  )  ,CORREL<2  •  2  4.  5  .210) 
.  .  APGC'3)  ,  NCOR  (  2.5.  21  0)  .  IP  AD  ,  ZSTA-< 

IT IM= IABS(M-K ) +1 

I =31 GN(  1  . 0 , M-K-0 .  1  ) 


C      V  *: 

C 


CCMPJTE     PAFTICLE     DISPLACEMENTS     AND     VLLDCITY     DEVIA 
TIONS 


ARG  (  I  )  =  I*  (  X  (  1  ,  I  ,  M  J-ZST  AP. ) 
ARG ( 2 ) = I * ( X(  1  »  2  ,  M  )  -X ( 1   . 2 , K )  > 
ARG(3)=I*X(1  ,2,K)  M  X  (  1  ,  3  ,  M  )  -  X  {  1  ,  3  ,  <  )  ) 

ARG (4 )=V( 1 .1 ,1)~VtL(X( 1 , 2, M) > 

A  P.  G  (  5  )  =  V  (  1  ,  2  ,  M  ) 

A  R  G  (  6  )  -  V  (  1  ,  3  ,  M  )  *  X  {  1  .  2  .  K  ) 

I  =  l  +  (  1  +  I  )  /2 

NCU«<  1  .  I  RAD,  IT  IM)=NCOR  <  I .  I  FAD,  IT  IM)+1 

I  .3  =  0 

no    36     11=1,3 

DO    3  6     12=1,3 

IF     (12     .LT.      II)     GOTO     3b 

13=13+1 

***     COMPUTE     CONTRIBUTION     TO    CORRELATIONS     1      THROUGH    6 


35 
C 

C  *** 

c 


c 

c  *** 

c 


CORRELC  I  .  13  ,  I  RAD,  IT  IM  )  =CGRREL  (  I  .  13,  HAD,  I  TIM  ) 

+  ARG(  I  1  )  tA%(  12  ) 
14  =  3*1  1+12 

COMPUTE     CONTRIBUTION     TO    CORRELATIONS    7     THROUGH     15 

CO"'EL(  I  ,  14+3,  IS  AD.ITI  M)=CCRREL(  I  ,14+3,  IRAQ,  IT  IM) 

+  A!!G(  I  1  )  *Afi  G(   12  +  3  ) 

COMPUTE  CONTRIBUTION  TO  CORRELATIONS  16  THROUGH  24 

CORRELC  I  ,  14  +  12 ,  Im AD, IT IM) =CORREL(  I.I4  +  12,IRAD,ITIM) 

+AR6( 12+3 ) *ARG( 11+6) 
SETU-^N 

END 
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cl 
cl 
c| 


**     THIS    FUNCTION    COMPUTES     THE     MEAN     AXIAL     VELOCITY     US  IN; 
AN    EQUATION     FORM     SUGGESTED     dY     PAL.         THE     FCRM     IS 

UZ/UZtMAX    =     Al* ( R/RO )**2+A2* (R/RC )**<2*MJ 

FUNCTION     VEL(R) 
'< 1 =  2 «R 

VEL  = 1 -0 . J  62  5  e  *R 1  *  ft  I 

IF( ALQG L0(R1)      .LT.     -0.1)     GOTO     10 
VEL  =  VEL-0. 637421^  Jt*j SO 
10     VEL=47.1435*VEL 
■ETU-N 
END 
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